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1 Introduction 



The purpose of these notes is to review one of the branches of modern string 
theory: the theory of matrix models with the emphasize on their intrinsic 
integrable structure. We begin with a brief description of the field and its 
place in the closest environement within entire string theory. 

The main content of the string theory Q is the study of symmetries in 
the broadest possible sense of the word by methods of the quantum field 
theory. The usual scheme is to start from some symmetry and construct a 
field-theoretical model (usually 2-dimensional, for the reason that we do not 
discuss here), which possesses this symmetry in some simple sense (e.g. as 
Noether symmetry or as a chiral algebra). The main idea at this stage is to 
find a model which is exactly solvable (if nothing but the symmetry is given 
this is a nice principle to restrict dynamics). The next step is to study the 
hidden symmetries of the model, which are somehow responsible for its exact 
solvability and are usually much larger than the original symmetry. 

This "inverse" step: model — > symmetry can be made at least with three 
different ideas in mind. 

One can look for some hidden local (gauge) symmetry of the model, which 
is fixed or spontaneously broken, i.e. identify it with some other model which 
had more fields - auxiliary from the point of view of the smaller model and 
gauge - from that of the larger one. (Examples: gauged Wess-Zumino-No- 
vikov- Witten (WZNW) model, topological theories in BRST formalism etc.) 

One can take for a new (full) symmetry of the model just its operator 
algebra (algebra of observables) (see 0, || and also M for the first results 
in this direction). It deserves mentioning that gauging of entire algebra 
of observables gives rise to a "string field theory", associated with original 
model (considered clS cL string model). 

One can construct effective action of the theory by exact evaluation of 
the functional integral. 

As to direct step symmetry — > model, one can take as example the best 
understood case, when original symmetry is just a Lie algebra. Then the 
quantum mechanical model can be constructed by geometrical quantization 
technique (see || for the most important example of Kac- Moody algebra and 
the WZNW model). 

2 See jj] for a general review 
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From the mathematical point of view the two elements of the above 
scheme look like algebra (theory of symmetries) and Analysis & Geometry 
(field-theoretical models). The idea of constructing models with a given sym- 
metry (and nothing else relevant for the dynamics) can be identified with the 
mathematical concept of "universal objects". 

| SYMMETRY | = | ALGEBRA | 

T I 

| MODEL | = | ANALYSIS | 

with this symmetry \ | GEOMETRY 

Theory of everything 

The sequence of iterations of the two arrows in the picture leads to 
a deaper understanding, enlargening and generalization of all the notions 
involved: symmetry, exact solvability, field theory, geometrical structures, 
quantization etc, thus stimulating considerable progress both in physics and 
mathematics. If this iterative process can somehow converge, the limit point 
will deserve the name of the theory of everything, which will indeed unify all 
the possible field theoretical models by embedding them into a huge, but well 
structuired theory, which will be also exactly solvable in some, yet unspecified 
sense of the word. We refer to |T| for more details about this semiphilosoph- 
ical programm, known under the name of (modern) string theory, and now 
turn to a more narrow subject: the theory of matrix models. 

At the moment it is mainly associated with the theory of effective ac- 
tions, at least so far this is where the main results of the modern theory of 
matrix models find their applications. This technique is especially suited for 
the study of effective actions, obtained after integration over 2-dimensional 
geometries (including the sum over genera) and it produces non-perturbative 
(exact) partition functions of particular string models. The main result of 
these studies points out that these partition functions exhibit two remarkable 
(though expected |J) properties: 

First, effective action for a given model is essentially the same as for 
any other model. In fact effective action is a function of coupling constants 
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("sources" in the old-fashioned terms), which are nothing but coordinates 
in the space of various models (configuration space of entire string theory): 
variation of couplings change one model for another. 

Second, effective action possesses a huge additional symmetry, which is 
somewhat similar to the general covariance in the space of all models (the 
above mentioned configuration space) and in the simplest examples, which 
have been studied so far, can be expressed in terms of integrable hierarchies. 
(This "general covariance" in the configuration space can after all turn into 
the main dynamical principle of the string theory.) 

Both these features seem to be very general, arising whenever the largest 
possible Lagrangian with a given symmetry is considered (without restric- 
tions on the possible counterterms, imposed by requirements of renormaliz- 
ability or by locality- minimality "principles" - this is why this phenomenon 
is not widely known to field theorists). An example of highly-nontrivial cal- 
culations leading to similar conclusions can be found in ref.0. 

We hope that these remarks will become more clearer after some specific 
examples will be considered below. Still they deserve being fromulated in 
the full generality, not only to intrigue the reader, but also because they can 
serve for better understanding of the ideas and outcomes of generic string 
theory. 

The " corner" of the string theory, associated with matrix models, can be 
described by the following picture: see Fig.l. 

The big blocks within the body of string theory, which are directly re- 
lated to matrix models are: theory of conformal models, that of the N = 2 
supersymmetry and the (loop-equation version of) the Yang-Mills theory (in 
any dimension). Also the Einstein gravity should be related to the subject 
in a way, similar to Yang-Mills theory, but these links are yet not clarified. 

Both conformal theory and N = 2 supersymmetry are sources of the con- 



cept of "topological models" ||-[jnj]. These arise after gauging of all contin- 
uous symmetries of the WZNW models and/or as models with BRST-exact 
stress tensors, naturally appearing in the context of iV = 2 supersymmetry. 
If formulated in a self-consisted way in the "universal module space" (uni- 
fication of module spaces of all finite-genus Riemann surfaces and bundles 
over them) these models turn into those of "topological gravity" . Generating 
functionals of topological gravity models in fact generate infinite sequences 
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of topological invariants of certain spaces (inverse definition is also possible 
in some cases [§], though the universal (generic) algorithm for the operation 
topology of some space — > topological gravity is not yet formulated). 

Alternative models of 2c? quantum gravity arise straightforwardly from 
conformal models through a procedure of "summation over geometries". 
There are two essentially different approaches to the problem. One ("Polyakov 
approach") is to make use of the complex structure, intrinsic for conformal 
theory |12| and sum over Riemann surfaces, what involves integration over 



module spaces and sum over genera. The main technique used in this ap- 
proach is the theory of free fields on Riemann surfaces [O], |TJ| and bosoniza- 



tion formalism for conformal field theories |L5|, [T^|. This approach requires 



solution of Liouville theory, which still remains a problem under intensive 
investigation (in turn related to conformal field theory). Further progress in 
this direction should be related (or can be expressed in terms of) the adequate 
theory of the universal module space, handle- gluing operators etc. Similar ob- 



jects arise in the field-theoretical approach to the topological gravity (see |L7 
for a recent review). 

Alternative approach to summation over geometries does not refer at all 
to the complex structure and instead involves a sum over random equilateral 
triangulations 118|-||20||. Q This is the place where matrix models first appear 



in the context of string theory. The random triangulation approach is by no 
means specific for conformal models (since it ignores the complex structure) 
and can be applied in many other situations - for example, to Yang-Mills 
(YM) theories in any dimensions (where instead of summation over geome- 
tries one needs "simply" to sum over ordinary Feynman diagramms). 

Applications of the matrix-model method usually involve two steps: for- 
mulation and the study of "discrete" model and then taking its "continuum 
limit", giving rise to a new - "continuous matrix model", which sometimes 
can again be represented in a form of some matrix integral. 

One of the main discoveries in the field of matrix models is that continuous 



3 Its relation to the Polyakov approach is a separate very interesting, important and 
badly understood problem, which allows a non-trivial reformulation in terms of num- 
ber theory (see pi]). The main puzzle here is that equilateral triangulations are in fact 
arithmetic Riemann surfaces - a dense discrete subset in the entire module space, with 
interesting and deap algebraic properties. Equivalence of the two approaches to 2d quan- 
tum gravity should imply the existence of some number-theoretical background behind 
the scene, which would be very nice to discover in full purity. 
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models arising finally from the random-equilateral-triangulation description 
of the simplest (minimal with c < 1) string models coincide with the sim- 
plest (CP 1 Landau-Ginzburg) models of topological gravity P,P^|-P^|: two 



(classes of) theories are identical (this is not yet proved in full detail, but is 
more than plausible). 

So far continuous models are actually found and somehow understood 
only for string models, based on the c < 1 minimal conformal theories (more- 
over, only for q = 1 in the (p, g)-seria). Conformal models with c > 1, which 
are relevant for description of gauge theories in space-time dimension d > 2 
(which possess particles, rather than only topological degrees of freedom), 
should give rise to the discrete matrix models with "non-factorizable" inte- 
gration over " angular variables" , of which the simplest (solvable) example is 
Kazakov-Migdal model [p5|1 . The issue of continuum limit for such models is 
yet not understood (at least in terms of integrable structures, which should 
probably generalize the familiar theory of Toda hierarchies). 

The goal of the study of matrix models is three-fold. First of all, one can 
look for the non-perturbative (exact) answers for the physical amplitudes in 
the given model. This is the subject which attracts most attention in the 
literature (for several obvious reasons). However, it is equally (and, perhaps, 
even more) important to understand the mathematical structure behind the 
matrix models (which involves topics like general theory of integrable hierar- 
chies, geometrical quanization, Duistermaat-Heckman theorem ("localization 
theory") etc). Also important for the purposes of string theory is to use the 
results of the study of matrix models in order to unify a priori different mod- 
els (according to the above mentioned principle: non-perturbative partition 
functions for different models differ by a change of variables in the space 
of coupling constants). Matrix models already played an important role in 
making this principle more clear and acceptable for many string theorists. 

Let us make the next step and look even closer at the field of matrix 
models, especially, at its mostly studied domain, associated with the d < 2 
string models. Then the following structure will be seen: Fig.2. 

The sample example of matrix model is that of 1-matrix integral 

Z N {t} = c N [ dHe£*=o tkTvH \ (1.1) 

JjVxJV 
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where the integral is over N x N Hermitian matrix H and dH = Flij dHij. 
There are three directions in which one can proceed starting from (|1.1|) . 

The first one |2B| is to look for an invariant formulation of properties 
of the functional Z^{t). It appears to satisfy the infinite set of differential 
equations (in fact these are just Ward identities for the functional integral 



13) [0): 

L n Z N {t} = 0, n > -1 

9 _i_ V 

dt k+n T ^k=0 dt k dt n . k ' 



which is known under the name of "discrete Virasoro constraints". Z^{t} 
can be represented as a correlator of screening operators in some auxiliary 
conformal model (of one free field on the "spectral surface"), and Virasoro 
constraints ( |1.2| ) are of course related to the Virasoro algebra in that con- 
formal model. Also Z^{t} is some r-function of integrable "Toda-chain" 
hierarchy (in fact this statement should be a corollary of the Virasoro con- 
straints, but this relation is still not very well understood). 

The most straightforward further developement [p6fl , ||28|| is to take the 
continuum limit of the Toda-chain hierarchy. In the specially adjusted ("double- 



scaling") limit p0|) it gives rise to the KdV-hierarchy, and the correspond- 
ing r-function appears subjected to the slightly different constraints ||29||,p8| 
(which again form a Borel subalgebra of some other "continuous Virasoro 
algebra" ) : 

C 2n Z cont {T} = 0, n > -1, 

C-2n = \ Eodd k=l k(Tk + T k ) dT ^ +2n + | Eodd \=l dT k dT 2n - k + 

+ ^n,0 + i( T l+ r l) 2 ^--l' (1-3) 

where r k = —\5kz- In fact 



16 

2 
3 

rcont 



z contr T y _ ^ Z N {t} 

*2fe+l— U 

and T are related to t by linear transformation 



;i.4) 



Tk = r(! + D' *° dd; 

9m = mt 2m , m > 1; g = 2N. (1.5) 
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This z cont {T} can be again represented in the form of a matrix integral (over 
nxn Hermitian matrix) ||22|| , p0| - p3|| : 



Z cont {T} = Z V {T} (1.6) 



with V(X) = 4f , where 



Z V {T} ~ F Vtn {L} = I dXe- tvV{x)+tTLX (1.7) 

Jnxn 

and 

T fc = -trZT*/ 2 , k odd. (1.8) 

The function Zy{T} (but not jFy n {L}) is in fact independent of n: the 
only thing that happens for finite values of n is that the r.h.s. of ( p..7|) can 
not describe 2y{T} at arbitrary points in the T-space, in accordance with 
(|1.8|) . Continuous Virasoro constraints ( |1.3|) are in fact equivalent to the 
trivial matrix-valued Ward identity 

(v'(^-)-L]^ n {L} = 0. (1.9) 



dL 



tr 



Another direction to proceed from the discrete 1-matrix model is to 
rewrite it identically in the form of Kontsevich model: this time with V(X) = 
X 2 and additional factor of (det X) N aunder the integral in Ty^ n \L\ [0^1 . 
Then the double scaling limit can be studied in internal terms of Kontsevich 



models 36 



The third direction is towards multimatrix models. In continuous version 
they should provide r-functions of reduced KP-hierarchies |37| (KdV is the 
p = 2 reduction), which are subjected to "continuous H^-constraints" P5| ]. 
Matrix models of such r-functions are Kontsevich models with V(X) ~ X p+1 
At discrete level, however, things are not so simple. The most 



popular discrete multimatrix models |34| are defined as the multiple matrix 
integrals of the form: 



Z N {t^} 

c- 1 / dHU...dH^ TT e^-^ H U ff e ^H^ ( L1 °) 

JNxN 



a=l a=l 
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(the form of the "interaction term" TrH^ a ' is restricted by the "solv- 
ability" principle, but not unambiguously). In fact these models are particu- 
lar examples of the " scalar-product eigenvalue models" and are not really dis- 
tinguished except for the 1-matrix (p = 2) and 2-matrix (p = 3) cases. This 
is reflected in the absense of any reasonable Ward identites and integrable 
structures for these models, which would somehow involve their dependence 
on the variables t^ a ' with 2 < a < p — 2. Therefore the " mult i- scaling con- 
tinuum limit" of these models can hardly be investigated at any degree of 
rigourousness. (It is not so much important for "physical" applications to 
have any discrete models associated with the continuum ones, but this is an 
interesting problem for the "science for science".) For the 2-matrix (jp = 3) 
case the Ward identities can be expressed in the form of " W-constraints" 



[ |38|j and look like |30 



wtt 1] {t}z N {t,t} = n m+n w^{i}z N {t,i} (l.n) 

(here t and t stand for t^ and t^ 2 \ and m, n are any non-negative integers). 

The really interesting set of discrete multimatrix models does exist, but 
it is somewhat different from ( |1.10 ). These theories will be refered to as 



"conformal matrix models", since they arise straightforwardly as general- 
ization of the " CFT-formulation" of the 1-matrix model P5| : it is enough 



to substitute discrete Virasoro constraints in the theory of one free field by 
the Wp-constraints in the theory of p — 1 free fields. Matrix integral for- 
mulation then involves an "interaction term" Det 

instead of e TrH( - a) H(a+1) ^ which is not very easy to guess a priori, but so 
defined models and their continuum limits can be examined in a manner, 
quite parallel to the one-matrix case (though not all is already done in 
this direction). Also this approach opens a possibility to formulate discrete 
models for any set of constraints, e.g. assosiated with the more exotic W- 
algebras and with quantum groups (i.e. can help to solve the inverse problem: 
constraints — > discrete matrix model). This is an option which also 
deserves further investigation. Another natural name for this set of theoris 
is "multi-component eigenvalue models". 

Kontsevich models should be also related to topological models of Landau- 
Ginzburg gravity (LGG), though this relation is not yet clarified in full detail 
(see, however |40|, [p|| ) . 



Among the main unresolved puzzles in this whole field is the description of 
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generic (p, g)-models. Formally, Generalized Kontsevich model ( |1.7| ) provides 
this description, but in fact the partition function (r-function) gets singular 
when the "phase transition" point where the q changes is approached, and 
Kontsevich model with V(X) = polinomial of degree p + 1 provides a nice 
description only of (p, l)-models. Generically, Kontsevich integral describes 
a duality transformation between (p,q) and (q,p) models: (p,q) — > (l,p) 
but not any of these models separately (The only exclusion are (p, 1)- 
models because they are related by Kontsevich transformation to the (l,p) 
models which are completely trivial.) 

In fact continuous models have two different sets of " time- variables" . 
Thus far we introduced T, which are essentially expansion parameters of the 
generating functional for correlation functions. More exact, these parameters 
T depend on the particular model (vacuum), around which the perturbation 
expansion is performed, and they differ slightly from the model-independent 
T. Another set of "times", r k = k ^_ k ^ Kes (V'(fjL)) p dfx, parametrizes the 
shape of the polinomial "potential" V P (X) (of degree p+ 1) and it describes 
the coordinates in the space of (matrix) models. These two types of variables 

- parameters of the generating functional and those labeling the shape of the 
Lagrangian - are almost the same (in fact they would be just the same, if 
there were no loop (quantum) effects). This similarity between T's and r's is 
reflected in the remarkable property of partition function of the (p, 1) model 

- it essentially depends only on the sum of "times" T and r ||40||: 



{T} = f p (r\f k + r k )r p {f k + r k } (1.12) 



with some simple (and explicitly known) function f p . (In eq.( |1.8| ) above for 

X 3 rp rp 1 , 

-g", -tfc — J-k — p 



monomial cubic potential V^{x) — T k = T k = \tiL fc / 2 , while r k 



-14. 
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3- 



The last thing to be mentioned in this general description of the field 
of the matrix model theory is their relation to group theory. Generalized 
Kontsevich model (|1.7|) is intimately connected to the "integrable nature" 
of group chracters and the coadjoint orbit integrals (characters of all the 
irreducible representations of U(N) are usually KP r-functions ||101|| ). In 
fact some "discrete (or quantum) version" of Kontsevich integral is a sum 
over all unitary irreps of U(n) ("integral" over a model of U(n), or over the 
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set of all coadjoint orbits): 



F?{G} = d R XR{G)e~ E *°=° VkCk{R) , (1-13) 



where da, xr and Ck{R) stand for dimension, character and the fc-th Casimir 
of irreducible representation R of U(n). Time variables 7*. ~ \trG , while 
potential V(X) = YlkLo s kX h . This expression can be further generalized to 

J^{G) =Y,XR(G) X R(G)e-^=o^ k (R)_ (i.i4) 

Properties of these "quantum" Kontsevich models deserve further investiga- 
tion (objects like ( |1.13| ) are also known to arise in the localization theory, in 



particular, in the study of the d = 2 YM theory, see, for example, H3) and 



These notes are essentially a review of the views and results of the group, 
working at Moscow (and Kiev). Since references will not be given every time, 
I present here the list of people involved into these investigations: 

L.Chekhov, A.Gerasimov, A.Losev, S.Kharchev, Yu.Makeenko, A.Mar- 
shakov, A.Mikhailov, A.Mironov, A.Orlov, S.Pakuliak, I.Polyubin, A.Zab- 
rodin. 

I also apologize for the somewhat sporadic references to the works of other 
groups. 



2 Ward identities for the simplest matrix mod- 
els 

2.1 Ward identites versus equations of motion 

We begin systematic consideration of matrix models from their simplest and 
at the same time the most basic property: the Ward identites (WI) for 
partition functions. Partition function is by definition a functional of the 
coupling constants in the Lagrangian and WI will be understood here as 
(differentiual or finite-difference) equations, imposed on this functional. If 
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partition function is represented in the form of a matrix integral, f\ the WI 
are usually implied by its invariance (or, better to say - covariance) under 
the change of the integration variables (thus the name "WI"). 

In ordinary field theory we are usually dealing with models, where WI 
either do not exist at all, or at most there is finite number of them - then 
they are interpreted as reflecting the symmetry of the theory. However, by 
no means the finite set of these WI provides a complete description of dy- 
namics of the theory: the number of (quantum) equations of motion (EqM) 
is usually infinite and their solutions are never fixed by the WI. In fact this 
difference between WI and EqM arises because the Lagrangians, considered 
in the ordinary field theory are not of the most general form: they are usu- 
ally severely restricted by "principles" like renormalizability or minimality. 
Because of this there is simply not many enough coupling constants in the 
Lagrangian to describe the result of any variation of integration variables as 
that of the variation of coupling constants, and thus not every equation of 
motion can be represented as (differential) equation for the partition function. 
In other words, by restricting the shape of Lagrangian for "non-symmetric" 
reasons one breaks the original huge "symmetry" (covariance) of the model, 
which was enough to describe all the dynamics (all EqM) as dictated by 
symmetry- and a broader view is necessary in order to recognize EqM as the 
WI, associated with that original high symmetry. This symmetry (it is not of 
Noether type, of course) is peculiar property of all the quantum mechanical 
partition functions, since these usually arise from the procedure of functional 
integration. 

It happened so that matrix models appeared to be the first class of 
quantum-mechanical systems (functional integrals) for which this identity: 

all EqM = all WI 

To avoid confusion we emphasize that such representation does not need to exist, at 
least in any simple form. The more the theory of matrix models develops, the less it has 
to do with matrices and matrix integrals. However (as in the case of entire string theory) 
original name has a tendency to survive. Anyhow, the main content of the theory of matrix 
models (at least of its branch, analyzed in these notes) is the search for invariant formu- 
lations of the properties of partition functions, while matrix integrals (if at all existing) 
are considered as their particular realizations (representations). Moreover, there can exist 
very different matrix integral representations of the same partition function, the simplest 
example being just the basic discrete 1-matrix model, which can be also represented in 
the form of Kontsevich integral (see below). 
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was not simply observed as a curious phenomenon, but became a subject 
of intensive investigation and is identified as the source of exact solvability 
(integrability) of the theory. Of course, significance of this observation (and 
its implications) is quite universal, by no means restricted to the field of 
matrix models themselves, however, it is yet not enough appreciated by the 
experts in other fields. In any case, we are going to deal only with matrix 
models in these notes. 

We proceed to consideration of the WI according to the following plan 
(not all the arrows will be actually discussed): see Fig.3. 



2.2 Virasoro constraints for the discrete 1-matrix model 



The basic example |26| , |27| which illustrates the arguments from the previous 
subsection is provided by the 1-matrix model 

Z N {t} = c N [ dHe^=o t ^ Hk . (2.1) 

JNxN 

This integral is invariant under any change of variables H — > f(H). It is 
convenient to choose the special basis in the space of such transformations: 

5H = e n H n+1 . (2.2) 

Here e n is some infinitesimal matrix and, of course, n > — 1. The value of 
integral can not change under the change of integration variable, and we 
obtain the identity: 



/ dH^**™" = [d(H + e n H n+1 )e£Zo^ H +^ Hn+1 

JNxN J 



r 

> 



i.e. 



dHe^o^ Hk V kt k TrH k+n + Tr%- = 0. (2.3) 
Vfc=0 bH ) 

In order to evaluate the Jacobian Tr * 5 ^"^ 1 let us restore the matrix indices: 



SH 

n 

n—k\ 



fc=0 fc=0 
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In Tr ^^ +1 one should take I = i and m = j, so that 

rrrn+l n 

Tr— — - = TrH k TrH n - k . (2.4) 

k=0 

Now we can note that since we started from Lagrangian of the most 
general form (consistent with the symmetry H — > UHU^), any correlation 
function can be obtained as variation of the coupling constants (all possible 
sources are included as counterterms) . In our particular example this is just 
a trivial remark: 

< TvH ai ...TiH an > = JdHe^o t ^ Hk TrH a \..TiH a - = 

= dt a ®ldt an Z N{t}- (2-5) 

We can use this relation together with (2T4 - ) in order to rewrite ( [2.3|) as: 

L n Z N {t} = 0, n > -1 (2.6) 

with 

oo o n rp 

L n = Y. kt k — h X] qT-qI ■ (2-7) 

k=0 ub k+n ^. = Q ub k ulj n—k 



Note that according to the definition ( p7l 

<9 



-Zn — NZn- 
dt 



Several remarks are now in order. 

First of all, expression in brackets in (|2.3j) represents just all the equations 



of motion for the model (|2.1|) , and ( |2.6j ) is nothing but another way to rep 



resent the same set of equations. This is an example of the above-mentioned 
identification of EqM and WI. 

Second, commutator of any two operators L n apearing in ( [2.6|) should also 
annihilate Z^{t}. It is another indication (not a convincing one, however) 
that we already got a complete set of constraints, that L n 's form a closed 
(Virasoro) algebra: 

[L n , L m ) = (n - m)L n+m , n,m>-l. (2.8) 
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Third, ( |2.6| ) can be considered as invariant formulation of what is Z^: it is 
a solution of this set of compatible differential equations. From this point 
of view eq.(|2.1D is rather a particular representation of Z N and it is sensible 
to look for other representations as well (we shall later discuss two of them: 
one in terms of CFT, another in terms of Kontsevich integrals). 

Fourth, one can try to analyze the uniqueness of the solutions to ( |2.6|) . 
If there are not too many of them the set of constraints can be considered 
complete. A natural approach to classification of solutions to the algebra of 



constraints is in terms of the orbits of the corresponding group |45fl . Let us 
consider an oversimplified example, which can still be usefull to understand 
implications of the complete set of WI as well as clarify the meaning of classes 
of universality and of integrability. 

Imagine, that instead of ( |2.6|) with L n 's defined in (|2.7|) we would obtain 
somewhat simpler equations: Q 

oo g 

l n Z = 0, n > with l n = V, ktk 



Then operator l\ can be interpreted as generating the shifts 

£2 — ► t 2 + eiii, 
*3 — > t 3 + 2eit 2 , 

We can use it to shift t 2 to zero, and eq. l\Z = then implies that 

Z{t 1 ,t 2 ,t 3 ,...) = Z{t 1 ,0,t 3 ,...) 

(t k = t k -&=&&=±, k>3). 

Next, operator l 2 generates the shifts 

^3 > ^3 + e 2^1; 

^4 — ^ ^4 + 2e2t2) 



5 One can call them "classical" approximation to ( |2.6| ), since they would arise if the 
variation of measure (i.e. a "quantum effect") was not taken into account in the derivation 
of (2.6). Though this concept is often used in physics it does not have much sense in the 
present context, when we are analyzing exact properties of functional (matrix) integrals. 
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and does not affect t 2 . We can now use eq. l 2 Z = to argue that 
Z(ti,t2, t 3 , i 4 , ...) = Z (ti, 0, t^, £4, ...) = Z (ti, 0, 0, £4, ...) 

etc. Assuming that Z is not very much dependent on tk with k — > 00, we 
can conclude, that 

Zfa, t 3 , t 3 , ...) = Z(ti, 0, 0, ...) = Z(l, 0, 0, ...) 

(at the last step we also used the equation l Z = to rescale t\ to unity). 

All this reasoning was correct provided t\ 7^ 0. Otherwise we would get 
Z(0, 1, 0, 0, ...), if tt = 0, t 2 ^ 0, or Z(0, 0, 1, 0, ...), if t x = t 2 = 0, t 3 ^ etc. 
In other words, we obtain classes of universality (such that the value of parti- 
tion function is just the same in the whole class), which in this oversimplified 
example are labeled just by the first non-vanishing time-variable. Analysis 
of the orbit structure for the actually important realizations of groups, like 
that connected to eq. (|2.7| ) has never been performed in the context of matrix 
model theory. It may deserve emphasizing that the constraints, as we saw, 
can actually allow one to eliminate (solve exactly) all the dependence on the 
time-variables, in less trivial examples they somehow imply the integrabil- 
ity structure, which is just a slightly more complicated version of the same 
solvability phenomenon. 



2.3 CFT formulation of matrix models 

Given a complete set of the constraints on partition function of infinitely 
many variables which form some closed algebra we can now ask an inverse 
question: how these equations can be solved or what is the integral repre- 
sentation of partition function. One approach to this problem is analysis of 
orbits, briefly mentioned at the end of the previous section. Now we turn 
to another technique 139], which makes use of the knowledge from conformal 



field theory. This constructions can have some meaning from the "physi- 
cal" point of view, which implies certain duality between the 2-dimensional 
world surfaces and the spectral surfaces, associated to configuration space of 
the string theory. However, our goal now is more formal: to use the means 
of CFT for solution of the constraint equations. 



6 This, by the way, is hardly correct in this particular example, when the group has no 
compact orbits. 
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This is very natural in the case when the algebra of constraints is Virasoro 
algebra, as in the case of the 1-matrix model, or some other algebra if it is 
known to arise naturally as chiral algebra in some simple conformal models. 
In fact the approach which will be now discussed is rather general and can 
be applied to construction of matrix models, associated with many different 
algebraic structures. 

We begin from the set of equations (|2.6|) which we shall further refer to as 
" discrete Virasoro constraints" . The CFT formulation of interest should pro- 
vide the solution to these equations in the form of some correlation function 
in some conformal field theory. Of course, it becomes natural, if we somehow 
identify the operators L n , which form Virasoro algebra with the harmonics 
of the stress-tensor T n , which satisfy the same algebra, and manage to relate 
the constraint that L n annihilate the correlator to the statement that T n 
annihilate the vacuum state. Thus the procedure is naturally split into two 
steps: First we should find a t-dependent operator (" Hamiltonian" ) H(t), 
such that 

L n (t)(e H ^... = (e H ^T n ... (2.9) 

This will relate differential operators L n to T n 's expressed through the fields 
of conformal model. Second we need to enumerate the states, that are anni- 
hilated by the operators T n with n > — 1, i.e. solve equation 

T n | G) = (2.10) 

for the ket-states, what is an internal problem of conformal field theory. If 
both ingredients H{t) and | G) are found, solution to the problem is given 
by 

(e m | G). (2.11) 

To be more explicit, for the case of the discrete Virasoro constraints we 
can just look for solutions in terms of the simplest possible conformal model: 
that of a one holomorphic scalar field 

00) = q + p\ogz + ~ir zk 

[Jn, Jm) = n$n+m,Q, [q,P] = 1- (2-12) 
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Then the procedure is as follows: Define vacuum states 



J*|0)=0, (JV|J_* = 0, k>0 
p|0) = 0, (N\p = N(N\, (2.13) 



the stress-tensor 



T{z) = \[d4>{z)Y = Y, T nZ~ n ~\ T n = Y,J-kJ k+n + \ E J aJb, (2.14) 

fc>0 a + b=n 

a,b>0 

and the Hamiltonian 

H(t) = = ^^WC*) 

U(z) = J2 k>0 t k z k , J(z) = d(j)(z). (2.15) 

It can be easily checked now that 

L n (N\e H ® ... = (N\e H ®T n ... (2.16) 

and 

T n |0)=0, n>-l. (2.17) 

As an immediate consequence, any correlator of the form 

Z N {t | G} = (N\e H ®G\0) (2.18) 

gives a solution to (|2.6|) provided 

[T n ,G}=0, n>-l. (2.19) 

In fact operators G that commute with the stress tensor are well known: 
these are just any functions of the "screening charges" f\ 

Q ± = jj ± = je^. (2.20) 

The correlator (|2.18| ) will be non-vanishing only if the matching condition 
for zero-modes of is satisfied. If we demand the operator to depend only 

7 For notational simplicity we omit the normal ordering signs, in fact the involved 
operators are : e H : and : e ±v/ ^* : 
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on Q + , this implies that only one term of the expansion in powers of Q. 



will contribute to ( 2.18Q , so that the result is essentially independent on the 



choice of the function G(Q + ), we can for example take G(Q + ) = e Q+ and 
obtain: 

Z N {t}~±(N\eP®(Q + f\0). (2.21) 

This correlator is easy to evaluate using the Wick theorem and the propagator 
<f>(z)(/)(z') ~ log(z — z') and finally we get 

Z N {t} = U N V e^^ U(z)mz) : lUikdZi : e^*> :| 0) = 

= m nili § Ci d^ U(zi) I&j(zi ~ z f (2.22) 

in the form of a multiple integral, which can in fact be directly related to the 
matrix integral in ( |2.1| ), see and the next section. 



Thus in the simplest case we resolved the inverse problem: reconstructed 
the integral representation from the set of discrete Virasoro constraints. How- 
ever, the answer we got seems a little more general than ([2.1|): the r.h.s. of 



eq. (|2.22j ) still depends on the contours of integration. Moreover, we can also 
recall that the operator G above could depend not only on Q + , but also on 
Q_. The most general formula is a little more complicated than Q2.22 ): 



z N {t | Q,c r } ~ {N + M)lMl ( N \^ m (Q + ) N+M (Q-) M \o) = 

-WTMmUi^Ri^- (2 - 23) 



nf +M nf(^-^ 2 



r J 



We refer to the papers for discussion of the issue of the contour-dependence. 
In certain sense all these different integrals can be considered as branches of 
the same analytical function Z^j{t}. Dependence on M is essentially elimi- 
nated by Cauchy integration around the poles in denominator in ( |2.23|) . 

Above construction can be straightforwardly applied to other algebras of 
constraints, provided: 

(i) The free-field representation of the algebra is known in the CFT- 
framework, such that the generators are polinomials in the fields <fr (only i n 
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such case it is straightforward to construct a Hamiltonian H, which relates 
CFT-realization of the algebra to that in terms of differential operators w.r.to 
the t-variables; in fact under this condition H is usually linear in t's and 
(/>'s). There are examples (like Frenkel-Kac representation of level k = 1 
simply-laced Kac-Moody algebras or generic reductions of the WZNW 
model |1|,|48|]-[|)1|]) when generators are exponents of free fields, then this 



construction should be slightly modified. 

(ii) It is easy to find vacuum, annihilated by the relevant generators (here, 
for example, is the problem with application of this approach to the case 
of "continuous" Virasoro and W^-constraints). The resolution to this prob- 
lem involves consideration of correlates on Riemann surfaces with non-trivial 
topologies, often - of infinite genus. 

(iii) The free-field representation of the "screening charges" , i.e. operators 
that commute with the generators of the group within the conformal model, 
is explicitly known. 

These conditions are fulfilled in many case in CFT, including conventional 
W-algebras |}2] and Af = 1 Q supersymmetric models 

For illustration purposes we present here several formulas from the last 
paper of ref.[39] for the case of the W r+ i-constraints, associated with the 
simply-laced algebras A of rank r. 

Partition function in such "conformal multimatrix model" is a function 
of "time- variables" tjj, , k = O...00, A = l...r = rankA, and also depends on 
the integer- valued r- vector N = {Ni...n r }. The W r+ i-constraints imposed 
on partition function are: 

W { n a \t)Zj^{t} = 0, n>l-a, a = 2...r + l. (2.24) 

The form of the PF-operators is somewhat complicated, for example, in the 
case of r + 1 = 3 (i.e. for A = A 2 (SL(3))) 

+ a ^ = }dt a dt b + di a diJ 



8 In the case of Af = 2 supersymmetry a problem arises because of the lack of reasonable 
screening charges. At the most naive level the relevant operator to be integrated over 
superspace (over dzd N 9) in order to produce screening charge has dimension 1 — ^A^, 
which vanishes when Af = 2. 
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wi 3) = e (km 



k,l>0 



+2E 

k>0 



E (kt 



dtk+n+l 
2 



- kt k lti 







a+b=n+k 



dt a dt b 



kt, 



dtk+n+l 
d 2 

di a di b 



- 2kt k lii — 







2kU 



dtk+n+l 

2 



)+ 



dt a dt b 



+ 



[2.26) 



+^ E 



5 3 



d 3 



a+fc+c=n 



dt a dt b dt c dt a dt b dt, 



■), 



and two types of time- variables, denoted through t k and t k - are ossociated 
with two orthogonal directions in the Cartan plane of A 2 : e = ^A, e = -^p- 



All other formulas, however, are very simple: Conformal model is usually 
that of the r free fields, S ~ / d<pd<pd 2 z, which is used to describe represen- 
tation of the level one Kac-Moody algebra, associated with A. Hamiltonian 



r+l 



A=l fc>0 



(2.27) 



where {/x A } are associated with "fundamental weight" vectors V\ in Cartan 
hyperplane and in the simplest case of A = A r (SX(r + 1)) satisfy 



r+l 



XX' 



r + 



p E = o, 



A=l 



thus only r of the time variables t^ . . .t^ r+1 ^ are linearly independent. Re- 
lation between differential operators (t) and operators in the CFT 
is now defined by 

Wi a \N\e H V... = (N\e H ®W ( : a) ... , 
a — 2, ... ,p; i > 1 — a, 



(2.28) 



where 



wl a) = & ^^Ui 



(2.29) 



9 Such orthogonal basis is especially convenient for discussion of integrability properties 
of the model, these t and i are linear combinations of time-variables t\ appearing in eqs. 
d^27|) and (|2^). 
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are spin-a generators of the algebra. The screening charges, that com- 

mute with all the W^ a \z) are given by 



Q(a) = I 



(2.30) 



{a} being roots of finite-dimensional simply laced Lie algebra A. 
Thus partition function arises in the form: 

Zfait} = (N\e H ^G{Q^}\0) 



(2.31) 



where G is an exponential function of screening charges. Evaluation of the 
free-feild correlator gives: 



z A N {t} - /n 



(A), 



i=l 



v A;/c>0 



X 



X 



nniK 

(a,P) i=lj=l 



(2.32) 



In fact this expression can be rewritten in terms of an r-matrix integral - a 
"conformal multimatrix model": 



Zir{t (a) } = ^n 1 I dH^...dH^ TT e Er=o4 a) ^) . 

■ [] Det ® I - I ® H {a+1) ) af3 (2. 



.33) 



K/3) 



In the simplest case of W3 algebra eq. ( |2.32|) with insertion of only two (of 
the six) screenings Q ai and Q a2 turns into 

Zfr uN2 (t,t) = j^^N^HQ^Y'iQ^fm = 

1 r 2 ' r _ (2-34) 

]J / dx^H / ^e c/ ^)A(x)A(x ! y)A( 2/ ), 



iVi!iV 2 ! 



where A(x,y) = A(x) A(y) Ylij(xi — yj). This model is associated with the 
algebra A = A2 (SL(3)), while the orig inal 1-matrix model (^2lD - (^23|) - 
with A = A l (SL(2)). 
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The whole series of models Q2.32H2.33| ) for A = A r (SL(r + 1)) is distin- 



guished by its relation to the level k = 1 simply-laced Kac-Moody algebras. 
In this particular situation the underlying conformal model has integer cen- 
tral charge c = r = rank.4 and can be "fermionized".^ The main feature 
of this formulation is that the Kac-Moody currents (which after integration 
turn into "screening charges" in the above construction) are quadratic in 
fermionic fields, while they are represented by exponents in the free-boson 
formulation. 

In fact fermionic (spinor) model naturally possesses GL(r+l) rather than 
SL(r+l) symmetry (other simply-laced algebras can be embedded into larger 
GL-algebras and this provides fermionic descriprion for them in the case of 
k = 1). The model contains r + 1 spin-1/2 fields ipi an d their conjugate ipi 
(b, c-systems); 



r+l 

s = E / 4>Md 2 z, 



r+l 

E 

central charge c = r + 1 , and operator algebra is 

ipj(z)tp k (z') = (z - z')5 jk : tpj(z)ip k (z') : + (1 - 5 jk ) : ij)j(z)il) k (z') : 
^j{z)i) k {z') ={z- z')5 jk : $j{z)$ k {z') : + (1 - S jk ) : ^j(z)ijj k (z') : 

The Kac-Moody currents of level k = 1 GL(r + 1) are just Jj k =: ipjip k : 
j, k = 1 . . . r+l, and screening charges are = iE^ § : ipjipk where E^ 
are representatives of the roots ct in the matrix representation of GL(r + 1). 
Cartan subalgebra is represented by Jjj, while positive and negative Borel 
subalgebras - by Jj k with j < k and j > k respectively. In eq.( p.23| ) Q + = 
z/fefe, Q- = i§^Mi while in eq.(£3§ Q (Ql) = i/fcfe, Q {a2) = 
i/^a (andQ^ =i§4> 2 i)3, = i § ifal)i, Q M = i $ Mi, Q K) = 

% § ^3^2) ■ Q( a ^ can be substituted instead of in ( |2.34j ) without changing 

10 This is possible only for very special Kac-Moody algebras, and such formulation is 
important in order to deal with conventional formulation of integrability, which usually 
involves commuting Hamiltonian flows (not just a closed algebra of flows) and fermionic 
realization of the universal module space (universal Grassmannian). In fact these restric- 
tions are quite arbitrary and can be removed (though this is not yet done in full details), 
sec section 4 below for more detailed discussion. 
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the answer. For generic r the similar choice of "adjacent" (not simple!) roots 
(such that their scalar products are +1 or 0) leads to selection of the following 
r screening operators — % § 4>i*p2 — — i§ ih^z-, — i§ i>zi)±, ■ ■ 
i.e. = i § for odd j and = — i § for even j. 

2.4 Gross-Newman equation 

We turn now to consideration of the WI for another sort of matrix mod- 
els. This subject concerns at least two important classes: the conventional 
discrete two-matrix models and Kontsevich models. As it was explained 
in the Introduction the theories of the second type arise in consderation of 
the (p, 1) continuous matrix models, as well as in the study of topological 
Landau- Ginzburg theories, while the two-matrix model is believed to exhibit 
a rich pattern of continuous limits and is capable to provide representatives 
of all the (p, q) universality classes (this line of reasoning, however, has never 
been really developed and we shall not discuss it in these notes). 

The starting point and the basic example is provided by the integral 

T v>n {L}= f dXe- trV(x ^ LX (2.35) 

Jnxn 

over nxn Hermitean matrix, which we shall further refer to as "Kontsevich 
integral" , keeping in mind its most important application (though this obvi- 
ous quantity has been of course considered by many other people). It may 
seem that the action in this integral is not of the most general type and we 
can no longer perform arbitrary change of variables X — > f(X) without 
changing the functional form of the integral. In fact this is incorrect, because 
"external field" L is matrix valued and coupled to X linearly, and therefore 
any correlator of X-fields can be represented through L-derivatives. Con- 
sider again the shift X — > X + e n X n+1 , n > — 1. Invariance of the integral 
implies: 

f dXe~ trnx)+trLX tr e n (-X n+1 V'(X) + LX n+1 + £ XhrX^A = 0, 
J V k=o J 



which can be rewritten as 




1 The obvious relation is used here: X 7 se tlLX — -]=- etrLX ■ Note that the order of 
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= tr e « (-afc) n+1 (V {at) - L ) = ( 2 -36) 

This system is in fact equivalent to a single matrix-valued equation 

d ' 



V 



dL t 



L T V {L] = 0. 



(2.37) 



As well as I know this equation was first written down in [53|, therefore it will 
be refered to as the Gross-Newman (GN) equation. It was rediscovered and 



implications for the theory of matrix models were investigated in p4T|,||30|| ,[|31j 

There are essentially two types of corollaries, which will be discussed in 
the next two subsections. First, GN equation can be used to characterize 
the function jFy{L} itself. This will lead us to consideration of Kontsevich 
models. Second, it can be used to derive equations for the 2-matrix model, 
which arises after Ty{L\ is further integrated with some weight over L. 



2.5 Ward identities for Generalized Kontsevich Model 

Being just the complete set of equation of motion the GN equation ( |2.37j ) 
provides complete information about the function Ty\L\. However, this 
statement needs to be formulated more carefully. A need for this comes, for 
example, from the observation that operators 

tlLm ( v '&- L ) < 2 - 38 ' 

do not form a closed algebra: their commutators have some different func- 
tional form. One of the reasons for these complications is that eq. ( [2.37 ) does 



not account explicitly for a very important property of Ty{L\. this func- 
tion in fact depends only on the eigenvalues of L. This information should 
be still added somehow to the GN equation. We shall analyze this issue of 
eigenvalue-dependence in more details in next sections. For our current pur- 
poses this argument implies that one should try to express equation ( 2.37] ) in 



matrix indices is reversed at the r.h.s. as compared to the l.h.s., i.e. derivatives are in 
fact w.r.to transported matrix L: f(X)e tlLX = /(^ _ )e trix (at least for any function 
f(x), which can be represented as a formal seria in integer powers of X). 
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terms of eigenvalues. Here, however, one should be carefull again. Clearly, 
not only Ty{L] depends on eigenvalues, it depends on their "symmetric" 
(Weyl-group invariant) combinations, i.e. it rather depends on quantities 
like trL a then on particular eigenvalues. Moreover, powers a here should be 
negative and fractional. 

Indeed, integrals like ( |2.35| ) are usually understood as analytical contin- 



uation from some values of parameters in the potential V, when integral is 
convergent. They can be also related to the formal (perturbation) seria aris- 
ing when integrand is expanded around a stationary point. To begin with it 
is reasonable to take n = 1 i.e. consider just an ordinary integral. For the 

sake of simplicity also take particular V(x) = —^f- Then the stationary 

i 

point is at x = \p and 

p+i _i p+i 

dxe- X — +lx ~ r^e^ P C >J p - ( 2 -39) 

fe>0 

It is now easy to understand what should be done in the general situa- 
tion with matrices and arbitrary potentials. First of all, one needs to solve 
equation for the stationary point, V'(X) = L. For this purpose it is most 
convenient to introduce a new matrix variable A instead of L, which by def- 
inition satisfies V'(A) = L. Then stationary point is just X = A. Second, 
one should separate the analogue of the complicated prefactor (quasiclassical 
contribution) : 

tr(AV"(A)-V(A)) 

Cy{A} = (2tt) b 12 (2.40) 
^detV'(A) 

Then the function that describes pure "quantum" contribution 
Z V {T} = Cv{A}" 1 ^V{y'(A)} (2.41) 

12 The "classical action" in fl2.40p can be also represented as tr (AV'(A) — V(A)) = 
tr J AdV'(A). Determinant of quadratic fluctuations is defined as 

(2tt)" 2 / 2 (detT/"(A)r 1/2 ~ / dYe- trV ^ Y \ 



where V^A, Y) = lim e _ i (V(A + eY) - V(A) - eV'(A)Y). For V(A) = we have 

V"(A) = (5^fc=o ^ k ® A p_fc_1 ). One could easily choose an "opposite" parametrization 



in eq.( 2.42| ): = —jtrA k . Though not quite obvious, this never influences any results 



(see section 2.10 for an example). Our choice of signs is motivated by simplification of 
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to be refered as partition function of the Generalized Kontsevich Model 
(GKM) [0O|1 ) can be represented as a formal (perturbation) series expansion 
in variables 

T k = ^trA~ fc . (2.42) 

GN ( [2.37D equation can be now rewritten as a set of differential equations 
for Z V {T}. Indeed, we already have: 

C v (V (^-) ~ l) C V Z V {T} = 0, (2.43) 



dL 



tr 



but it is still necessary to express the operator at the l.h.s. in terms of T. 
This is in fact possible to do, using the relation: 

^-m^H (2.44) 



dL tr k 0L tr <9Tfc 

and substituting the traces of A-matrices, which can arise in the process of 
calculation, by T's. It is important only that A's usually appear in negative 
powers: this is already achieved by the choice of a proper normalization 
factor Cy{A}. For monomial potential V p (X) = this is especialy simple: 
L = A p and %f*- = -±A"P- fc . 

OLtr P 

This reasoning allows one to rewrite eq. (|2.43 ) identically in the form 



J2A' l O l (T)Z v {T} = 0, (2.45) 
i 

where Oi are some differential operators, depending on the shape of V, but 
independent on the size n of the matrix (as all the above reasoning never 
refered to particular values of n, except for a sample example at the very 
beginning). It remains to use the fact that matrix L can be arbitrary large 
and have arbitrarily many independent entries, in order to conclude that we 
derived a set of constraints on Zy in the form 

O^ZyiT} = 0. (2.46) 



formulas for the GKM, including the relation between L and A. Instead, some sign factors 
appear in formulas, related to Toda-like representations of partition functions and those 
involving VK-operators. 
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For potential V of degree p + 1 these appear to be exactly the "continu- 
ous Virasoro constraints". See refs.[24] and [ 30 ] for detailed analysis of the 
Virasoro case p = 2 (associated with the pure topological gravity and with 
the double-scaling limit of the 1-matrix model), and |55| for the exhaustive 
presentation of the case of p = 3. 



2.6 Discrete Virasoro constraints for the Gaussian Kont- 
sevich model 

As a simplest illustration of the technique, described in the previous subsec- 
tion, we derive now the constraints for the Gaussian Kontsevich model 
with potential V(X) = \X 2 : 



(detL) 



N 



dX(detXy y e 



-+LX 



(2.47) 



In this case L = V'(A) = A, and the time- variables are just 



T, 



-trA _fc = -trL k . 
k k 



(2.48) 



To make the model non-trivial an extra "zero-time" variable |3g] is in- 
troduced, which was not included into the previous definition ( |2.41| ). Now 
note that the A^-dependence of Kontsevich integral (|2.35|) can be described 
simply as an extra term in the potential: V(X) — > V(X) = V(X) — NlogX 
(though this can not be done neither in the quasiclassical factor Cy nor in 
the definition of time- variables T). Since the GN equation depends only on 
Kontsevich equation, we can use it with V substituted by V . Then we have 



instead of ( ggg) : 



T * 

e- trI r 







n+l 



{detL) N \dL tr/ 







dL 



(detL) N e +tI -Zx2{N,T} = 0. 



N 







-i 



tr 



dL 



tr , 



L 



(2.49) 



In order to get rid of the integral operator (^-) 1 one should take here n > 
rather than n > — 1. In fact all the equations with n > follow from the one 
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with n = 0, and we restrict our consideration to the last one. For n = we 
obtain from (|2.49 ): 








2N 




or 

d 



dL tr 

and it remains to substitute: 

1 az . 




0, (2.50) 



dZ 



d 2 z 



dL tr ~ ^k=0L*+i d T k i 

Eoo 1 ( / f r 1 A 9-Z I y-*m— 1 

m=-l L'"+ 2 \^fc>max(m,0) \ lL L k ~ m ) 8T k ^—'k=X dT k dT m _ k 

and finally obtain: 



d 2 z 



oo 
m=— 1 



27V^- + jV 2 5 m , - iV (tri) 5 mi _x) Z 




= E™=_i ^e^L m (T + r)e~ NT °Z = 0. (2.51) 
Here L m {t) are just the generators ( |2.7| ) of the discrete Virasoro algebra 

e^L m (t)e"^ = e Nt ° ( £ H fc ^- + £ ) e""* (2.52) 

and at the r.h.s. of (|2.51 ) 

13 This small correction is a manifestation of a very general phenomenon: from the 
point of view of symmetries (Ward identities) it is more natural to consider Zy not as a 
function of T- variables, but of some more complicated combination Tk+rk, depending on 
the shape of potential V. If V is a polinomial of degree p+ 1, f k = |tr(V {Xj)~ k/p , while 
l— — 

Tk = fc(p-fc) ^ cs (^'(z 1 )) p ^Z 1 - F° r monomial potentials these expressions become very 
simple: Tk = Tk and — — ^^-Sk, p +i- See Q and section 4.9 below for more details. 
In most places in these notes we prefer to use invariant potential-independent times Tk, 
instead of Tk, but then Ward identites acquire some extra terms with ru (which in fact 
will be very simple in our examples, which are all given for monomial potentials). 
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Thus we found that the WI of the Gaussian Kontsevich model ( p. 47] ) 
coincide with those of the ordinary 1-matrix model, moreover the size of the 
matrix N in the latter model is associated with the "zero-time" in the former 



one. This result [56] of course implies, that the two models are identical: 



e- Nl0 Z^ {N, T 1; T 2 , . . .} ~ Z N {T , T x , T 2 , . . .}. (2.53) 

2 

We shall discuss direct connection between these two matrix integrals (|2.1| ) 
and ( 2.47|) in the next section, after some more details will be presented 



about the structure of "eigenvalue" matrix models. 

2.7 Continuous Virasoro constraints for the V = ^- 
Kontsevich model 

This example is a little more complicated than that in the previous subsec- 



tion, and we do not present calculations in full details (see |24| and pOfl). Our 
goal is to demonstrate that the constraints which arise in this model, though 
still form (Borel subalgebra of) some Virasoro algebra, are different from 
(|2.6|) . From the point of view of the CFT-formulation the relevant model is 
that of the twisted (in this particular case - antiperiodic) free fields. These so 
called "continuous Virasoro constraints" give the simplest illustration of the 
difference between discrete and continuous matrix models: this is essentially 
the difference between "homogeneous" (Kac-Frenkel) and "principal" (soli- 
ton vertex operator) representation of the level k = 1 Kac-Moody algebra. 
From the point of view of integrable hierarchies this is the difference between 
Toda-chain-like and KP-like hierarchies. We shall come back to a more de- 
tailed discussion of this difference later, when the "multi-scaling continuum 
limit" will be considered. 

Another (historical) aspect of the same relation also deserves mentioning, 
since it also illustrates the interrelation between different models. The dis- 
crete 1-matrix model arises naturally in description of quantum 2d gravity as 
sum over 2-geometries in the formalism of random equilateral triangulations. 
The model, however, decribes only lattice approximation to 2d gravity and 
(double-scaling) coninuum limit should be taken in order to obtain the real 
(continuous) theory of 2c? gravity. This limit was originally formulated in 
terms of the contraint algebra (equations of motion or "loop" or "Schwinger- 
Dyson" equations - terminology is taste-dependent), leaving open the prob- 
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lem of what is the form of partition function Z cont {T} of continuous theory. 
Since the relevant algebra appeared to be just the WI for Kontsevich model 
(with V(X) = this proves that the latter one is exactly the continuous 
theory of pure 2d gravity. At the same time, Kontsevich model itself can be 
naturally introduced as a theory of topological gravity (in fact this is how 
the model was originally discovered in f2~2||). From this point of view the 
constraint algebra, to be discussed below in this subsection, plays central 
role in the proof of equivalence between pure 2d quantum gravity and pure 
topological gravity (in both cases "pure" means that "matter" fields are not 
introduced). 

After these introductory remarks we proceed to calculations. Actually 
they just repeat those for the Gaussian model, performed in the previous 
subsection, though formulas get somewhat more complicated. This time we 
do not include zero-time N and just use eq.( [2.37| ) with V(X) = 4p. Now 
it is also much more tricky (though possible) to work in matrix notations 
(because fractional powers of L will be involved) and we rewrite everything 
in terms of the eigenvalues of L. 



We substitute 



C 



3 



i-r 2 X 3 / 2 



n 7 ,5(v / ^5 + v^) 



7 

and introduce a special notation for 



d 2 \ d 2 v 1 (_d d_ 



v -r-i 9 r 9 rr 1 1 v 



Then ( p.37| ) turns into 



D \ 2 ^ 1 ( V 




Er-rb-^ 2^{T} = 0. (2.54) 



Now we need explicit expression for T: 

1 
k 



T k = (2.55) 
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and as we already know from the previous subsection we also need 

2 



-5, 



fc,3- 



(2.56) 



It will not be explained untill we turn to consideration of integrable struc- 
ture of Kontsevich model in the following sections, but Zx±{T} is in fact 

3 

independent of all time- variables with even numbers (subscripts). Therefore 
we can take only k = 2a + 1 in (|2.55 ), 



T 2a +i — £<5 A 5 
T"2a+l — ~~ 1^0,1 



(2.57) 



and 







<9A 7 ~ 
d 2 



dT 2a+1 dZ 

d\-y dT 2a + 1 



1 y^°° \ a 2 dZ 
■ 2 2^a=0 A 7 QT 2a+1 > 



{T} 



1 v^°° 



A. 



-a— 6— 3 



dZ 



- v c 

4 Z^ a ,fe=0 /v 7 0T 2 a+idT 2M 



These expressions should be now substituted into (|2.54|) and we obtain: 

1 00 87. 
4 2^ A 7 

41 a,fe=0 



+ E 

a=0 



dT 2a+ i dT 2 b+ 1 
1 



+ 



Z r, Z Z A^v — Ag 



a=0 



1 



A 7 2 — A^ 2 



<9Z 



7 ^ <5^ 7 -y/A 7 (vA5 + -y/A 7 ) 



(2.5* 



<9T ; 



2a+l 



\n+2 
n=-l A 7 



31 



with 



C 



2n 



a=Q 



E « + 



x 9 | 1 E 

dT 2a +2n+l 4 a+6=n _ x 

a,fe>0 



(T 2a+ i + r 2a +i) x 

1 „ 1 . 



- E k ( T k + r k ) 



d 



2n-l 



dn 



d 2 1 1 



odd fe=l ^ J fc+2n ^ odd fc 

Factor | in front of the first term at the r.h.s. in (|2.59|) is important for £ 2 „ 
to satisfy the properly normalized Virasoro algebra:^ 



[C2ni & 



2 m 



(n - m)£ 2n +2m- 



Coefficient \ in front of the second term can be eliminated by rescaling of 
time- variables: T — * |T, then the last term turns into j^Tfd n -i. 

We shall not actually discuss evaluation of the coefficient in front of Z 
(with no derivatives), which is denoted by [ 



in 



sec 



2| and 



In fact almost all the terms in original complicated expression cancel, giving 
finally 



1 



16A2 



+ 



It 
4Ar 



and this is represented by the terms with S nt o and 5 n _i in expressions ( |2.59| ) 
for the Virasoro generators £ 2n . 

The term with the double T-derivative in fl2.58| ) is already of the necessary 
form. Of intermidiate complexity is evaluation of the coefficient in front of 
- in ( |2.58| ), which we shall briefly describe now. First of all, rewrite this 
coefficient, reordering the items: 



(a + -)A 7 ' 



+ 



E 



a: 



a; 



-a; 



4 



E 

5^7 



\-a-2 



-a-l 



(2.60) 



14 



Therefore it could be reasonable to use a different notation: C n instead of Cm- We 
prefer because it emphasises the property of the model to be 2-reduction of KP 

hierarchy (to KdV), see section 4 below. 
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The first two terms together are equal to the sum over all j (including j = i): 



a+4 ,a+| 



2 s - \ ' 2 i A-,, — A 

\ a+2 - , , 

v 7 ^7 ^8 



1 ^A^-AlAf 1 1 



E 



9\a+2 \ _ \ o+I ' 

ZA 7 5 A 7 A S \ s 2 



Similarly, the next two terms can be rewritten as 



1 A~ a ~ 2 1 \ K 



1 3 

1 \ 2 \ a+ 2 \ a + 2 1 

1 \ - A 7 Ag — A 8 1 

_ 2A a + 2 ^ A - Aa 

ZA 7 <5 A 7 A <5 A 5 2 

The sum of these two expressions is equal to 

1_ A»+ 2 -A£+ 2 1 

ia+2 2^ 



ZA 7 5 A 7 A <5 A 5 2 

Note that powers a + 2 are already integer and the remaining ratio can be 
represented as a sum of a + 2 terms. Adding also the last term from the l.h.s. 
of fl2.60|) , we finally obtain: 

1 1 a 1 1 
— + ~ V — — V - 

A?; +1 2 „t^ A™+ 2 V \°-' - 



7 -»=-i^r ^ a: ^ 2 



9 E -^(2a-2n + l)(T + r) 2a-2n+l 
Z n=-l A 7 



in accordance with (|2.58|) and fl2.59|) . 



2.8 VK-constraints for the asymmetric 2-matrix model 

We turn now to a very different application |38[] of the GN equation (|2.37 ). 
Namely, we shall now consider jFy n {L} as a building block in construction 
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of conventional discrete two-matrix model 

Z N {t,t} = c 2 N J dHdHe^ {t ^ Hk+l ^ Rk ^ HR 



J dLe^^Tu^L}. (2.61) 



Now L plays the role of H and U{H) = £ fe t k H k . 

We can now use GN equation to derive a relation for Z]y{t,i}. Take 



[u(-^-) + LjF i j )N {L} = 0, (2.62) 

multiply it by e 1 * 11 ^ = e^k tkTvLk and integrate over L. In order to express 
this relation in terms of t-derivatives of z it is necessary to have some "scalar" 
rather than matrix equations, therefore we" 11 actually need to take trace of 
). However, in order not to loose any information, we first multiply 



( |2.62| ) by L n and then take the trace. In this way we obtain: 

J dLe^ tkTrLk TrL n (u(-^-) + Lj T d {L\ = 0. 
Integration by parts gives: 

J dLF {L}Tr (fj(-JL) + A L n e^ tkTrLk . (2.63) 



Now we need to introduce a new class of operators p8[ . Consider the 
ion of Tr^ 
of terms like 



action of Tr^^-L n on e TlU ^ = e^^ tkTrLk . It gives some linear combination 



tiL a \..tiL a 'e tvU{L) = - d \ e~ trU{L) 

ot ai ...dt ai 

i.e. we obtain a combination of differential operators with t-derivatives, to 
be denoted W(t): 

W^\t)e^ = Tr-^L n e^ L \ m, n > 0. (2.64) 
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For example, 



oo Q n pp. 

WV = Ekt k -^- + Y,aA^, n>-l; 

oo o oo 

*?'=Etw*.flr- + E«* E a^r+ (2 ' 65) 

fc,i=l ^"<fc+Z+n fe=1 a +fe=fc+n UL a^ L b 

j-V^ >T ^ V g? ■ (" + !)(" + 2) g 



fc=l a+b=n+l dtadtb+k-1 a +b+c=n dt a 8t b dt c 2 <9t n 



Note, that while and W^ 2 ) are just the ordinary ([/(l)-Kac Moody and 
Virasoro operators respectively, the higher WA^-operators do not coincide 
with the generators of the W-algebras: already 

OO O OO !^)2 

w® + w m = J2 mu-^— + 2 E kt h £ ( 



/,./_!. dtk+l+n k= i a+b=k+n ®t a dt b 

E 

a+6+c=n 



3 _ , rri dt a dt b dt c 



W-operators (in variance with ordinary VT-operators) satisfy recurrent rela- 
tion: 



wt +1) = Ekt k Kt+ E w'™™> n ^~ m - (2 - 66) 

fc=l k=0 ° l k 



Actually not too much is already known about the W operators and the 
structure of W-algebras (in particular it remains unclear whether the neg- 
ative harmonics W^ 71 ^ with n < — m can be introduced in any reasonable 
way), see |38| for some preliminary results. 

Equation (|2.63|) can now be represented in terms of the VF-operators: 



d_ 



Y.i-r'khWZ-k + ] Z N {t,t} = 0. (2.67) 

,fc>i 
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This relation is highly asymmetric in t and t, and in fact it provides a suitable 
description of the WI only in the somewhat peculiar case when potential 
0(H) is a polinomial of finite degree. See refs. [pwj] and |J8] for discussion of 
such asymmetric models. 

2.9 VK-constraints for generic 2-matrix model 



When both potentials U and U in fl2.61| ) are generic formal seria, eqs.( |2.67| ) 



represent only a one-parametric subset of the 2-parametric family of WI. 
Before we describe the whole set, let us emphasize that the two-matrix model 
( |2.61|) is the one, where the action is not of the most general form, consistent 



with some symmetry. Therefore it is not covariant under arbitrary change 
of variables H, H — > f(H, H), f(H, H), and our usual method of derivation 
of Ward identities does not work. The reason why generic 2-matrix model 
with action containing all the possible combinations Tr (H 0-1 H bl H 0-2 H b2 ...) 
is never considered seriously is essentially our poor understanding of the 
unitary- matrix interals for "non-eigenvalue" theories, to which class such 
generic model belongs. For reasons to be explained in the next section such 
problems do not arise for the models of the form ( [2.6 1|) or ( [2.331 ), an d this is 
why they attracted most attention so far. Hopefully the problems with the 
unitary-matrix integrals are temporal and this restricted class of multimatrix 
models will be unlarged, this should be especially easy to do in the part of 
the theory dealing with constraint algebras, but this subject is beyond the 
scope of the present notes. 

In order to derive the complete set of WI for the model ([2.611) , we apply 



the following semi-artificial trick. Note that exponential e TrHH satisfies: 

(0777 CYn \ 

TrH n ^- T T H m ^L- e TlHR = 0. (2.68) 

Let us integrate this identity over H and H with the weight e TrU ( H )+ TrU ( H ) 
and then integrate by parts. We obtain an identity: 




(2.69) 

IV (__|_) H»-Tr(-^L-\ H m )e T ^ +T ^ 
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which can be represented in terms of W operators [B(J : 

wi™+ 1 \t)Z{tJ} = (-) m - n W^n\i)Z{t,t}, forallm,n>0. (2.70) 

This is the full(?) set of WI for the 2-matrix model. When one of poten- 
tials (say, U(t) is polinomial of finite degree, the most of this symmetry is 



"spontaneously broken", the surviving part being described by eqs. ([2.67|) . 

Among other things eq. fl2.7Q|) reveals an amusing automorphsim of the 
Wqo algebra: 

W^^W%*\ m,n>0, (2.71) 

for example, Virasoro's Borel subalgebra is formed not only by operators 
but also by W^ 2 \ n > — 1 (while U(l) Borel subalgebra - not only 
by = ^, but also by W^ n+l \ n > 0). 

One can attempt to apply the same procedure and derive V^-identites for 
the conventional (p — l)-matrix models with p — 1 > 2. In principle, this 
is possible, but unfortunately the arising equations neither have a nice form 
nor is there many enough of them. However, for illustrational purposes we 
shall scetch some relevant fromulas in the rest of this subsection. 

Consider the multimatrix integral 



Z = J dH l ...dH p _y 



. e TrC/i(Hi)+...+Trf/ J1 _ 1 (if J ,_ x ) e Tr(H 1 H2+H 2 H a +...+H p ^ 2 H p ^ 1 ) 



Acting on Z, operator W^H^\t^) produces an insertion of TrH™ 
at the position of ... in ( 2.72 ). Integration by parts gives: 



(2.72) 







TrH? (--^L-\ — ► (-) m Tr#?#™ = {-) m TrH™H? 



tr 



15 Relations ( 2.68| ) and thus ( 2.70 ) are in the obvious sense associated with TrH n H m . 
Of course there are similar relations, in the same sense associated with any object like 
Tr(H ai H bl H a2 H b2 ...) and with products of such traces: it is enough to substitute all 
H — > g ® t to obtain the l.h.s. of the equation and substitute all H — > ^ to obtain its 
r.h.s. (one should only remember that such substitution is possible, say in the l.h.s. if all 
the H are put to the right of all H, in order to restore the matrix-product form of the 
relation, one should carefully take into account all the commutators, arising when s ® is 
carried back to original position of the corresponding H). All such relations can appear 



to be just implications of eq.(2.7C) 
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In the case of p — 1 = 2, that we discussed above, this can be rewriten as 

/ d* \ n ( 1) V 

-) m TrH™ — — and integration by parts gives (-) mr M 2 m ' > 



\dH 2 ,tr J V ®H2,tr J 

what is equivalent to the action of (— ) m+n wjn-n\t^) on Z: we reproduce 
equation fl2.70| ). 

However, for p— 1 > 2 things are more complicated. Insertion of Tr H™H™ 

is equivalent to that of TrH™ — H 3 ] , which after integration by 

\0H2jr 

parts and acts on e U2( - H2 ^ and gives: 



M " 1 - oL - H > 



V k 



dH> 



( 2 -73) 



2,tr 



Derivatives remaining at the r.h.s. should be carried through the first bracket 
and than act on e u ^ H ^ etc. After all we get some linear combination of terms 
like TrH^H^H^H^ 2 . . . with t^-dependent coefficients. 

Now, if we are dealing with the p — 1 = 3- matrix model, every H 2 standing 
to the right of all H^s can be substituted by 9 ^ , otherwise one should also 

include terms with commutators when this tttt — is carried back to the place 
where H 2 was standing. In this way we obtain a combination of insertions of 
the form 



(2.74) 



The resulting operator can be expressed through W(t^) and we obtain an 
identity, saying that some algebraic combination of W(t^') and W(t^) with 
t^) -dependent coefficients annihilates the partition function. 

For p — 1 > 3 insertion of H 2 is equivalent to that of -rrS rather 

than - , and the procedure should be repeated again and again. Finaly 
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one arrives at constraints where the operators are algebraic combinations of 
W(t w ) and W(t ip -^) with the coefficients, which depend on t&, . . . ,&~ 2) 
(moreover these are infinite seria in W operators, unless all the intermediate 
potentials U 2 , ■ ■ ■ , U p -2 are polinomials of finite degree. 

This is of course not a too illuminating procedure and in fact it was 
never worked through to get concrete identities in any nice form. Instead 
it can serve to illustrate the problems, peculiar for the class of conventional 
multimatrix models (at least for p — 1 > 2). It can also emphasize the beauty 
of conformed multimatrix models, which have clear advantages already at 
the level of Ward identites. 



2.10 VK-operators in Kontsevich model 

After VT-operators are introduced, we can also rewrite the GN equation ( j2.43j ) 
for Kontsevich models in terms of \Vs. Namely, we shall prove the following 
identity fl38| : 




±)^^A-i-iw££ 1 \T)Z{T h }, (2.75) 



valid for any function Z which depends on = =F^trA k , k > 1 and 
T = itrlogA with n x n matrix A. Application of the identity ( [2.75 ) is 



most straightforward in the Gaussian model ( |2.47| ), e.g. for transformation 
of eq. (|2.50|) into eq. Q2.51| ) (we remind that L = A in this case). In other 



cases calculations with the use of identity ( 2.75|) , accounting for the quasi- 



classical factor Cy{L} and the difference between L = V'(A) and A become 
somewhat more involved, though still seem enough straightforward. Also for 
particular potentials V(X) partition function Z V {T} is actually independent 
of certain (combinations of) time- variables (for example, if V(X) = it 
is independent of all T p k, k G Z + ), and this is important for appearence of 
the constraints in the standard from like eqs. [2.58|) , ( 2.59f) , i.e. for certain 



reduction of iy-constraints to the ordinary VF-constraints. This relation 
between W- and iy-operators deserves further investigation. 

The proof of eq.( [2.75| ) is provided by the following trick. Let us make a 
sort of Fourier transformation: 

Z{T} = JdH g{H}e^=o T ^ H \ (2.76) 



39 



where integral is over NxN Hermitean matrix if .R Then it is clear that once 
the identity (|2Tg ) is established for Z{T} substituted by e TvU(H \ U(H) = 
Y^T=o TkTrH k , with any matrix H, it is valid for any function Z{T}. The 
advantage of such substitution is that we can now make use of the definition 
( |2.64|) of the W operators in order to rewrite ( [2.75|) in a very explicit form: 

/ Q \ m + 1 OO 

° p TrU(H) _ mrofl a -Z-lT^(m+l)/ T x TH7(H) _ 



OA 



tr / l>0 




e TrU(H)_ (277) 



Now expression for T's in terms of A should be used. Then 

e TrU(H) = Det ±l( A &I-I&H) 

and substituting this into (|2.77|) we see that ( ^.75[ ) is equivalent to 

rj \ m+1 / f) \ m 1 

' (±) m+1 I-Tr' ' ' 



OAtrJ y ' \dH tr J A&I-I&H 

■Det^iA ® I - I ® H) =0 

Here " Tr" stands for the trace in the if-space only, while Det = Det ® det 
- for determinant in both H and A spaces. After one A-derivative is taken 
explicitly, we get: 



(J® Tr) 




(2.78) 



This is already a matrix identity, valid for any A and H of the sizes n x n 
and NxN respectively. For example, if m = (W^-c&se), it is obviously 

16 Here it is for the first time that we encounter an important idea: matrix models - 
the ordinary 1-matrix model (2.1) in this case - can be considered as defining integral 
transformations. This view on matrix models can to large extent define their role in the 
future developement of string theory. 
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satisfied. If both n = N = 1, it is also trivially true, though for different 
reasons for different choice of signs: for the upper signs, the ratio at the l.h.s. 
is just unity and all derivatives vansih; for the lower signs we have: 



dh 



E 



Va + b— m— 1 
a,b>0 



dX 



dh 



dX + dh 



and this obviously vanishes since + Jr)/(A — h) = for any /(x). 

If m > and A, if are indeed matrices, direct evaluation becomes much 
more sophisticated. We present the first two nontrivial examples: m = 1 and 



m 



2. The following relartions will be usefull. Let Q = 

Det ±l Q^-Det^Q = ± [(I ® Tr)Q] ; 
Det ±l Q^-Det^Q = t [(tr ® J)Q] ; 



1 



Then 



9 



9A tr 
I J® 



d \ 



OH, 



= -[(tr® I)Q)Q; 
Q = [(/ ® Tr)Q] Q. 



(2.79) 



if , 



This is already enough for the proof in the case of m — 1. Indeed: 



Det ±x Q 



d 







QDet Tl Q 



= {- [(tr ® I)Q] Q ± [(I ® Tr)Q] Q} =f 
T{[(I ® Tr)Q] Q =f [(tr ® /)Q] Q} = 0. 

The first two terms at the r.h.s. come arise from A-, while the last two - from 
if -derivatives. 

In the case of m = 2 one should take derivatives once again. This is 
a little more tricky, and the same compact notation are not sufficient. In 
addition to ( |2.79| ) we now need: 



'_d_ 
dA 



® / [(tr ® I)Q] Q = - [(tr ® J)Q] 2 Q - 8. 



(2.80) 



fr 
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Here 



[(tr ® I)Q] 2 = [(tr <g> /) [(tr <g> J)Q] Q] , (2.81) 



while in order to write *B explicitly we need to restore matrix indices (Greek 
for the A-sector and Latin - for the H one). The (m, 7&;)-component of ( |2.80| ) 
looks like: 



d 



9A.0O, 



6 im QliQ% = -Q%Q^Q% - Q%Q l LQ% (2.82) 



and appearence of the second term at the r.h.s. implies, that B 1 * = QspQasQpy 
Further, 

(ait ® J ) t( J ® Tr )^ Q = 

- [(/ ® Tr) [(tr ® I)Q] Q] Q - [(/ ® Tr) [(/ ® Tr)Q] Q] Q; 

( J ®db) [(tr®/)Q]Q = 

+ [(tr <g> /) [(/ <g> Tr)Q] Q] Q + [(/ ® Tr) [(tr <g> J)Q] Q] Q; 

(/ ® dfc) [(/ ® Tr)Q] Q = + [(J ® Tr) [(/ ® Tr)Q] Q]Q + B. (2.83) 

It is important that B that appears in the last relation in the form of B^L = 

QasQspQpi * s exac tly the same B as in eq.( |2.80| ). 
Now we can prove ( 2.78| ) for m = 2: 

Det±1 Q {{atY ® J " J ® (aft) 2 ) QDet Tl Q = 
= {± [(/ ® Tr)Q] (- [(tr ® I)Q] Q ± [(J ® Tr)Q) Q) - 

- (- [(tr ® I) [(tr ® I)Q] Q]Q-B)± 

± (- [(/ <g> Tr) [(tr (g) J)Q] Q] Q - [(tr <g> /) [(/ <g> Tr)Q] Q] Q)} - 

- {t [(tr ® I)Q] ([(/ ® Tr)Q] Q =p [(tr ® J)Q] Q) + 
+ ([(J ® Tr) [(/ ® Tr)Q] Q] Q + B) =f 

=F ([(tr ® /) [(/ (g) Tr)Q] Q] Q + [(/ ® Tr) [(tr ® I)Q] Q] Q)} (2.84) 

where the terms 1,2,3,4,5,6 in the first braces cancel the terms 1,3,2,4,6,5 in 
the second braces and identity (|2.81|) and its counterpart with (tr ® I) — > 
(/ ® Tr) is used. 

Explicit proof of eq. (|2.78|) for generic m is unknown. 
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3 Eigenvalue models 



3.1 What are eigenvalue models 



Given the present state of knowledge we need to consider in most cases only 
the narrow class of the "eigenvalue" models. These models have the property 
of being associated with conventional integrable hierarchies (of (multicom- 
ponent) KP and Toda type), where integrable flows just commute (instead 
of forming less trivial closed algebras), and thus with the level-1 Kac-Moody 
algebras (by artificial tricks, familiar from the bosonization formalism in 
conformal field theory |58| these can be sometimes generalized to particular 
other levels like k = 2). This means that the models are essentially associated 
with abelian Cartan subalgebras rather than with full matrix algebras. [[] In 
CFT-formulation (see below) this means that the eigenvalue models can be 
represented in terms of the free fields, which bosonize the Cartan subalgebra 
of the whole group in the WZNW model (the remaining (/3, 7)-fields |16[] be- 
ing (almost) neglected - their remnants are observed in the form of "cocycle" 



factors in the Frenkel-Kac formulas 47], see |pq| ). In the matrix-integral 



representations the integrals for the eigenvalue models are in fact reduced to 
those over diagonal matrices (consisting of eigenvalues of original matrices, 
thus the name "eigenvalue models"). 

Most important, from the physical point of view eigenvalue models de- 
scribe only topological (discrete) degrees of freedom, but not any propagat- 
ing particles. Q This can be understood if one notes that matrix models 
usually possess gauge symmetry, associated with the unitary rotation of ma- 
trices, M a — > U^ a M a U a) i.e. matrix models are usually gauge theories. 
In the case of eigenvalue models this symmetry is realized without "gauge 
fields" V a p, which would depend on pairs of indices a, (3 and transform like 

17 Groups, arising in the theory of matrix models and integrable hierarchies are not just 
those of matrices, appearing in the integral representations: the latter ones are in the best 
case related to the zero-modes of the former. Moreover, even this relation is not usually 
simple to reveal. This remark is important to avoid confusion in the next paragraphs. 

18 Particles are always related to the "angular (unitary-) matrix" integrals (as everybody 
knows from example of the Wilson lattice QCD) which are highly less trivial to deal with, 
though these are also integrable in some broader sense of the word - within the (yet non- 
existing) generalization of integrable hierarchies from the fields in the Cartan subalgebra 
to the entire WZNW model. 
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V a /3 — > U^VapUp. In other words, eigenvalue models are gauge theories 
without gauge fields, i.e. are pure topological. Thus it is not a surprise that 
they usually live in the space time of dimension d < 2, Q since for d > 2 there 
should be particles, associated with the gauge fields. At the "boundary" lies 
the model of "d = 2 (c = 1) string", which has one particle-like degree of 
freedom (dilaton, which becomes tachyon in the d > 2 models). This very 
interesting model is much worse understood than the d < 2 models, at least 
its properties are already somewhat different from other eigenvalue models 
(especially in the most interesting "compactified" case), and it will not be 
discussed in these notes. Later we" 11 return to the subject of non-eigenvalue 
(d > 2) theories, though not too much is yet known about them, but now we 
are going to concentrate on the eigenvalue models. 



3.2 1-matrix model 

Hermitean matrix integrals are usually transformed to the eigenvalue form 
by separation of angular and eigenvalue variables. As usually, the simplest 
is the case of the 1-matrix model 

Z N {t} = c N f dHe^=o tkTvHk , (3.1) 

JNxN 

where this separations does not involve any information about unitary-matrix 
integrals. Take 

H = U ] DU, (3.2) 

where U is a unitary matrix and diagonal matrix D = diag(/ti.../tjv) has 
eigenvalues of H as its entries. Then integration measure 

N \dU] N 

dH= J] dR l3 = J 1] dhA 2 (h), (3.3) 

i,j=l [dUCartanl £_i 

19 Let us remind that in the Polyakov formulation which is the least counterintuitive 
for interpretation of what happens in the space-time (target space) , string models usually 
involve Liouville field, identified as a time-variable in the target-space formalism. (Note 
that for this reason there is usually (at least one) time in the string theory, while space 
can be of any dimension (at least between and 25), not obligatory integer.) Because of 
this extra Liouville field the space-time dimension d usually differs by 1 from the central 
charge of the CFT model, which is coupled to 2<i-gravity to from a string model: d = c+ 1 
and d < 2 is the same as c < 1. 
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where " Van-der-Monde determinant" A(/i) = det^h 3 ^ = YifyAhi — hj) and 
[dU] is Haar measure of integration over unitary matrices. 

The way to derive eq. ( |3.3j ) is to consider the norm of infinitesimal variation 

II 5H || 2 = ES=i I | 2 = Y&^SHijSHa = Ti(5H) 2 = 
= Tr (-W5UWDU + ET+IWtf + W5Duf = 
= Ti(5D) 2 + 2iTr5u[5D, D] + 2Tr (-5uD5uD + (5m) 2 D 2 ) , 

where <5m = j5UW = 5v) and 5D = diag(5/ti . . . 8h^). The second term 
at the r.h.s. vanishes because both D and 5D are diagonal and commute. 
Therefore 

N N 

|| 5H || 2 = E(^) 2 + E (Wij(8u)ji(hi ~ hf- 

i=l i>i=i 

Now it remains to recall the basic relation between the infinitesimal norm 
and the measure: if || 51 || 2 = G a bSl a 5l b then [dl] = y/det a bG a b I\ a dl a , to 
obtain eq.(|3.3|) with Haar measure [dU\ = n« duij being associated with the 
infinitesimal norm 

N N 

|| Su || 2 = Tr(Su) 2 = SuijSuji = | Suij | 2 

i,j=l i,j=l 

and [dUcartan] = U4L1 dun- 

Coming back to the 1-matrix model it remains to note that the "action" 
TrU (H) = YlkLo tk~£^H k with H substituted in the form ( |3.2| ) is independent 
of U: 

N 

TtU(H) =Y / U(h l ). 

8=1 

Thus 

1 N r N 
MO = M II / dhe ^ Ufa - hjf = 

i=l i>j 

1 N r 

= ^J[jdhe u ^A 2 (h), (3.4) 
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provided cn is chosen to be 



where the volume of unitary group in Haar measure is equal to 

Vo W) = rriV • ( 3 - 6 ) 

1 lfc=l K - 



A simple way to derive eq. (|3.6| ) will be described at the end of this section, 
as an example of application of orthogonal polinomials technique. 

3.3 Itzykson-Zuber and Kontsevich integrals 

Let us proceed now to Kontsevich integral, 

F v>n {L} = f dXe -trV(X)+trLX (3J) 
Jnxn 

We shall see shortly that it in fact depends only on the eigenvalues of the 
matrix L (this fact was already used in the previous section), however, this 
time somewhat more sophisticated unitary matrix integrals will be involved. 

Substitute X = U X D X U X ; L = U ] L D L U L in Q and denote U = U X U ] L . 
Then 

^V,n{L} = 

= m=xl d Xi e- v ^A\x) /nxnI - I ga_exp(E™, =1 a; 7 / 5 | U lS | 2 ) . (3.8) 

In order to proceed further we need to evaluate the integral over unitary 
matrices, which appeared at the r.h.s. 

This integral can actually be represented in two different ways: 

I n {X,L}=[ J&L-e*™"* = (3.9) 

Jnxn [aUcartan] 

M e E; fel °M*l£Vl 2 (3.io) 

n [dU Car tan] 
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(the U's in the two integrals are related by the transformation U — > UxUU^ L 
and Haar measure is both left and right invariant). Formula ( |3.9|) implies 
that I n {X, L} satisfies a set of simple equations 



tr 



tr 



d 



dX tT/ 
d \ 



dL 



- trL fc ] I n {X, L} = 0, k>0, 
trX k I I n {X, L} = 0, k > 0, 



(3.11) 



tr 



which by themselves are not very restrictive. However, another formula, 
(|3.10|) , implies that I n {X,L} in fact depends only on the eigenvalues of X 
and L, and for such I n {X, L} = I{x 7 ,ls} eqs.( |3.11| ) become very restrictive 
and allow to determine I{x 1 ,ls} unambigously (at least if I{x 1 ,lg} is 
expandable in a formal power seria in x 7 and Is). The final solution is 



In{X,L} 



n{n — 1) 

(2ir) 2 det 7 5e 



1 l S 



n\ A(x)A(l) 

Normalization constant can be defined by taking L = 0, when 



(3.12) 



I n {X,L = 0} 



Vo\ uin) _ (2tt) n(n 2 1 



(Vol 



and using the fact that 

det jS f 7 (k) 



A(l) 



{h=0} 



fn-\ i 

n 

\k=0 



k\ 



det^sd 5 - 1 / 7 (0) 



Eq.( |3.12 ) is usually refered to as the Itzykson-Zuber formula [60]. In 
mathematical literature it was earlier derived by Kharish- Chandra ]^TJ, and 

20 When acting on /, which depends only on eigenvalues, matrix derivatives turn into: 



tr- 



d 



8X tr 

d 2 



i =E 



7 dx-y ' 



tr ° f - y Jlf + v 1 f g g_\ f. 

<9X 2 ^ 7 S^^^ij^js dxs) 1 ' 



etc. 
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in fact the integral ( |3.9| is the basic example of the coadjoint orbit integrals 
32]- 0, which can be exactly evaluated with the help of the Duistermaat- 



Heckmann theorem JB5|] , ||43|| , [j44|| . This calculation is the simplest exam- 
ple of the very important technique of exact evaluation of non — Gaussian 
unitary- matrix integrals, which is now doing its first steps (see |}6|-[|68|) and 
will be discussed at the end of these notes. 

Now we turn back to the eigenvalue formulation of the GKM. Substitution 
of (ETT3) into (p|) gives: 

n(n — 1) 



Fv,n{L} = (2? ^ m 2 fl / dx 5 e- v ^/\{x)-det lS e x ^ = 

n(n — 1) 

= ft / dx s e- v ^ +Xsls A{x), (3.13) 



5=1 ■ 



where we used antisymmetry of A(x) under permutations of x 7 's in order to 
change ^def^e^ 7 ' 15 for e^s Xsls under the sign of the x$ integration. 

We can now use the fact that A(x) = det^x^ -1 in order to rewrite the 
r.h.s. of ( EU3D : 

JV,n\L\ = (2vr) 2 — -j- , (3.14) 

where 

= J dxx^- 1 e~ v W +lx , 7 > 1. (3.15) 
These functions y3(/) satisfy a simple recurrent relation: 



7-1 

01 \01 

with 



* (3-16) 



$(/) = 0,(1) = J dxe- y{x)+lx . (3.17) 

Note also that if the "zero-time" N is introduced (see subsection 2.6 above 
and f36|), then 

J- V ,n{N | L,\ = J-V{X)~N\ogX,n{L\ = {2tt) * (3.18) 



48 



with just the same y3 7 (7) and 7, S — 1 ... n. If we divide by the quasiclassical 
factor Cy{A}(detA) iV , L = V'(A), in order to transform Kontsevich integral 
into Kontsevich model (see section 2.5), we get: 

7 r ¥T i _ 1 det 7(5 y7 7+ jv(AcO ( o-, Q \ 
Zv{N > T} ~ (deUF ' A(A) • (3 ' 19) 

The role of Cy{A} is to convert (p(l) into the properly normalized expansions 
in the negative iteger powers of A: 



e -AV'(A)+V(A) /^7U 

^(A) = V y 7 (V" (A)) = A^(l + 0(A )), (3.20) 

V27T 

and to change A(/) = A(V'(A)) in the denominator of (|3.18 ) for A(A) in 
( |3.19 ). Instead of the simple recurrent relations ( p. 16 ) for tp the normalized 
functions <p satisfy: 

(p 7 (X) = A(Pi-i(\) = A 7_1 $(A), (3.21) 
where $(A) = y?i(A) an d operator 

1 d 1 V"'(X) 

now depends on the potential V(x). 



3.4 Conventional Multimatrix models 

The multimatrix integrals of the form 

Z N {t^} = 

= C^" 1 / dH^...dH^ IT eS^o^ 5 ^) TT e TrH^Hi^ (3 23) 

can be rewritten in the eigenvalue form using the same Itzykson-Zuber for- 
mula (13TT2]) . Indeed, substituting #( a ) = f/W 1 D^U a) and then defyining 
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jj{a) ^(a+ljt - f/(a) ) wg obtain: 

z N {t^} = 

= w\ n n / dht ] e- y ^ q) )a 2 (/^)) n /^w,^ 1 )} 

i v . 1- i >^ 1 



a=l i=l o=l 
P-l AT , , , p-2 



(3.24) 



inn/ d/4 a) e- v( ^ a)) n ^^a^wja^; 

i V ■ 1 ■ i «/ 1 



iv • a =l i=l o=l 



,(a),(a+l) 

where the same trick is done with the substitution of -pdetjje i * for 

EAT ,(a),(a+l) ( a \ 
*=i 4 * under the sign of the /i- -integration (step by step: first - 

for a = 1, then for a = 2 and so on). Note that all the Van-der- Monde 
determinants disappeared from the final formula at the r.h.s. of eq. (|3.24|) , 
except for those at the ends of the matrix chain (at a = 1 and a = p — 1). 

If the chain was closed rather than open, i.e. there was an additional 
factor of e Tr ' ff(p 1>H<1) under the integral in (|3.23|) , then the trick with separa- 
tion of all angular-variable (unitary-matrix) integrations would not work so 
simply: in addition to the Itzykson-Zuber integral the much more involved 
quantities would be required, like 



In\Xi, X 2 \ L} = 

[dUi] [dU 2 ] 

<n [dUi t Cartan\ [dU 2 £artan\ 

exp (irXxUxLUl + XxX 2 U 2 LU\ + trX 1 (U 1 Ul)X 2 (U 2 Ul] 



(3.25) 



This (so far unresolved) closed chain model (lattice Potts model) is an exam- 
ple of non-eigenvalue models, in the p = oo case it turns into "compactified" 
c = 1 model. This theory is more complicated then the so far simplest class 
of non-eigenvalue models of "induces Yang-Mills theory", known as Kazakov- 
Migdal models. 



3.5 Determinant formulas for eigenvalue models 

We are now prepared to make the crucial step towards understanding of 
mathematical structure behind eigenvalue models, which distinguishes their 
partition functions in the entire variety of arbitrary iV-fold integrals. This 
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structure expresses itself in the form of determinantal formulas, which we are 
now going to discuss. In the next section 4 these formulas will be identified 
as examples of r-functions of KP and Toda hierarchies. 

Looking at the relevant integrals Q3.4|) , ( [3.24 ) one can notice that inte- 



grals over different eigenvalues with non-trivial measures which depend on 
the shape of potentials U or V, are almost separated, the only "interac- 
tion" between different eigenvalues being defined by universal (potential- 
independent) quantities, made from the Van-der-Monde determinants. This 
feature is intimately related both to its origin (decoupling of angular variables 
in original matrix integral) and to its most important implication (integra- 
bility). The main property of the Van-der-Monde determinant is that it is 
at the same time a Pfaffian (and it is in this quality that it arises from 
matrix integrals) and a determinant (and this is the feature that implies 
integrability) : 

J[(hi - hj) = A(h) = det^/ir 1 . (3.26) 

i>j 



We already used this property above, when going from eq. ( |3.13|) to eq. (|3.14|) 



which as we shall see later is the crucial step in the proof of integrability of 
Kontsevich model. In that case determinantal formula (|3.14|) for partition 



function was trivial to derive, because the integrand was linear in Van-der- 
Monde determinants. Now we turn to slightly more complicated situations, 
involving products of Van-der-Monde determinants. 
Consider an eigenvalue model of the form: 

z N=^I[j dfi hk - hk A(h)A(h), (3.27) 

' k=l 

to be refered to as "scalar-product" model. All conventional multimatrix 
models (|3.23|) belong to this class. In the case of the 1-matrix model (|3.4|) 



dn h>h = dhdhe u{h) 5(h - h) , (3.28) 
while for conventional multimatrix models (|3.24j) 



P~2 , p-l p-2 



dA*Ad), h (p-i) = dhWdhW II / dh {a) II e Uaihia)) II e h(a)h(Q+1) (3.29) 



a=2 a=l a=l 
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If dfihfi = 5{h — h)dhd[ih we call this measure local. The main feature of 
local measure is that operator of multiplication by H (or any function of h) 
is Hermitean. Thus measure is local in the 1-matrix model, but is non-local 
for all p — 1 > 1. In the latter case the measure is defined to depend only on 
h = h^ 1 ' and h = h^' 1 ', all other hS a \ a = 2 . . .p — 2 being integrated out, 
what makes the "interaction" between h and h more complicated than just 
S(h — h) in the one-matrix (p = 2) and e hh in the two- matrix (p = 3) cases. 
In no sense the set of particular formulas ( |3.29| ) for p > 3 is distinguished 
among other scalar-product models, and from now on we shall not consider 
conventional multimatrix models with p—1 > 2 as a separate class of theories. 
Eqs. (|3~27l ) and (|3~26D imply together that 

1 N f -•_] 
Z * = ^ II / d^rh^tikh^Detjkhi 1 = 

k=l 

= DeUj J d/i^h^h^ 1 = Delicti" 1 \ h j ~ x ), (3.30) 
where an obvious notation is introduced for the scalar product 

(f(h) | g(h)} = f dn h - h f(h)g(h). 



We can now be a little more specific and introduce time- variables and 
ik so, that 

oo oo 

U{h) = Yl hh k , U(h) = £ hh\ (3.31) 

k=— oo fc=— oo 

and d[i h ^ is already independent of h and h. If we now denote H f (t, t) 
{ I | 1), then 



and 



Det N HL (3.33) 
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where Detjv stands for determinant of the N x N matrix Tii-ij-i (which 
is defined itself for any integers with i,j = 0, . . . , N — 1. Characteristic 
property of TC{j is its peculiar time-dependence: 

dh = Hi+k >i ] dh = H ^ +k - (3 ' 34) 

Eq. ( |3.33| ) provides the determinantal formula for all scalar-product mod- 
els. The case of the local measure - for 1-matrix model - is a little spe- 
cial. In this case U(h) contains the whole information about the measure: 
d^hh = ~~ h)dfih, dfih = e u ^dh, and there is no U(h) (or t simply 
coincide with t). Then ( |3.33| ) is still valid, but 

Hli = (h l | h j ) = (h l+j ) = -^—H f {t) = 
3 X ' ' dfi ht% X 1 d^ h dt i+j v ' 

d ^ 1+3 



if i,j>0 



- H' { t). (3.35) 



The same formula ( |3.35| ) can be also derived as a limit of eq.( |3.14 ) for Kont- 
sevich integral. Indeed, 

Z N {t} = c N f dHe™ {H) = lim T UN {L} = 
JNxN l^o 

= ^o ^'C^ = ^^F^^ = ^«U S -, (3-36) 



where this time 

«*-«-. g= S g ^ C = 0) ' = 4 ' (I)"" 2 *'" 1 !'=» • (3.37) 



Now we note, that the action of 4 on ^ u ^(l) = j dxe u ^ +lx is equivalent to 
that of f^")> since this is no longer a matrix integral, and thus 



d x 1+3 



i.e. H f (t) = $^>(0). 



w£=(^-l ^ } (°)» ( 3 - 38 ) 
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Con formal multimatrix models were introduced in section 2.3 above just 
as eigenvalue models. For the A v _\ series partition functions are defined to 
be 

^ JN a xN a 

p-2 (3 39) 

• J] Det H {a+l) ) = 

P-1 1 N a , , , p-2 

= n ttt n / dh^e^^A^) n ik^ - 4 a+1) )- 

This expression does not have a form of eq.( |3.27|) , thus conformal matrix 
models for p— 1 > 1 are not of the "scalar- product" type. We shall sometimes 
call them (p — l)-component models, because they are related to the multi- 
component integrable hierarchies. The simplest way to proceed with their 
investigation is to use the same trick with Kontsevich integral, which was 
just applied in the 1-matrix case. 

Let us start from a very general (p — l)-component model: 

Z = TT / dHWe TtU «^K(HW . . . H^-V). (3.40) 

^JN a xN a 



It can be also represented in terms of Kontsevich integrals: 



d d 



p-i 



J tr ULl tr f «=1 



(3.41) 



L( a )=0 



This representation is not very usefull, since the limit L — > is not easy 
to take, unless if is a polinomial in the eigenvalues of all its arguments. 
However, this is exactly the case for our confromal models ( |3.39| ). Indeed, 

^ ,= iH^ 0/ - /0 5iM' (3 - 42) 
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Still this is not very convenient, because representation ( 5711 ) for T contain 
A(L) in denominator, which are not very pleasant to differentiate. Simpli- 
fication can be achieved if instead we rewrite the original expression at the 
r.h.s. of ( |3.39|) as follows: 



- a fd\ n p-2 A (e a \ A ( a 



X 



(3.43) 



where A(/i, />') ee n£#i"M II&K -/*',) iCi^-^)- This formula 
already takes the specific form of .fT into account. The product of integrals 
in brackets at the r.h.s. of (|3.43|) is equal (for every fixed a) to 



N. 



1 N <* 



(3.44) 



3=1 



(compare with eq.( 



In order to simplify the notation we shall further denote $^ c/q ^(/) = 



J dxe Ua ^ +lx through $ a (/), and 
- through d k $ a {l). Thus 



2(")=0 



(3.45) 



If p — 1 = 1, differential operator is just a square of determinant A(d/dl) 
and we can use the relation 



Ep Det 



*? rt p(j) - 








1 

hp(\) 
h 2 

a p(i) 


hp(2) 
h 2 

a P{2) 
^P(2) 


h 2 

^P(3) • • 
'<P(3) • • 


i.JVi-1 "I 

■ a P(N 1 ) 

h Nl 

■ a P(iVi) 

h N 1+ l 

■ "P(iVi) 


h m-i 
_ a p(i) 


h Ni 
n P(2) 


' 4 P(3) • • 


h 2N!-2 
■ Ai P(JVi) 



(3.46) 



55 



where the sum is over all the N\ permutations P of N elements 1 . . . N, 
in order to conclude that (|3.45| ) reproduces our old formula (|3.33|) , (|3.38|) : 

Z N = Detijd^- 2 ®. 
c p — 1 = 2 w< 

A(h)A(h,h')A(h' 



N — ^hj*- 

For p — 1 = 2 we need to use a more complicated analogue of Q3.46| ) : 



1 

hp(i) 
h 2 

a p(i) 


h P {2) ■ ■ 
h 2 

h 3 


1 -1 

n 'P{N 1 ) 

n P{N l ) 
h N l+ l 


1 

hp(i) 
h 2 - 


hp(2) ■ ■ 

h 2 - 

h 3 - 


02-1 

h N - 2 

1.N2+I 
n P(N 2 ) 


. n P{i) 


0i 


iM+Nr-2 
■ P(jVi) 
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where M = E„=i N a . Making use of this formula, we conclude that the r.h.s. 



of (|3.45|) for p — 1 = 2 is also representable in the form of determinant: 



Det 



<9$ <9 2 <l... d Nl ® dl> 



9$ . . . d"*- 1 ® 

<9 2 4... 9^4 



(3.47) 



here $ = $1, $ = $2 and all arguments l^ a ' = 0. It is especially easy to 
check formula (|3.47|) in the simplest case of Nx = N2 = 1. Then it just 

1 1 . Analogues expressions for p 



says that h — h = Det 



h h 



1 > 2 are 



more involved, they are no longer just determinants: this is obvious already 
from consideration of the simplest case of N\ = . . . — iV p _i = 1, the product 
[]a=i(/i' ' — h^ a+1 ^>) is no longer determinant of any nice matrix. 



3.6 Orthogonal polinomials 

Formalism of orthogonal polinomials was intensively used at the early days 
of the theory of matrix models. It is applicable to scalar-product eigenvalue 
models and allows to further transform (diagonalize) the remaining determi- 
nants into products. In variance with both reduction from original iV 2 -fold 
matrix integrals to the eigenvalue problem, which (when possible) reflects 
a physical phenomenon - decoupling of angular (unitary-matrix) degrees of 
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freedom (associated with <i-dimensional gauge bosons), - and with occurence 
of determinant formulas which reflects integrability of the model, orthogonal 
polinomials appear more as a technical device. Essentially orthogonal poli- 
nomials are necessary if wants to explicitly separate dependence on the the 
size N of the matrix in the matrix integral ("zero-time") from dependencies 
on all other time- variables and to explicitly construct variables, which satisfy 
Toda-like equations. However, modern description of integrable hierarchies 
in terms of r-functions does not require explicit separation of the zero-time 
and treats it more or less on the equal fooring with all other variables, thus 
making the use of orthogonal polinomials unnecessary. Still this technique 
remains in the arsenal of the matrix model theory^] and we now briefly ex- 
plain what it is about. At the end of this section two simple applications 
will be also described: one to evaluation of the volume of the unitary group, 
another - to direct proof of equivalence of the ordinary 1-matrix model and 
the Gaussian Kontsevich model. Both these examples make use of explicitly 
known orthogonal Hermite polinomials and in this sense are not quite rep- 
resentative: usually orthogonal polinomials are not known explicitly. Some 
applications of such "abstract" theory of orthogonal polinomials to the study 
of matrix models will be mentioned in the following sections. 

In the context of the theory of scalar-product matrix models orthogonal 
polinomials naturally arise when one notes that after partition functions 
appears in a simple determinantal form of eq.( |3.30|) , any linear change of 



basises h % — > Qi(h) = J2k^-ikhk, h? — > Qj(h) = J2iBjih l can be easily 
performed and Z — > Z -detA-detB. In particular, if A and B are triangular 
with units at diagonals, their determinants are just unities and Z does not 
change at all. This freedom is, however, enough, to diagonalize the scalar 
product and choose polinomials Qi and Qj so that 

(Qi(h) | Qjih)) = e%. (3.48) 

Qi and Qj defined in this way up to normalization are called orthogonal 
polinomials. (Note that Q does not need to be a complex conjugate of Q: 
"bar" does not mean complex conjugation.) Because of above restriction on 



21 Of course, one can also use this link just with the aim to put the rich and beau- 
tifull mathematical theory of orthogonal polinomials into the general context of string 
theory. Among interesting problems here is the matrix-model description of g-orthogonal 
polinomials. 



57 



the form of matrices A and B these polinomials are normalized so that 

Q i (h) = h i + ...', Q j (h)=h? + ... 

i.e. the leading power enters with the unit coefficient. From ([T 
it follows that 

N 

Z N = l[e**-K (3.49) 

8=1 

This formula is essentially the main outcome of orthogonal polinomials theory 
fro matrix models: it provides complete separation of the iV-dependence of Z 
(on the size of the matrix) from that on all other parameters (which specify 
the shape of potential, i.e. the measure dfi h ^), this information is encoded in 
a rather complicated fashion in 0j. As was already mentioned, any feature of 
matrix model can be examined already at the level of eq. fl3.30p , which does 
not refer to orthogonal polinomials and thus they are not really relevant for 
the subject. 

We can, however, reverse the problem and ask, what can matrix models 
provide for the theory of orthogonal polinomials.0 The first question to ask 
in the theory of orthogonal polinomials is: 

Given the measure d\i h ^ what are the corresponding orthogonal polino- 
mials? 

Usually the answer to this type of questions is not at all straightforward. 
Its complexity, however, depends on what one agrees to accept as a suit- 
able answer. Of particular interest for our purposes below would be integral 
representations. It would be very helpfull to have just an integral trans- 
formation, converting the set of orthogonal polinomials for given d/i^h m t° 
some standard set, like Qf = x % . Unfortunately, such transformation is 
rarely available, though there are important examples: classical orthogonal 
polinomials and their g-analogues (expressed through the (g-)hypergeometric 
functions, which usually possess integral representation of a simple form, see 
[p5[| for an introductory review of such integral formulas, which are in fact 

22 0f course, we can hardly get anything new for that theory, but the purpose is to 
see, which features are immediate consequences of the "physically-inspired" approach. As 
usually this can help to somehow organize the existing knowledge in appropriate system. 
This is, however, not our goal in these notes: only a very simple example will be mentioned, 
which will be also of use in our futher considerations. 



) and (143) 
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well known in CFT). The simplest example of this kind, which will be used 
below is the set of Hermite polinomials: 

Re k (h) = -^e^ [ (ix) k e-^- ixh dx = (h - —) k ■ 1 = 
\/2n J ah 



III d k _h? 1 h£ (I _ ;,~ 



= e T(-_). e -T = 5ieT ( k -2 ;j5 )V- T =»* + ..., (3.50) 

orthogonal with the local measure dfih — e 2 . 

For generic measure the answer of this type does not exist in any universal 
form. However, matrix models still provide a somewhat peculiar integral 
representation for any measure, with the number of integrations depending 
on the number of polinomial. In order to obtain this expression, let us 
consider a slight generalization of formula (|3.27 ) 

I N 

Zn{\} = j^UJ dfi hkrhk A(h)A(h) [](A 7 - h k ). (3.51) 

k=l fc,7 

Then A(h) nfc, 7 (A 7 — hk) = A(h, A)/A(A), and A 7 can be just considered as 
h,N+-y, which are not integrated over in ( |3.51| ). Then it is clear that 

A(M)=Det( ^ifo.) ^+7-1(^)^1 (3>52) 

\ Qi-l(Xs) Qn+~/-1{X6) J 

while A(h) = DetjkQj-i(h~k)- Since all the QN+~j-i{hk) ar e orthogonal to all 
Qj-i(hk) (because N + 7 — 1 7^ j — 1), we obtain: 

det 7 5<5jv+ 7 -i(A5) , , 

^N\^S\ = ~A~7~\\ {6.06) 

l\{\) 



In particular, 



Qn(X) = (3.54) 



where both the numerator and denominator can be represented by iV x N- 
matrix integrals. 

Inverse "main question" of the theory of orthogonal polinomials is: 
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Given a set of polinomials 

Qi(h) =hi + 
Qj(h) = hj + 

what is the measure dfi h ^ w.r.to which they form an orthogonal system? 

We shall not discuss the complete answer to this question and consider 
only the case of the local measure, when Qi = Qi. Then usually the answer 
does not exist at all: not every system of polinomials is orthogonal w.r.to 
some local measure. It is easy to find the necessary (and in fact sufficient) 
condition. As was mentioned above, the local measure is distinguished by the 
property that multiplication by (any function of) h is Hermitean operator: 

(hf(h) | g(h)) = (f(h) | hg(h)), if dii h - h ~ 5(h - h). (3.55) 

This property implies, that the coefficients Cjj in the recurrent relation 

i 

hQ t (h) = Q i+1 {h) + £ CijQjih) (3.56) 

3=0 

are almost all vanishing. Indeed: for j < i 



c 'j 



(hQi{h) 


1 Qj{h)) 


(Qj{h) I 


QAh)) 



(Qi(h) 


1 Qi(h)) 


(Qj(h) 


1 Qi(h)) 



(Qi(h) 


hQ 3 {h)) 


(Qj{h) 




- = 





(3.57) 



In other words, polinomials, orthogonal w.r.to a local measure are obliged to 
satisfy the " 3-term recurrent relation" : 

hQi{h) = Q i+ i(h) + CiQi(h) + RiQi-i(h) (3.58) 

(the coefficient in front of Qi+i can be of course changed by the change of 
normalization). Parameter Cj vanishes if the measure is even (symmetric 
under the change h — > — h), then polinomials are split into two orthogonal 
subsets: even and odd in h. Partition function ( |3.49| ) of the one-component 
model can be expressed through parameters Ri = e^ 1- ^ -1 of the 3-term 
relation: 

N-l 

Z N = Z 1 1 [[ Rf~\ (3.59) 
i=i 
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thus denning a one-component matrix model (i.e. particular shape of poten- 
tial), associated with any system of orthogonal polinomials. 

Our "inverse main question" in the case of the local measure should be 
now formulated as follows: Given a set of orthogonal polinomials Qi(h) = 
h l + . . . which satisfy the 3-term relation fl3.58Q , what is the measure &\x\?. 

As every complete orthogonal system of functions, orthogonal polinomials 
satisfy the completeness relation: 

oo 

Y.^QiQ^Qiih) = 6^{h,h), (3.60) 

8=0 

where 5-function, associated with the measure d^h is defined so that 

f(h)8^(h,h')dix h - h = f(h') (3.61) 



for any function f(h). Since for the local measure dfih = e u ^dh the 5- 
function is just 8^(h, h) = e- u ^5(h - h), as an answer to our question we 
can take a representation of U(h) in terms of the corresponding orthogonal 
polinomials: 

e- u ^5(h ~h) = ± 90®AV . (3.62) 

k=Q \Qk | Qk) 

As usually this relation should be understood as analytical continuation. The 
squared norms || Qk 1 1 2 in denominator are expressed through the coefficients 
Ri of the 3-term relation up to an overall constant: || Qk || 2 = lTi=i R% I 

Qoll 2 - 

For example, in the case of Hermite polinomials ( |3.50|) we have: 



Ee k+1 (h) = (h- J^)He fc (/i) = KKe k (h) - J^He fc (/i) = 

= KRe k {h) - kRe k -i(h) (3.63) 

(the last equality holds because and h — J- play the role of annihilation 
and creation operators respectively). This means that the 3-term relation 
is satisfied with R k = k and thus II Het || 2 =|| Hen II 2 k\ We shall use the 
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normalization condition || Heo || 2 = \^2n. Then for e u( - h ^ we get: 

e- u{h) 5(h -h) = 

» He fc (fr)He fc (fr) 1^1,. d fc , d fc 
^ 11 JTe fc || 2 v^F^A;! 1 dh } 1 d/i J 



fc=0 II II V ^" fc=0 

1 ■' ^ 00 1 , d 2 , , h 2 h 2 



=e 2 



27r fc=0 



— ( = ) e 2 2 = e 2 2 edhdhe 2 

k=n kVdhdh' v/2tt 



1 /■ /■ dada n 2 s 2 a jL +5l A. _i£_h£ 
zlm / / e e 2 ^ 2 e dft^ <&e 2 2 



V2^ J J 2vr 

1 l m / f dad(X c ~ I ("+") 2 c ~ I («+«) (fe+fc) c - \ (a-S) (h-h) 



2vr J J 2ti 

= e 2 — h). 

3.7 Scalar-product models in Miwa parametrization 

We shall now make the first step towards clarification of the interrelation 
between the scalar-product and Kontsevich models. We already know that 
in the latter case the important role is played by representation of time- 
variables in the form of 

T k = itrA" fc , (3.64) 
with nxn matrix A, which will be further refered to as Miwa parametrization 



(expressions of some similar form were first introduced in [70]). Let us now 
perfrom such transformation in the case of the scalar-product model. Let us 
use eq.( |3.31|) to define the time-dependence of the measure, only ignore the 
t-variables. Namely, introduce d[x h ^ = e u ^dP h j l (i.e. dv h j l = e u ^dfi h ^). 
Substitute 

t k = T (j trA ~ fc + r k) ( 3 - 65 ) 



and obtain: 

e u(h) = e-vw^tr^mY = e ^(")^ 



detA , , 



detA 



e -v(h) _u_ (3-66) 
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where V(h) = ±J2k r kh k . Let us choose upper signs in these formulas. Then 
we can use eqs. fl3.5l|) and ( |3.53|) to conclude that in Miwa parametrization 



7 {dn} _ 1 rjjdu} r ^ x _ ^{du} det 7< 5Qj V+7 _ 1 (A^ 



(detA)" JV ° J JV A(A)(detA) J 

where dv^ = c~ v!kK) dv^ and Qk are the corresponding orthogonal polino- 
mials. In other words, we reduced the model with potential U(h) to another 
model with potential, —V(h), and expressed the difference in terms of or- 
thogonal polinomials Qk'- 

_ 1 det 7 aQjy +7 _i(A^) , . 

Z { f } ~ (detA) N ' A(A) ' 1 ' 

If V(h) is adjusted to give rise to some simple orthogonal polinomials, (i.e. 
if the new model zjf u ^ is easy to solve), this representation can considerably 
simplify the original model. 

Another interpretation of this formula is that we obtained a GKM-like 
representation of the from of ( p.!9|) for the discrete scalar-product model. 
The only difference is that (pjY^ in ( p. 19 ) are changed for Q 7 -i in (|3.67|) . 



This is an important difference, because in GKM are defined to by 

integral formulas like ( |3.15| ), = ((x 7-1 )) or, alternatively, satisfy the 
recursive relations like Q3.21|) . Moreover, generic (py^ are infinite formal 



series in A -1 , while Q 7 _i are orthogonal polinomials. This discreapancy is one 
of important stimuli for further developement of the concept of Generalized 
Kontsevich model, as well as for search for convenient integral representations 
for orthogonal polinomials. 

There is, however, at least one interesting situation when the two formulas 
indeed coincide. This is the case of Gaussian potentials V and V, when both 
and Q-y-i are represented by orthogonal Hermite polinomials, which 
possess integral representation, exactly adequate in the context of GKM. 
This is the subject of our consideration in the next subsection. 



3.8 Equivalence of the discrete 1-matrix and Gaussian 
Kontsevich models 

Let us take the ordinary 1-matrix model with the local measure d\Xh = 
e u ^dh to be the scalar-product model, considered in the previous subsection 
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and take Miwa parametrization with upper signs and with r k = — \8k,2 (as 
we did in the section 2.6). Then V(h) = J2k r kh k = — = The rel- 

evant orthogonal polinomials Q are just Hermite polynomials of imaginary 

h 2 

argument:^ Q k 2 = i~ k ~Kek(ih) = h k + . . .. These polinomials possess 



an integral representation (|3.50| ) 



i 1 - fe He ifc _ 1 (z/i) = -Ue-^ / x k - x e~^+ xh dx ( = 0) ^\h). (3.68) 

V 27T J 



Using (|3.67|) and ( p,19| ) we obtain a remarkable relation between the two 



matrix models: 

Z N {t = 0;t k = -itrA- fc + ±5 fcj2 } J NxN dHe^T=o^ Hk 



Z N {tk - |4,2} f NxN dHei H2 

f dX(detX) N e~ tr2 ^ +AX = Z^{N,t}, (3.69) 

N Jnxn 2 



_ tr — 

e tr 2 



(27r)^(detA) 



where Z N {t k = |<5/c,2} = (— 27r)^ c^. This relation can be also regarded as 
an identity 



f NxN dHe2 TrH Det(A ® / - / ® H) 

Li|A , U . (3.70) 



J NxN dHe^ 
I nxn dXe-^ x2 det N (X + A) 



Inxn dXe ~ 



2 



valid for any A. Note that integrals are of differents sizes: AT x AT at the 
l.h.s. and tixnat the r.h.s. While A^-dependence is explicit at both sides of 
the equation, the n-dependence at the l.h.s. enters only implicitly: through 
the allowed domain of variation of variables t k = — ^trA~ fc + \5k,2- (This can 
serve as an illustration to the general statement that the shape of Kontsevich 

23 Note that this system of functions (p k = i~ k Hek(ih) looks like ipo = 1, ipx = h, if2 — 
h 2 + 1, . . ., and does not resemble any set of orthogonal polinomials with a local measure 
(for example the product <po-tp2 = h 2 + l may seem positive definite, this being inconsistent 
with orthogonality requirement (ipo | tp 2 ) =0). The thing is that integration at the l.h.s. 



of eq.(3.69) is well defined only along the imaginary axis, while integrals along the real 



axis are understood as analytical continuation. 
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partition function Zy, considered as a function of T's rather than L or A, is 
independent of the matrix size n.) Identity ( 3.69Q was anticipated from the 
study of Ward identites for the Gaussian Kontsevich model in []5B] (see eq. 
(2.53) in the section 2.6 above), and it was derived in the present form in 
ref.||. 

Eq. ( |3.69| ) can be used to perfrom analytical continuation in N and define 
what is Zjy for N, which are not positive integers. Since = for all 
negative integers (see eq.( |3.77| ) below), the same is true for Z^. In the next 
section 4 we shall see that it is characteristic property of r-functions of forced 
hierarchies. 



3.9 Volume of unitary group 

Formalism of orthogonal polinomials provides also a simple derivation of 
eq.( |3.6| ) for the volume of unitary group. Consider eq.( |3.4|) with U(H) = H 2 . 



Then Gaussian matrix integral can be easily evaluated: 



c N I dHe 2 

INxN 



n N 
c N Ylj dHue-^n d 2 H l3 e-\ H ^ = (2tt 

i=l i<j 



N 2 /2 



while according to eqs.(|3.48|) and (|3.49|) the same integral is equal to 



1 N r+oo N N 

jnU e-Hn(^-^) 2 =ni|He^ 

7V! i=l J -°° i>j 3=1 



Here 1 1 Hej_i 1 1 stand for the norms of orthogonal Hermite polinomials (|3.50|) , 
l|He fc || 2 = _ 
get: 



2nk\. Comparing the two expressions for the same integral we 



N -n- 1 



'N 



1 = (2-)- n 



k=0 



(27T) 



JV(JV-l) 



2nk\ 



ntL-o 1 k\ 



(3.71) 



According to 



-N 



Vol 



U(N) 



(Vol^d)) 



N 
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and Vol[/(i) = 2n. Thus we obtain eq.(|3[ 



Volc " K) = "nE^- 

An example of somewhat more sophisticated (quantum) group-theoretical 
quantity, arising from Gaussian matrix models, is provided by the following 
formula for the q- factorial [FT] (see also |72| ): 

1 -TT 1 SS NxN dH[dV]e-<»™+™™* 
fa Q)n ~ L\ 1 - t Vol u(N) J NxN dHe~™ 2 ™ 2 ■ 

Integral in the numerator is over Hermitean (H) and unitary (U) N x N 
matrices, and q = m 2 — yrnF—1. 

Explicit expression ( |3.71| ) can be used to prove that cn = for all negative 
integer iV |3^| . Eq. (|3.71|) defines cn only for positive integer N as a finite 
product. There is an obvious prescription for analytical continuation of such 
products, provided continuation of the items is known (it can be considered 
as implied by the similar formula for integrals with the varying upper limit): 
Let 

N 

F(N)= £ f(k). (3.73) 

k=~ oo 

Then 

N 

S(N) = J2f(k)^F(N)-m (3-74) 

k=l 

and, obviously, F(0) - F(-N) = £Li-jv f( k )i so that 

JV-l 

S (- N ) = F (-N) - F(0) /(-*)■ (3.75) 

fc=0 

Exponentiation of this formula gives the rule for the products. In the case 
of cat one can treat factorials in ( |3.71D as Gamma-functions, 

N 

{2tx) E ^ R c n = nr(fc), (3.76) 

fc=i 
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and obtain: 

/N-l 



{2n) -mm c _ N = I fir(-k)) = o, (3.77) 

\fc=0 / 



because of the poles of T-functions. 



4 Integrable structure of eigenvalue models 
4.1 The concept of integrability 

Integrable structure of dynamical system implies that all the dynamical char- 
acteristics - solutions of equations of motion for a classical system and func- 
tional integrals for a quantum one - can be found exactly. According to this 
description the notion of integrability is not very concrete, and in fact it 
evolves with time, including more and more classes of theories into the class 
of integable systems. Nowadays we consider the following types of theories 
as clearly belonging to this class: 

- Free motion (classical or quantum) on group manifolds and homogeneous 
spaces; 

- 2-dimensional conformal theories and their "integrable massive defor- 
mations" ; 

- Integrable hierarchies of the (multicomponent) KP and Toda type and 
their reductions; 

- Functional integrals, subjected to conditions of (generalized) Duistermaat- 
Heckman theorem; 

- (Eigenvalue) matrix models; 

- Topological theories; 

- Many supersymmetric models (at least those allowing for Nicolai trans- 
formation and/or Duistermaat-Heckman-like description); 

- Systems with (infinitely) many local integrals of motion. 

This list (nothing to say about the order of items) is rather arbitrary. Also 
different items are not really different and (as it should be) can be considered 
as different descriptions of the same reality. Now we discuss very briefly at 
least some of the most important views on the concept of integrability. 



67 



Often the notion of integrability is related to occurence of " many enough" 
integrals of motion ("many enough" means equal to the number of degrees 
of freedom). This is, however, not such a rigid definition as one can think. 
In fact, in classical mechanics there is usually a complete set of integrals of 
motion available: just initial conditions in the phase space (or, to be more 
sophisticated, angle- action variables). The problem is, however, that 

a) these obvious integrals are very complicated (non-local and multi- 
valued) functionals of the current coordinates, and 

b) in general situation they are very "unstable" under a small change of 
current coordinates ("divergency of trajectories"). 

In order to avoid these problems one usually imposes a "locality" con- 
dition on equations of motion. While this is a reasonable thing to do for 
particular classes of theories (e.g. possessing a well defined kinetic term, 
which is quadratic in momenta), this is not a nice decription in general situ- 
ation, since "locality" is not invariant under arbitrary (including non-local) 
change of variables. In practice, when approached from this side, integrabil- 
ity implies a kind of "regular" behaviour of trajectories and some more or 
less nicel defined transformation from "natural" (or, better to say, "original") 
coordinates to the action-angle variables. 

Situation becomes even less clear when quantum theory is considered, 
since " chaotic behaviour" no longer implies anything really " chaotic" for the 
quantum system. Again, very much depends on what kind of observables 
one wants to consider, and any notion of "regularity" is not enough under 
arbitrary change of variables. 

This can be made even more transparent, if one recalls the idea of univer- 
sality classes, so important in the modern theory. The idea is that even in the 
cases when behaviour of the system seems absolutely haotic from any naive 
point of view (like in the cases of turbulence or quantum gravity), one can 
and should introduce new variables (which can be very complicated functions 
of original ones), which have smooth and well defined correlation functions. 
In most cases one is not attempting to find a complete set of such variables 
(and thus some information is lost), but this reflects nothing but the current 
state of knowledge, and in fact in studies of 2d quantum gravity the goal of 
complete descrption is already clearly fromulated. 

Despite these comments, the "definition" of integrability in terms of 
"many enough" local integrals of motion should be put at the first place 
in our discussion, because most of the systems which were so far considered 
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as integrable more or less naturally get into this class, allowing for some 
prefered choice of dynamical variables ("more or less" appears because some 
"minor" non- locality is usually present in any interesting examples, where 
angle-action variables are not obvious from the very beginning). 

This " definition" is so unclear because we attempted to look for a generic 
description of integrability. Most interesting approaches, however, go from 
another direction. One starts from some simple system and then perform a 
change of variables, which makes it looking much more complicated (being 
still simple in its essense). This appears to be a much more fruitfull view on 
the problem and in fact all the other items of our list above are describable 
in this kind of terms. 

A trivial, but surprisingly representative example of this approach is pro- 
vided by a free particle, moving in flat D-dimensional space. The eigen- 
functions of Laplace operator are just plain waves or, equivalently, spherical 
harmonics. The radial part of the j-th harmonic is already a not very simple 
function, satisfying the equation 



d 2 D - 1 d C 2 {jY 



dr 2 r dr r 2 



ip(r) = Eip(r). (4.1) 



This equation is of course less trivial than the original Laplace equation, but 
their solutions are related in a simple way. In order to find a solution of ( |4.1|) , 
say, for j = 0, one should just take an angular average of a plane wave: 



Mr) = y e^t^-V; | V \= 1. (4.2) 

This integral representation expresses the solutions of QL1|) through Bessel 
functions, and this is in fact the proper way to derive the well-known formula: 

Mr) = (y) J^kr). (4.3) 

If one expands the exponent in the integral in a series, the standard expansion 
for the Bessel function arises. 

A slightly more involved example is the quantum mechanical model of a 
particle in the potential e~ q , i.e. the theory of equation 

1 ' h + e"« ) i(>(q) = (4.4) 



dq 



2 



69 



(one of course recognizes a simplified version of Toda models). It can be 
solved by projection of the simple Sshrodinger equation for a particle, moving 
on the upper part of the hyperboloid Xq — x\ — x\ = 1; x > J73[] . If 

x = cosh | + -z 2 e q ^ 2 ; x\ = sinh ^ — -z 2 e q ^ 2 ; x<i = ze q ^ 2 , 



then q = log(xo + #1), Laplace operator on hyperboloid is 

d 2 Id 1 _ q d 2 
L= 8ql~2dq~ + A eq d^ ^ 

and average of the wave function ip\(q, z) provides the following expression 
for solutions of (|4.4j) : 

^Jq) = e lXq / t 2lX - l e- {t+eq ^dt. (4.6) 
Jo 



This idea, which is sometimes refered to as "projection method" (see |73 
for a broad review) , reveals hidden symmetries of some complicated systems 
(which do not possess any symmetry at all in the usual, Noether-like, sence 
of the word), by considering them as embedded into wider theories with 
more degrees of freedom. Quantum mechanical examples of applicability of 
the method are by no means exhausted by the two systems above, one can 
consider various projections, starting from (exactly solvable problem of) the 
free motion on any group manifold, and in general this gives rise to a very 
important theory of "zonal spherical functions", which nowadays is increas- 
ingly attracting attention because of its obvious links to integrability theory 
and quantum geometry (see [?4[] for discussion of the latter relation and [|75| 



where also relations with orthogonal polinomials and Generalized Kontsevich 
model are partly revealed). An extremely important example of free motion 
on a group manifold (in the infinite-dimensional - Kac-Moody - case) is pro- 
vided by 2-dimensional WZNW (Wess-Zumino-Novikov-Witten) model and 
the corresponding version of projection method is known as Hamiltonian re- 
ductions in conformal field theory. Again the resulting theories (like minimal 
conformal models) do not possess any kind of symmetry in the usual sense 
of the word, but still they are very simple and exactly solvable, remembering 
their origin in the theory of free fields. 

In principle the theory which is reduced, i.e. complemented by constraints 
(initial conditions), does not need to be absolutely symmetric, i.e. to have 



70 



Casimir operator or even zero (as in the WZNW case) as its Hamiltonian. 
It can be in fact possible to use the projection method to gain a lot of infor- 
mation about reductions of theories with more sophisticated Hamiltonians, 
which are non-trivial functions of group generators. The simplest exam- 
ple is provided by the theory of quantum-mechanical " quasi-exactly-solvable 
models" [|76],|77j and its CFT-generalizations [f7j,|78[]. A more elaborated 



technique has the name of "localization theory" ,0 (known also as geometrical 
quantization, Fourier analysis on group manifolds and Duistermaat-Heckman 
theory), it provides a very wide generalization of the above averaging proce- 
dure, which maped plain waves into Bessel functions. The classical sample 
system to illustrate all the aspects of integrability, starting from free motion 
and ending with anyonic statistics, VFoo-algebras and 2d Yang-Mills theory, 
is the Calogero-Sutherland system, which can be associated in a uniform 
way with any simple Lie algebras and in an "intermediately involved" form 
looks like a multiparticle theory in 1+1 dimension with interaction potential 
g 2 sin -2 e(xi — Xj). (See |73| for the introduction to the theory of Calogero- 
type models and |]B0|], [PIJ] for the new developements.) 

All this discussion was necessary just to illustrate a very simple idea: the 
theory of free particles, though trivial, is in fact unexhaustively deep. It is 
enough to impose sophisticated initial conditions or perform a sophisticated 
change of variables in order to obtain very complicated dynamical systems, 
which after they are studied per se appear to be surprisingly system, the 
reason for this simplicity being that the real underlying dynamics is just 
trivial - that of the free particles,- though it may be very hard probelm to 
reveal this simplicity when the system is given. It is advantage of the general 
theory, that one can begin from the proper side: just from the theory of free 
particles and just start making it more and more complicated, by introducing 
different kind of variables, considering correlators of sophisticated operators 
and so on. Everything what can be obtained in this way is by definition 
trivially integrable, though it may be not so simple to guess for somebody 
who did not know where the particular system at the end of this procedure 
appeared from. 

We now proceed to discussion of particularly important realization of this 

24 For various views and approaches to this theory see Q,[^2), @,fj, @ . , Q , 
fli||,||79|]. (So far there are no connections with Andersson localization in the solid state 
physics.) 
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idea: the theory of <9-operators in \q dimension (i.e. the theory of free holo- 
morphic fields in 2 real dimensions). When considered as functions of moduli 
of the bundles over Riemann surfaces (i.e. of boundary conditions, imposed 
on 2d free fields), these simple objects (known as "r=functions") start look- 
ing a little involved and after all appear related to sophisticated non-linear 
equations (but of course integrable) in 2 and 3 dimensions (like KdV or 
Kadomtsev-Petvishvili (KP) equation). We do not attempt to present an ex- 
haustive theory of r-functions and integrable hierarchies (besides being still 
uncompleted, this is a very big field), but instead concentrate on the very core 
of it, which is just simple determinant formulas for the simplest r-functions 
(namely, associated with free-fermion theory and level k = 1 Kac-Moody 
algebras). This issue will be discussed in some details, because besides being 
the basis of integrable hierarchies theory, it is also exactly the place, where 
the links with the matrix models are found. 

4.2 The notion of r-function 

There are several different definitions of r-functions, but all of them are 
particular realizations of the following idea: r-function is a generating func- 
tional of all the correlation functions in the theory of free particles in 1+1 
dimensions. This basic quantity is a kind of "det D" , where " D" is a time- 
evolution operator (continuous or discrete) and "det" is a sort of a product 
over eigenvalues of " D" , which is usually expressed in the form of a functional 
integral, associated with free particles (it is not a priori Gaussian in original 
variables). This quantity is the most general definition of r-function. 

In practice one is usually more specific. The mostly well studied version 
of r-function arises if one thinks about free particles of a peculiar type: 
free fermions with quadratic Hamiltonian and continuous time evolution, 
i.e. one considers a theory of spin-1/2 b, c-system (fermions) ip(z, z),ip(z, z), 
described by the functional integral 

t{A} ~ Det(d + A) ~ 

~ J D4>Dip exp ^ ^dipj exp 

(LL 

A(z, z)5(z — z)ip(z)ip(z) ) 

where z plays the role of time and A — A(z, z)5(z — z)dzdz is some (|, 1; |, 1)- 
bidifferential (i.e. contains a factor of dz 1 l 2 dzdz 1 l 2 dz). 
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Of course, one can think about more general r-functions, involving many 
fermions (this is often done), and more general b, c- and (3, 7-systems, in 
particular, arising in the context of WZNW model, associated to any Kac- 
Moody algebra of any level. 25 Also of interset is consideration of discrete 
time evolution (described by difference rather than differential equations), 
though, as usual in the 2d theories, this is not a really independent problem. 

In the language of matrix models the restriction to free — fermion r- 
functions is essentially equivalent to restriction to eigenvalue models. Serious 
consideration of non-eigenvalue models, aimed at revealing their integrable 
(solvable) structure will certainly involve the theory of generic r-functions, 
but both these things are matters of the future research, and we" 11 not go 
into details about them in these notes. 



4.3 r-function, associated with the free fermions 

Because of specific form of the Lagrangian in (47 ) the functional integral can 
be easily represented in Hamiltonian form, provided topology of the 2-surface 
on which z, z are coordinates, is trivial (genus 0: sphere or annulus). Namely, 
consider ip and ip as operator- valued functions of z only (not of the time z). 
Then the only reminiscent of kinetic term J d2z tpdtp is canonical commutation 
relation: 



$(z),if>(z)]+ = 5{z - z)dz 1/2 dz 1/2 . (4.7) 

Then 

t{A} ~ (0 I exp (J> £ A{z,z)^{z)i)(z)\ I 0). (4.8) 
Now it is usual to expand around z = 0: 

^)=EneZ^nZ n dz^ 2 ; $( Z ) = £ ngZ ^Z""" 1 dz 1 ' 2 ; 

25 The main technical difference between generic and " free- fermion" cases is that La- 
grangian of generic free field theory is not just quadratic in scalar fields <j>, but can also 
contain particular combinations of exponents eF . It also deserves noting that the most 
general expresiion, quadratic in scalar fields, if rewritten in terms of fermions is in fact 
quartic (but of course not a generic quartic interaction arises in this way). Integrable 
nature of certain quartic-fermion interactions is well known from the theory of Thirring 
models (in this class of models interactions are usually local). 
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ipm | 0) = for m < 0; ^ m | 0) = for m > 0; 
v4(z, z) = J2m,n£Z z m 1 z n A mn dz 1 / 2 dz 1 / 2 ; 



so that 



A(z, z)lp{z)lp(z) = A mni>rJPn- 
dz Jdz m,neZ 

In fact this expansion could be around any point zq and on a 2-surface of any 
topology: topological effects can be easily included as specific shifts of the 
functional A(z,z) - by combinations of the " hadle-gluing operators". Anal- 
ogous shifts can imitate the change of basic functions z n for z n+a and more 
complicated expressions (holomorphic 1/2-differentials with various bound- 
ary conditions on surfaces of various topologies). 

One can now wonder, whether local functionals A(z,z) = U(z)5(z — 
z)dz l l 2 dz 1 / 2 play any special role. The corresponding contribution to the 
Hamiltonian looks like0 



where 



Hcartan = f U { Z )ljj{z)ljj{z) = i U(z)J(z), (4.9) 
Jdz Jdz 



J(z) = ^{z)^{z) = ]T J n z~ n - X dz (4.10) 
nez 



is the U(l)k=i Kac-Moody current; 



J n= 

m+ni \Jmi Jn\ ^^m+n,0- 

(4.11) 

m£Z 



If scalar function (potential) U(z) is expanded as U(z) = Y^k&z% z ■> then 

Hcartan = ^^tkJk- (4-12) 

neZ 



26 Note that normalization factor here is different by a factor of from that in discus- 
sion of discrete models in sections 2.3, 2.7 and 2.8. This is not just a change of notations, 
since Miwa transformation can lead to different results when this normalization is changed. 
See a footnote in section 4.6 below for more detailed discussion. 
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This contribution to the whole Hamiltonian can be considered distinguished 
for the following reason. Let us return to original expression (|4.8|) and try 
to consider it as a generating functional for all the correlation functions of ip 
and if). Naively, variation w.r.to A(z, z) should produce bilinear combination 
if)(z)if){z) and this would solve the problem. However, things are not just 
so trivial, because operators involved do not commute (and in particular, 
the exponential operator in (|4.8|) should still be defined less symbolically, 
see next subsection). Things would be much simpler, if we can consider 
commuting set of operators: this is where abelian U(l) k=1 subgroup of the 
entire GL(oo) fc=1 (and even its purely commuting Borel subalgebra) enters 
the game. Remarkably, it is sufficient to deal with this abelian subgroup in 
order to reproduce all the correlation functions.^] The crucial point is the 
identity for free fermions (generalizable to any b, c-systems): 

:^AMA): = :exp(^ A j) : (4.13) 

which is widely known in the form of bosonization formulas:^] if J(z) = 
dcf>(z), 

4>(A) ~ : e+®> : ( : V(oo)$(A) : = : e^'^ : ) ; 

^(A) ~ : e- m : ( : ^(A)^(oo) : = : e Woo) ^ (A)) : ) . 

This identity implies that one can generate any bilinear combinations of if>- 
operators by variation of potential U (z) only, moreover this variation should 

27 We once again emphasize that this trick is specific for the free fermions and for the 
level k = 1 Kac-Moody algebras, which can be expressed entirely in terms of free fields, 
associated with Cartan generators (modulo some unpleasant details, related to "cocycle 
factors" in the Frenkel-Kac representations , which are in fact reminiscents of free fields 
associated with the non-Cartan generators (parafermions) [ p8| , - but can, however, be put 
under the carpet or/and taken into account "by hands" as "unpleasant but non-essential(?) 
sophistications) . 

28 Formulas in brackets are indeed correct, before them the usual symbolic relations are 
written. Using these formulas we get: 

: ^(A)^(A) : = : e *^~*W . = . e £> . = . X J : 

This identity can be of course obtained within fermionic theory, one should only take into 
account that ^-operators are nilpotent, so that exponent of a single ^-operator would be 
just a sum of two terms (polinomial) . 
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be of specific form: 



\kez ) ' ■ 



]^ z k Jk {z k ~z k 



i.e. 



^=U?-f1 (4 - l4) 



Note that this is not an infinitesimal variation and that it has exactly the 
form, consistent with Miwa parametrization used in the previous section 3. 

Since any bilinear combination can be generated in this way from U(z), 
it is clear that the entire Hamiltonian A mn il) m ip n can be also considered 
as resulting from some transformation of V (i.e. of "time- variables" In 
other words, 

r{A} = A [t]r{A = U}. 
These operators O a are naturally interpreted as elements of the group GL(oo), 



acting on the "Universal Grassmannian" P2| , , ||84| , parametrized by the 
matrices A mn modulo changes of coordinates z — > f(z)- This representa- 
tion for t{A} is, however, not very convenient, and usually one considers 
infinitesimal version of the transformation, which just shifts A 

r{t\A + 6 A} = d SA [t]r{t | A}, (4.15) 

note that this transformation clearly distinguishes between the dependencies 
of r on t and on all other components of A. The possibility of such rep- 
resentation with the privileged role of Cartan generators is the origin of all 
simplifications, arising in the case of free-fermion r-functions.0 

29 It is also the reason, why these are free-fermion tau-functions, that appear in the 
study of ordinary integrable hierarchies: the Hamiltonian flows, which describe evolution 
in different i-directions just commute, because i's are associated with the commuting 
Cartan generators of GL(oo). In the more general situation the flows would form closed 
but non — abelian algebra. 
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Relation ( |4.15| ) is the basis of the orbit interpretation of r-functions [^3| . 



It is also important to understand the role of the " string equation" and othert 
constraints, imposed on r-functions in the theory of matrix models. These 
arise as some particular subalgebras in the set of (9-operators, and their role is 
to specify particular points A in the Grassmannian, of which this subalgebra 
is a stabilizer.^] The simplest examples are in fact provided by formulas from 
the section 2.3 above, where combinations of the screening charges describe 
A's, which are stable points of discrete Virasoro- and W^-constraints (in the 
latter case the multi-fermioia system is used). 

The fact that the r-function at all the points A of Grassmannian can be 
obtained by the group action from t{0}, has an implication, known as Hirota 
equation. The idea is just that there are Casimir operators in the group, 
which commute with the group action and thus the eigenvalue of the Casimir 
operator is the same for t{A} at all A. In the free-fermion case the simplest 
example of Casimir operator is given by 

J Q = jj = j>W= (4.16) 

The eigenvalue of this operator for the vacuum state | 0) is an infinite subtrac- 
tion constant, and this makes equation JqOa | 0) = OaJq | 0) = const ■ Oa | 
0), or Jqt{A} = const • t{A} not very interesting. However, this operator is 
represented in bilinear form and in such cases the following trick is usually 
usefull. 



30 This relation is straightforward in the case of Virasoro constraints, since Virasoro 
algebra is just a subalgebra of the GL(oo) acting on r-functions, and thus is a symmetry 
(covariance) of associated integrable hierarchies |8^] . VF-constraints do not form Lie- 
subalgebra of this GL(oo), they arise after certain reduction, which in turn exists in a 
simple form not everywhere on the Grassmannian (in particular W is not a symmetry of 
entire KP hierarchy |35| ) : here we deal with a more sophisticated self-consistency relation, 
which remains to be understood in full details (e.g., it is unknown, whether reduction 
exists at least at any Virasoro-stable point, which would significantly simplify this kind 
of consideration). In fact, the entire relation between the constraints and r-functions is 
not exhaustively worked out: for example, there is still no clear and satisfactory proof, 
that the full set of Virasoro and/or VK-constraints implies that partition function is a 
r-function, which would be pure algebraic and not refer to the uniqueness of solutions to 



the constraints. Result, widely discussed in the literature, see 29 , is that string equation 
(the lowest Virasoro constraint L-\Z — 0), if imposed on Z, which is somehow known to 
be the properly reduced r-function, implies the entire set of Virasoro and W^-constraints 
(though even this proof can still have some loopholes). 
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If operator, which is bilinear in generators of the algebra, T a T a , commutes 
with the action of the group, so does T a <g> T a , if the group action on tensor 
product of representations is defined as | )<g> | ) — > Oa \ ) <8> Oa \ )■ (Indeed, 
then (T a (g) / + / <g) T a ) 2 commutes with the group action and so does T a <g) 

T a = I (( T a (g, J + J g, T a)2 _ yaya (g, / _ J (g, T a T a ).) If further, T a <g> T a 

annihilates the product of two vacuum states: 

(T a <g>T a ) | 0)® | 0) = 0, (4.17) 
then the same equation holds for all A: 

{T a ®T a ) | | A ) = 0. (4.18) 

Condition ( |4.17 ) is trivially valid in our case: 

E^I°)®^|°) = ' ( 4 - 19 ) 

nez 

since in every term in the sum one of the vacuum states is annihilated: the 
first one of n > and the second one if n < 0.0 Thus we obtain the relation 

E I O a ) ® ^ | O a ) = 0, (4.20) 

which can be now multiplied from the left by 

(0 | i/j(oo)e Hcartan{t) <g> (0 | ■ijj(oo)e Hcartan( - t ' ) 

(t' k do not need to coincide with tk) and after insertions of ^-operators are 
expressed as the shifts of times, we obtain: 

££>-r{t | A}®D+r{t\A} = 0, (4.21) 

neZ 



31 It is easy to verify directly that ^ n ip n ^ip n is indeed a Casimir operator in the tensor 
product: 

En i>n ® V'nj ^ ® Si,™ ^imV'iV'm + Ei,m A lm,1pli'm 8) i] = 
= En (V'n ® E m A nmi>m ~ Ez -^JnV'J ® ^n) = Ei Em ^/toW 7 / ® ^m ~ V/ ® ^m) = 0- 
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where J2n>o D^z n = exp (± J2k>o k^~&~) ■ This is particular form of Hirota 



equation Pq| , which is often used to define r-functions, associated with inte- 
grable hierarchies. If one takes (|4.8|) for the definition, as it is more natural to 
do in the general "theory of everything" and as we did above, eq. ( (4.21| ) is the 
starting point for the path, leading to hierarchies in conventional form of dif- 
ferential equations, the Lax and pseudodifferential representations naturally 
appearing on the way. We do not go along this path in these notes. 

The last remark to be made, before we proceed to more detailed formu- 
las, is that r-functions can be considered as determinants Detd of the d- 
operators acting on fields with some complicated boundary conditions (like 
ip{z) ~ exp (j2k>o tkZ' k ^j in the simplest case of t-dependencies) . Entire A- 
dependence is usually described in this context as that on the point of the 
"universal module space" , which once appeared in the study of string models 
on Riemann surface of arbitrary genus P7| . From this point of view more 
general r-functions are sections of the bundles over universal module space, 
associated with conformal models, more sophisticated then just the theory of 
free fermions (and b, c-systems). The WZNW model is, of course, the most 
important example to be studied in this context. 

The crucial feature of all the quantites, associated in this way to confor- 
mal models is applicability of Wick theorem, reducing multipoint correlation 
functions to pair correlators. In the free fermion case this is just a conse- 
quence of quadratic form of Lagrangian, in generic situation this follows from 
existence of holomorphic operator algebra, which allows to define the corre- 
lators by fixation of monodromy properties, dictated by pairwise collision 
of points. Wick theorem is the concrete source of determinant formulas for 
taw-functions, which are used in order to establish their relations with matrix 
models and other branches of string theory. 

After this discussion of the context where free-fermion r-functions can and 
do appear, we turn now to more detailed and exact formulas, relevant in this 
particular free-fermion case. They are mostly due to Japanese school [155(1 , 
though many other people contributed to this field after it was established. 
The only sophisticated part of the work with these formulas is accurate ac- 
counting for the normal ordering routine which will be mostly unnecessary for 



our purposes. We shall mostly follow the presentation of papers [ 3j| , j30j , [p€ 
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4.4 Basic determinant formula for the free-fermion cor- 
relator 

Let us consider the following matrix element: 

T N {t, i\G} = (N \e H G e R | N) (4.22) 



where 



G = exp I A mn1pm4>n 
\m,ndZ 

fc>0 fc>0 

J(z) = 1p(z)ljj(z) = J n Z~ n ~ l dz] J n = J^^k+n] 

n<=Z k 



(4.23) 



ipm \ N) = 0, m < N; (N | ^ m = 0, m > N; 
4> m | iV) = 0, m > N; (N \ ijj m = 0, m < TV; 
J m | /V) = 0, m > 0; (TV | J m = 0, m < 0. 

The " iV-th vacuum" | N) is defined as the Dirac sea, filled up to the level N: 

oo N-l 

i n) = n i oo) = n ^ i -°°> ^ 

j=7V i=— oo 

oo JV-1 

(TV |= (oo | J] A = (-oo | II ^ ( 4 - 24 ) 

j=7V i=— oo 

where the "empty" (bare) and "completely filled" vacua are defined so that: 

ijj m | -oo) = 0, (-oo \ip m = 0, 

i) m I oo) = 0, (oo | i) m = (4.25) 

for any m G Z. For the only reason that operators J, H, H and G are 
defined so that they have usually ip at the very right and tp at the very left, 
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we have also: 

J m | -oo) = 0, (-00 I J m — 0, 

G ±l I -00) =| -00); (-00 I G ±l = (-00 |; 

e ±H I -00) =| -00); (-00 I e ±H = (-00 I . (4.26) 

Now we can use all these formulas to rewrite our original correlator ( |4.22i ) 
as: 

(N I e H G e S | N) = 




= (-00 I _Q \ e G e \ J] A I -00) = 

\ i=— 00 / V i=— 00 / 

(-00 I e ~* f jf iPi) e H Ge R ( ft ) e~ s | -00) 



(4.27) 



N-1 N-1 



(-00 1 n nt\ n i -°°> 



= Det_ 00<iii<w (-oo I ^[t]^ G [t] I -00) = 
= Deki,j<oT~Li+N,j+N ■ 
The last two steps here were introduction of "GL(oo)-rotated" fermions, 

%[t] = e~ H ^e H - %[t] = e B ^e- R ; Vf[t\ = G^G" 1 , (4.28) 

and application of the Wick theorem to express multifermion correlation 
function through pair correlators 

n ..(t,t) = (_oo I ^[t]^ff[t\ I -00) = 

= (-00 I %[t] G I -00), 

(once again the fact that G^ 1 | —00) =| —00) was used). The only non-trivial 
dynamical information entered through applicability of the Wick theorem, 
and for that it was crucial that all the operators e H , e H , G are quadratic 
exponents, i.e. can only modify the shape of the propagator, but do not 
destroy the quadratic form of the action (fields remain free). This is ex- 
actly equivalent to the statement that "heisenberg" operators are just 
"rotations" of ip, i.e. that transformations ( |4.28| ) are linear. 
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We shall now describe these transformations in a little more explicit form. 
Namely, their entire time-dependence can be encoded in terms of " Shur poli- 
nomials" P n (t). These are defined to have a very simple generating function 
(which we already encountered many times in the theory of matrix models): 

J2Pn(t)z n = exp(jrt k z k ) (4.30) 

n>0 \fc=l / 

(i.e. P = 1, P\= ti, P 2 = \ + t 2 etc.), and satisfy the relation 



dt 
Since 



9P n 

- Pn-k. (4.31) 



oo 



1 



exp [J2t k z k =n J2—t n k k z knk 

\k=l J k>0 \n k >0 Uk - 

Shur polinomials can be also represented as 



Pn(t)= £ IlAC 1 ■ (4-32) 

{n fe} \k>0 ) 

Now, since 

e- B Ae B = A + [A, B] + 1[LA, B],B} + B],B],B] + . . . 

and 

[V'ij <4] = V'i+fc? [[V'i) <4i]; ^fc 2 ] = i ) i+k 1 +k 2 - > ■ ■ ■ i 

we have for every fixed k: 

n k >0 nk - 
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It remains to note that all the harmonics of J in H = J2k>otkJk commute 
with each other, to obtain: 



%(t) = e~ H ^e H = n e ~ tkJk mm n etkJk 

\k>0 J \k>0 

I \ 



n>0 



(4.32) 



n>0 l>i 

Similarly, relation [Jk,ipj] = ipk+j implies, that 

m>j 



n>0 



and finally^ 



T~tij = E i>i (-00 | G Vm | -OO) Pi-i(t)P m ^j(t) 
m>j 

= E Tl m Pl-i(t)P m -j(t), 

m>j 



which implies also that 



dt k 
dU K 



H 



i,j+k- 



Matrix 



- lm 



(-00 G ip m \ -oo) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



32 Eqs.(4.34) can be also interpreted as representations of Shur polinomials in terms of 
fermionic correlators in the bare vacuum: 



P m (t) = (-00 | 4>j +m e H ipj | -oo) ; 
P m (t) = (-oo | ipie H ip l+m | -oo) 



S3 



is the one which defines fermion rotations under the action of GL(oo) -group 
element G: 



Gijj m G 1 = Y,iezipiTim, 

G~ l 4)iG = J2mezTlm^m, OT G^iG' 1 = Y, m zz{T~ 1 )lm'4>m- (4.38) 

If G — 1, Ti m = <W- If all tk = tk = 0, Ttij = Tij. 

4.5 Toda-lattice r- function and linear reductions of Toda 
hierarchy 

In the previous subsection we derived a formula 

r N {t,i | G} = Det itj<0 H i+ N,j+N (4.39) 

for the basic correlator, which defines " Toda-lattice r-function" . For obvious 
reasons i are often refered to as negative-times, r-function can be normalized 
by dividing over the same quantity for all time-variables vanishing, but this 
is not always convenient. Eq. (|4.39| ) has generalizations - when similar matrix 



elements in a multifermion system is considered - this leads to "multicompo- 
nent Toda" (or AKNS) r-functions. Generalizations to arbitrary conformal 
models should be considered as well. It has also particular "reductions", of 
which the most important are: KP (Kadomtsev-Petviashvili), forced (semi- 
infinite) and Toda-chain r-functions. This is the subject to be discussed in 
this subsection. 

Idea of linear reduction is that the form of operator G, or, what is the 
same, of the matrix TJ m in eq.( |4.35| ), can be adjusted in such a way, that 
r N {t,t | G} becomes independent of some variables, i.e. equation(s) 

f E a W- + E a-sr + £ PkD N (k) + 7 ) r N {t, t | G} = (4.40) 

\ k k k k k J 

can be solved as equations for G for all the values of t, i and iV at once. 
(In ( |4.40j ) D N {k)f N = /jv+fe — In-) m this case the system of integrable 



equations (hierarchy), arising from Hirota equation for r, gets reduced and 
one usually speaks about "reduced hierarchy". Usually equation ( |4.4U| ) is 
imposed directly on matrix H^, of course than ( |4.40|) is just a corollary. 
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We shall refer to the situation when ( 4.40 ) is fulfilled for any t, i, N as 
to "strong reduction". It is often reasonable to consider also "weak reduc- 
tions", when Q4.40 ) is satisfied on particular infinite-dimensional hyperplanes 
in the space of time-variables. Weak reduction is usually a property of entire 
r-function as well, but not expressible in the from of a local linear equa- 
tion, satisfied identicall for all values of t, i, N. Now we proceed to concrete 
examples: 

Toda-chain hierarchy. This is a strong reduction. The corresponding 
constraint (f4.40|) is just 



dHij 



dTiij 



or, because of (|4.36|) , TC 



i+k,j 



dt k dt k 

7iij + k- It has an obvious solution: 



n 



n 



i+3 ' 



(4.41) 



(4.42) 



i.e. Tiij is expressed in terms of a one-index quantity 7ij. It is, however, not 
enough to say, what are restrictions on Hij - they should be fulfiled for all 
t and i at once, i.e. should be resolvable as equations for T[ m . In the case 
under consideration this is simple: T/ m should be such that 



T, 



III! 



■ l+m- 



(4.43) 



Indeed, then 

n 

and 



'j 



J2l,mTlmPl-i{t)P m -j(t) — J2l,mTl+mPl-i{t) P m -j(t) 



Y,n>oTn+i+j (Dfc=0 Pk(t)Pn-k(t)) 



n l = Y J Tn +l [Y. p k{t)Pn~k{t))- 

n>0 \fc=0 / 



(4.44) 



Volterra hierarchy. Toda-chain r-function can be further weakly reduced 
to satisfy the identity 



dr. 



2N 



dt 



2k+l 



0, for all k, 



(4.45) 



21+1 



=0} 
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i.e. t 2 n is requested to be even function of all odd-times £22+1 (this is an 
example of "global characterization" of the weak reduction). Note that Q4.45Q 
is imposed only on Toda-chain r-function with even values of zero-time. Then 
( |4.45|) will hold whenever Hi in ( |4.44|) are even (odd) functions of t odd for 
even (odd) values of i. Since Shur polinomials Pkif) are even (odd) functions 
of odd-times for even (odd) k, it is enough that the sum in ( |4.44| ) goes over 
even (odd) n when i is even (odd). In other words, the restriction on TJ m is 
that 

T lm = f l+m , and f 2k +i = for all k. (4.46) 

Forced hierarchies. This is another important example of strong re- 
duction. It also provides an example of singular r-functions, arising when 
G = exp (j2 Amn'ipm'ijjn^ blows up and normal ordered operators should be 
used to define regularized r-functions. Forced hierarchy appears when G can 
be represented in the form |SI| G = G P+, where projection operator P + is 
such that 

P + I N) =| N) for N > N , 

P + I N) = for iV < N Q . (4.47) 

Explicit expression for this operator is 33 

P+ = ■ exp [ - £ M] : = J] (1 - M) = II Ml- 

V KNo J KNo KNo 

Because of ( |4.47| ), P + \ —00) =), and the identity G \ —00) =| —00), which 
was essentially used in the derivation in (4.27), can be satisfied only if Go is 
singular and T} m = 00. In order to avoid this problem one usually introduces 
in the vicinity of such singular points in the universal module space a sort of 
normalized (forced) r-function rjy = One can check that now T/ m = 00 

for all l,m < No, and t' can be represented as determinant of a final- 



dimensional matrix |90|,[[39|: 



r f N = Det2v <ij<jvWj- for > N ; 



33 Normal ordering sign : : means that all operators ifj stand to the left of all operators 



tp. The product at the r.h.s. obviously implies both the property (4.47) and projection 
property P? = P + . 
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4 = i; 

4 = for iV<iVo. (4.48) 

For N > Nq we have now determinant of a /imte- dimensional (N — No) x 
(N — Nq) matrix. The choice of No is not really essential, therefore it is better 
to put iV = in order to simplify formulas, phraising and relation with the 
discrete matrix models (iVo is easily restored if everywhere N is substituted 
by iV — No). For forced hierarchies one can also represent f as 



r 



ft = Det <i, j<N didin f , (4.49) 



where 1-0 = Hq and d\ — di — For forced Toda- chain hierarchy 
this turns into even simpler expression: 



T 



f 



N = Det < i>j<N d{ +j H f , (4.50) 



while for the forced Volterra case we get a product of two Toda-chain r- 
functions with twice as small value of N pTf : 



t{ n = (pet < iJ<N dt j H f ) ■ (Det Q < ij<N d l 2 +1 (d 2 Hf)] 

= 4W\ ■ T^d 2 Hf}. (4.51) 

Forced can be always represented in the form of a scalar-product ma- 
trix model. Indeed, 

Wii = Y. T i™ P i-ti) P ™- = f je u{h)+0{h) WVT{h, h)dhdh, (4.52) 

where T{h,h) = Y^im^imhr^hr^ 1 , and e u ^ = e^>o tkhk = Ei> h l Pi(t). 
Then, since Deto< i j < ]yh l = A N (h) - this is where it is essential that the 
hierarchy is forced - 

Det <ij<^ = 11/ / e u ^ +0 ^T(hi, h^dhdh ■ A N (h)A N (h), (4.53) 

i 

i.e. we obtain a scalar-product model with 

dfi h - h = e u{h)+0Ch) T{h, h)dhdh. (4.54) 
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Inverse is also true: partition function of every scalar-product model is forced 
Toda-lattice r-function - see section 4.7 for more details. 

KP hierarchy. In this case we just ignore the dependence of r-function 
on times t. Every Toda-lattice r-function can be considered also as KP 
r-function: just operator G KP = Ge H (a point of Grassmannian) becomes in- 
dependent. Usually iV-dependence is also eliminated - this can be considered 
as a little more sophisticated change of G. When N is fixed, extra changes of 
field- variables are allowed, including transformation from Ramond to Neveu- 
Schwarz sector etc. Often KP hierarchy is from the very beginning formulated 
in terms of Neveu-Schwarz (antiperiodic) fermionic fields (associated with 
principal representations of Kac-Moody algebras), i.e. expansions in the first 
line of ( |4.23| ) are in semi-integer powers of z: ^ns( z ) = J2nez V^- 2 ™ '^dz 1 ' 2 . 



Given a KP r-function one can usually construct a Toda-lattice one with 
the same G, by introducing in appropriate way dependencies on i and N. 
For this purpose t kp should be represented in the form of ( |4.39|) : 

r KP {t | G} = Det id . M!) 1 '. (4.55) 

where 'Hf^ = J2iTijPi-i(t). Since XJ m is a function of G only, it does not 
change when we built up a Toda-lattice r-function: 

TN\t,t | G} = Detij <0 Hi + N,j+N', 

Hij = Y,i, m TimPi-i{t)Pm-A?) = LC^-i^l ( 4 -56) 



Then 



r KP {t | G} = r {t,0 I G}. (4.57) 



If we go in the opposite direction, when Toda-lattice r-function is considered 
as KP r-function, 

T {t,i\G} = T KP {t\G(t)}; 

n* p = J2n im p m -At) and 

m 

Tij{G{i)} = J2 T i™{G}P m -j(t). (4.58) 

m 

KP reduction in its turn has many further weak reductions (KdV and 
Boussinesq being the simplest examples). We shall mention them again in 
section 4.9 below, after Miwa transformation of representation fl4.39|) will be 
considered in the next subsection. 
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4.6 Fermion correlator in Miwa coordinates 



Let us now return to original correlator ( |4.22| ) and discuss in a little more 
details the implications of bosonization identity (|4.13|) . In order not to write 
down integrals of J, we introduce scalar field:Q 

<K Z ) = E ^z k + <f> + Jologz, (4.59) 
kez-o 



such that d(j)(z) = J(z). Then ( [4.131) states that: 

: V(A)V5(A) : = : e^-^ : (4.60) 

"Normal ordering" here means nothing more but the requirement to neglect 
all mutual contractions (or correlators) of operators in between : : when 
Wick theorem is applied to evaluate corrletion functions. One can also get 
rid of the normal ordering sign at the l.h.s. of ( |4.60D , then 

V>(A)V>(A) = : e^' x) : : e~ m : (4.61) 

In distinguished coordinates on a sphere, when the free field propagator is 
just log(,2 — z), one also has: 

i>(z)4>(z) = : i){z)i)(z) : 

Our task now is to express operators e H and e H through the field <fi. This 
is simple: 

H = jf U(z)J(z) = jf U(z)d(f)(z) = ~j <j>(z)dU(z). (4.62) 

Here as usual U(z) = J2k>o^k zk and integral is around z — 0. This is very 
similar to generic linear functional of (f>-(\) = — J2k>o \Jk^~ k , 

H = J 0_(A)/(A)dA, (4.63) 



3 One can consider (j> as introduced for simplicity of notation, but it should be kept 
in mind that the scalar-field representation is in fact more fundamental for generic t- 
functions, not related to the level k = 1 Kac-Moody algebras (this phenomenon is well 
known in conformal filed theory, see |16| for more details). 
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one should only require thatQ 

/(A) 



i.e. 



dU(z) = I ^rd\, 
J z — A 

U{z) = J log (l - ~) f(X)dX. (4.64) 

In terms of time-variables this means that 

t k = — J X~ k f{X)dX. (4.65) 

Here we required that U(z = 0) = 0, sometimes it can be more natural to 
introduce also 

t = J log Xf(X)d\. (4.66) 

This change from the time- variables to "time density" f(X) is known as Miwa 
transformation. In order to establish relation with fermionic representation 
and also with matrix models we shall need it in "discretized" form: 

h = | (E 7 A~ fc - E 7 A^ fc ) , 

to = -£ (E 7 log A 7 - E 7 log A 7 ) . (4.67) 



We changed integral over A for a discrete sum (i.e. the density function 
/(A) is a combination of (^-functions, picked at some points A 7 , A 7 . This is 
of course just another basis in the space of the linear functionals, but the 
change from one basis to another one is highly non-trivial. The thing is, 
that we selected the basis where amplitudes of different (^-functions are the 
same: parameter £ in ( 4.67 ) is independent of 7. Thus the real parameters 
are just positions of the points A 7 , A 7 , while the amplitude is defined by 
the density of these points in the integration (summation) domain. This 
domain does not need to be a priori specified: it can be real line, any other 
contour or - better - some Riemann surface.) Parameter £ is also unnecessary 

35 As it is usual nowadays, a factor of 2iri is assumed to be included into the definition 
of contour integral § . 
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to introduce, because basises with different £ are essentially equivalent. We 
shall soon put it equal to one, but not before Miwa transformation will be 
discussed in a little more detail. 

Our next steps will be as follows. Substitution of Q4.63|) into fl4.67|) , gives: 



H = -£^0_(A 7 ) +£^0_(A 7 ). (4.68) 

7 7 

In fact, what we need is not operator H itself, but the state which is created 
when e H acts on the vacuum state (N |. Then, since (N \ J m = for m < 0, 
(N | e ~^- (A) is essentially equivalent to (N | e"^ (A) with (j)-(X) substituted 
by entire 0(A). If £ = 1, e - *^ can be further changed for ip(X) and we obtain 
an expression for the correlator (|4.22| ) an expression where e H is substituted 
by a product of operators ^(A 7 ). The same is of course true for e . Then 
Wick therem can be applied and a new type of determinant formulas arises 
like, for example, 

t ~ ^W(A) det ^ (iV 1 ^ X ^ { ~ Xs) G 1 N) (4 ' 69) 



It can be also obtained directly from (|4.27 ), ( 4.2S ) and ( [4.35 ) by Miwa trans- 



formation. The rest of this subsection describes this derivation in somewhat 
more details. 

The first task is to substitute 0_ by 0. For this purpose we introduce 
operator 

oo 

t k J k = H + + H„, (4.70) 

fc=— oo 

where H + = J2k>o^kJk is just our old H, H_ = J2k>ot-kJk, and "negative 
times" t-k are defined by "analytical continuation" of the same formulas 
flOD and flUTTp : 



t-k = U X k f(X)dX = -f (E K ~ E A?) • (4-71) 



7 7 



Then 



£ t k J k = ff + + #_ = -M E 0(A 7 ) - E <K\) • ( 4 -72) 

k=-oo V 7 7 / 



91 



Further, 

e H ++ H - = e -^e H+ e H - = e^e H ~e H+ , (4.73) 

where 

S (t) = £ kt k t„ k = -e e I (e ( v - E (a* - ^ 

k>0 k>0 ft \ 7 8 



eiog 



S 1-^)1-^)' 



where prime means that the terms with 7 = 5 are excluded from the product 
in the numerator and accounted for in the infinite "constant", added at the 
r.h.s. In other words, 

»>M= const Vn 7 >,(^)(V-*^ a 



ilil^-^) 



Since (JV | J m = for all m < 0, we have (iV | e = (N |, and therefore 

(N\e H = (N\ e H + = (N | e H ~e H + = e~^{N | e H++H ~ . (4.76) 
From eg. (|4.72|) , 



e 



= const ■ J] = e-^ ] : : e*^> : (4.77) 



where "const" is exactly the same as in ( [4.751 ). If £ = 1, eq.(4.61) can be 
used to write :0 

(, ' |( ' = ^is) (i ' l ?«) (4-78) 



36 The choice of £ can be dictated by particular purposes. Here we impose the re- 
quirement on Miwa transform to represent e H = e HcaTtan as a product of dimension- 1 /2 
operators - this is most natural from the point of view of Hirota equations and simplifies 
the relation with integrable hierarchies. However, in section 2.7 and 2.8 we used another 
requirement (and there £ = rather than £ = 1). There the 1-matrix model was con- 
sidered, which is characterized by especially simple form of the full hamiltonian (product 
of dimension-zero operators), and it was more important to adjust operators, which arise 
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Similarly, 



e s i n) = n^) i n ) A(A,A) - ■ 

T T A2(A)A 2 (A) 



where 



(4.79) 



(4.80) 



and we used the fact that J m | AT) = for all m > 0. Finaly, 

A(A, A) A(A, A) 



r N {t, t\G} = (N \e H G e n \ N) 



A2(A)A2(A)A2(A)A2(A) {u 
N). 



■ (n i n ^(a 7 ) n ^(a 7 ) g n v>&) n 

7 7 <5 5 

Singularities at the coinciding points are completely eliminated from this 
expression, since poles and zeroes of the correlator are canceled by those 
coming from the Van-der-Monde determinants. 

Let us now put N = and define normalized r-function 



T {t,t\G} 



To{t,t 


G} 


r {0,0 


G} 



(4.82) 



i.e. divide r.h.s. of Q4.81Q by (0 | G | 0). Wick theorem now allows to rewrite 
the correlator at the r.h.s. as a determinant of the block matrix: 



det 



( (0 1 


^(A 7 )^(A 5 ) G 


0) 


(0 


V(A 7 ) G 4>{h) 


|0> 




(0 | G | 0) 






(0 | G | 0) 




(0 


| $(\s) G V(A 7 ) 


|0> 


(0 


G ^(XjMh) 


|o> 


V 


(0 | G | 0) 






(0 | G | 0) 





(4.83) 



/ 



from e Hcartan after Miwa transform to have simple correlators with e A ^ . When ana- 
lyzing the 1-matrix model from this point of view one should also keep in mind that it 
was actually represented in s.2.3 in terms of two complex fermions. Screening charges 
are Q<+) = § = § ■ipxfa = § e^ 1 "* 2 , = § = § = § e^ 2 "* 1 , while 

4> = 75(^1 _ fa)- The Hamiltonian is H Ca rtan = ■^Y.k t kJk = |IZk**( J k _ J k)> and 
Miwa transformation generators insertions of operators X1X2, where \i an d Xi have di- 
mension 1/8 (rather than 1/2 as in the one (complex)-fermion system, considered in this 
section). 
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Special choices of points A 7 , . . . , A<s can lead to simpler formulas. If A 7 — > 
A 7 , so that tk — > 0, the matrix elements at the right lower block in ( 4.83Q 
blow up, so that the off-diagonal blocks can be neglected. Then 



r {t,t | G} — > r KP {t | G} 
A(A,A) 



(0| 


e H G 


10) 


(0 


1 G | 


0) 



^det 



A 2 (A)A 2 (A) 



~/5 



(0 


^(A 7 )V(A ( 


') G 


0) 


(0 


G 


0) 



(4.84) 



This function no longer depends on t-times and is just a KP r-function. 
Matrix element 



^■ A)= (0IGI0) (4 ' 85) 
is singular, when A — ► A: <^(A, A) — > —^r. If now in ( [4.84|) all A — ► oo, 



T K P{t | G} = det^CA^ (4 g6) 



where 



<^(A) = (0 | V(A) («9 5 -^) (oo) G | 0) ~ A 5 " 1 (l + O Q)) . (4.87) 

This is the main determinant representation of KP r-function in Miwa parametriza- 
tion. 

Starting from representation ( |4.86| ) one can restore the corresponding 
matrix Hff in eq. flQgP @: 

r^ p {t} = Iz'ip^iz^k^dz, (4.88) 

i.e. 



= j> z l ^(z). (4.89) 

Then obviously = Ti. i+k j. Now we need to prove that the r-function 

is given at once 

by dCt A(x) Xs) and ^tnfj p {t}. In order to compare these two 
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expressions one should take tk = \ Y% -\ fc > so that 



exp ( tk zk 

,k>0 



7 =1 A 7 ~ Z 



n A 7 e 



R 7 A 7 (A) 
^ - A 7 A(A) ; 



(4.90) 



where 



A 7 (A) = J] (A Q - A^ 



A(A) 



a>/3 



(A a — A~ 



(4.91) 



and 



KP 



(4.92) 



As far as n is kept finite, determinant of the infinite-size matrix ( ^4.92|) , 
^ , — since it is obvious from (^4.92|) that the rank of 



Det i j <0 'Hf j p 



the matrix is equal to n. Therefore let us consider the maximal non-vanishing 
determinant, 



^Gt-n<i,j<o'Hij 
f n \n f ( 

J|A 7 det 



KP 



^ +1 A 7 (A) 



Uj-y 



A 7 A(A) 
det 7 j^(A 7 ) 



det 7 j(/3j(A n 



(4.93) 



A(A) 



We used here the fact that determinant of a matrix is a product of determi- 



95 



nants and reversed the signs of % and j. Also used were some simple relations: 



^ A 7 (A) _ 1 



4 A(A) A»(A)' 



A 7 X7 



1 1 



n 



-(n-1) 



thus 



a (! a) = n t — it = (-) n(n_i)/2 A(A) n ^ 



HA^H^-^n^^det, 1 ' 



A !A(A) 7 A 7 A(A) ' 



Since (|4.93| ) is true for any n, one can claim that in the limit n — > oo 
we recover the statement, that r ifP {t} = Detij <0 Hfj P with 7i^ P given by 



eq.( |4.92| ) (that formula does not refer directly to Miwa parametrization and 



is defined for any t and any j < and i). This relation between <£> 7 's and 
7~Cfj P can now be used to introduce negative times ik according to the rule 
(|4.58|) . Especially simple is the prescription for zero-time: TCij — > Tli + N,j+N, 
when expressed in terms of (p just implies that 

det (p^(X s ) det ip 1+N (X s ) 

A(A) ~* (detA) Ar A(A) ' 1 ' j 

Generalizations of (|4.88|) , like 



Hiji^t} = j> iz^iO | ij){z) G ${z) | 0)eEk^ zk +^dzdz, (4.95) 
also can be considered. 



4.7 Matrix models versus r-functions. 

We are now prepared to return to our main subject and discuss integra- 
bility properties of eigenvalue matrix models. The claim is that partition 
functions of all these models, when considered as functions of time-variables 
(parametrizing the shapes of potentials) are in fact r-functions of (perhaps, 
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multicomponent) Toda-lattice and/or KP type. (Interesting non-eigenvalue 
models are believed to be related to integrable systems of more general type, 
not restricted to level k — 1 Kac- Moody algebras.) 

Partition functions are, however, not generic Toda or KP r-functions: 
first, they usually belong to some reduced hierarchies, second, the relevent 
operators G (points of Grassmannian) are restricted to stay in peculiar do- 
mains of the universal module space, specified by "string equations". String 
equation is in fact nothing but the set of Ward-identities (Virasoro or W- 
constraints in the examples under investigation), which are now interpreted 
as equations on G. The very possibility of such interpretation is highly non- 
trivial and reflects some deep relation between the constraints and integrable 
structure. In the case of Virasoro constraints this is not a puzzle, because 
Virasoro algebra is a symmetry (covariance) of the hierarchy, the situation 
with other constraints is less clear (see the footnote at section 4.3). In fact, 
when applied to a r-function of appropriately reduced hierarchy, the infinitely 
many constraints usually become dependent and it is enough to impose only 
the lowest Virasoro constraint L_it = (or £- p t = 0, where p is the degree 
of reduction), in order to recover the entire set f29| . It is this lowest con- 
straint (or rather its ti-derivative, ^-(L_ir) = 0) that traditionally carries 
the name of " string equation" . It is often much simpler to deduce than the 
entire set of identites, what is important in practical applications (especially 
because determinant formulas, which imply integrability, can be also simpler 
to find in some situations than the Ward identities). 

In order to give a complete description of some sort of (matrix) models 
from the point of view of integrability theory it is enough to specify the hier- 
archy, to which it belongs (if partition function is interpreted as r-function, 

ZmodclW = r{t | Gmodcl} (4.96) 

and the string equation which serves to fix operator G - the point in the 
universal module space. Q After that it becomes an internal (yet unsolved) 
problem of integrability theory to explain, what is so special about the set of 

37 As we argued in the Introduction and in section 2.1, the word "matrix" can be 
probably omited if generic Lagrangians are considered in other models of quantum field 
theory. Also the universal module space (where moduli are - of bundles over spectral 
Riemann surfaces) can (and should) be treated as a " space of theories" . It is one of the 
great puzzles (and beauties) of string theory, that Riemann surfaces appear both in the 
world-sheet and spectral " dimensions" . See § for more discussion on this issue. 
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points {Gmodei} in this space. (We shall touch this problem in the next sub- 
section, devoted to Kac-Schwarz operators.) Alternatively, if there is nothing 
special, it is an (unsolved) problem of matrix model theory to find models, 
associated with any points G in the universal module space (or explain what 
is an obstacle, if any). 

We proceed now to description of particular matrix models from this 
point of view. As everywhere in these notes we consider only the most impor- 
tant classes of scalar-product, conformal (multicomponent) and Generalized 
Kontsevich models (GKM). All other examples (like models of complex, or- 
thogonal, unitary etc matrices) can be included into consideration with more 
or less effort (see |28| and |H] for the cases of complex and unitary models 
respectively), but they do not add much for the general theory that we are 
now considering. String equations will be discussed in the next subsection. 

Scalar-product models. These were exhaustively discussed in sections 3.5- 
3.7. We remind that all conventional multimatrix models (with inter-matrix 
interaction of the form exp(Trif ( a ) if( Q+1 )) belong to this class. The crucial 
formulas are: 

Zn = DetTvHjj = Deto<i,j<jv-i7Xy = 
= Det-N<i,j<o7~ii + N,j+Ni 

/L y dtidtj 11 \dtij \dtj ,L ■ 

Here 

K = I ~ hj ) = J dfae^WW 

Further, 

e U(h) = e J2 k > ^ k = Y i h l P l (t); 

i 

m 

and thus 

TL = ((h l I /> m », (4.100) 



(4.97) 
(4.98) 

(4.99) 
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where the scalar product (( | )) is w.r.to the measure dfi h ^ (while ( | ) is 
w.r.to dn Kh = fW+vMdfa). 

One would immediately recognize in these formulas representation ( |4.39| ) 
of Toda-lattice r-function, be there no additional restriction that determinant 
in ( 4.97|) is over finite-dimensional N x N matrix (indices are constrained: 



hj > —N). This can be automatically taken into account if we require that 

T/ m = oo for all I, m < 0, (4.101) 

and identify Zn as a r-function t* of forced Toda-lattice hierarchy (thus the 
supescript / carried by Ti and T). We conclude that partition functions of 
any scalar-product model is a r^-function of forced Toda-lattice hierarchy. 

Let us now consider them as KP r-functions. This means that the 
t-dependence is simply ignored. However, N will be preserved explicitly 
as a parameter, labeling KP r-function. After Miwa transformation = 
— \ J2-y ~ r k, described in section 3.7, we get: 

7 det y5 QN+-y-i(Xs) in9 ^ 
N - ~K{X) ' (4.11)^) 

where Q are orthogonal polinomials w.r.to the measure dv h h = e~ rfc/l * dfi h ^. 
We conclude that in the framework of KP hierarchy the scalar-product mod- 
els are distinguished by the fact, that the corresponding y? 7 (A) in ( f4.86p are 
polinomials rather than infinite series in powers of A -1 . 

1- Matrix model. This is particular example of scalar product model with 
a local measure dfi h ^ = e u ^ +u ^5(h — h)dhdh. In this case 

Hlj = (h* I V) = (h i+ *) = ^-H f = +3 H f . (4.103) 

Thus in this case we deal with the (forced) Toda-chain reduction of Toda- 
lattice hierarchy. In the end of this section 4 we use orthogonal polinomials to 
present a detailed description of 1-matrix models as a Toda-chain r-function. 

This model can be alternatively defined as Gaussian Kontsevich model: 
see section 3.8. The fact that partition function is a r-function follows then 
from the general statement for GKM, see below. The fact that it is forced 
r-function is related to the property C-n = 0, mentioned at the end of s.3.8 
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(and proved in s.3.9). Also reduction to Toda-chain hierarchy can be observed 
directly in terms of GKM: see ref.|36| for more details. 

Multicomponent (conformal) matrix models. These are related to multi- 
component hierarchies, with r-functions representable as correlators in multi- 
fermion systems. An example of determinant formula which substitutes 
(|4.39| ) in the 2-component case, is given at the end of section 3.5, where it is 



derived from consideration of the relevant matrix model |p9|| . For derivation 
of the same determinant formula in the theory of r-functions see ref . [|9^ 



Generic theory of multicomponent hierarchies is now making its first steps 



and we do not review it in these notes. See |93[] for the group-theory approach 
to the problem. 

Generalized Kontsevich model. Determinant formulas for this case are 
derived in Section 3.3. The most important expression is 

Zy{iV ' T} " (de^F W) ' ( } 

where 

Cf^X) = ^. e -A7'(A)+V(A)^) J x J-l e -V( X ) + V'(X) Xdx = 

= A^(l + 0(\- 1 )); (4.105) 

and 

y> 7 (A) = ^ 7 -i(A) = A 1 " 1 ®^). (4.106) 

For N = this is just the representation, peculiar for KP r-function in Miwa 
parametrization = j-trA~ k , see eq. (|4.86j ) above. Thus 



Z v {T} = t kp {T\G v }, (4.107) 

where it is operator G (the point in Grassmannian) which depends on the 
shape of potential V(X). We also remind that the only way in which Z 
depends on the size of the matrix n is through the domain of variation of 
the time variables T. If ( [1.104 ) is extended to full Toda-lattice r-function, 



by introduction of negative times, we get 



x f dX(detX) N exp ( -tiV(X) + trAX + J2 T ^ trX 



(4.108) 



— k 

L k UXW 

fe>0 



100 



When this extended partition function is considered as KP r-function, we 
have instead of ( M04| ): 

and relevant ^-functions are 

= A^+7-i(i + o(A- 1 )), (4.110) 

with 

V{x) =V{x)-N\ogx-T k>Q f k x- k , 

V(x) = V+(x) (4.111) 

(where V+(x) is the positive-power fragment of Laurent series V(x)). Func- 
tions <£> 7 (A) in (|4.105|) are just equal to (py^(X) 



T=0 

4.8 String equations and general concept of reduction 

The role of string equation is to fix the point G in the universal module 
space (UMS), associated with the particular matrix model, so that partition 
function, considered as a function of time variables, will appear as the cor- 
responding r-function of a fixed shape. In this sense the idea behind the 
string equation is exatly the same as reduction of integrable hierarchies. The 
difference is that linear reductions, as defined in the section 4.5 above, are 
not enough to fix G unambiguously: they just specify certain subsets in the 
Grassmannian, which are still infinite-dimensional. The reason, why these 
are usually linear reductions that are considered in the conventional theory of 
integrable hierarchies, is that they are associated with the simplest possible - 
Kac- Moody - subalgebras in the entire GL{p6). String equations, even their 
simplest examples, are usually fragments of more complicated - Virasoro and 
VF-algebras, and are in fact considerably more restrictive. Moreover, string 
equation is usually a distinguished fragment, because it usually belongs to the 
Virasoro component of Ward identities, and Virasoro is still a Lie subalgebra 
in GL(p6). This is what makes the problem of string equations very similar 
to the "classical" one with linear reduction. 
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More specifically, in order to include string equations (and in fact the en- 
tire set of Virasoro - but not W - constraints) into consideration of reduction 
it is enough to allow the coefficients in ( |4.4U| ) to depend on t and i, without 
changing the order of time-derivatives. Of course, there are no obvious rea- 
sons to think that any point G in the UMS can be selected by imposing this 
kind of linear-derivative constraints on r-function, - and further investigation 
can require essential generalization of such restricted notion of string equa- 
tion. However, some of the eigenvalue matrix models are already known to 
possess string equations of such simple type, associated with Virasoro sub- 
algebras of GL(oo). We do not go into details of the general theory - it is 
far not completed yet, - but instead present several examples of how string 
equations arise in particular matrix models. These examples can illustrate 
also the simplifications arising when only string equations and not the entire 
sets of Ward identities need be derived. In particular, it is clear that in cases 
when r is represented as Det^Tiij, & linear differential equation imposed on 
Tiij will give rise to a similar equation on r itself. Most of known string equa- 
tions can be derived with the help of this technical idea. They are usually 
associated with invariance of integrals under constant shifts of integration 
variables Sh = const in scalar-product and other discrete models, and with 
the action of operator tr^- in GKM. For somewhat more involved ideas, 



associated with string equations see [[54[ . 

Scalar-product models. String equation can be easily deduced for very 
specific type of measures dfi h j l . Since integral 



Hi 



Kh?e v ^ +d ^dii h - h 



(4.112) 



is invariant under the shift of integration variable 5h = const, 



dU(h) d , . 



0, (4.113) 



or 



d 

fe>0 Ct fc _i 




(4.114) 



String equation arises straightforwardly, when operator S is linear. This is 
true, if log(dp, h j l ) ~ hf(h) with any function f(h). If the measure dfi h h is 
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also required to be symmetric in h and h, we obtain conventional 2-matrix 
model as the only example: 

dfi h - h = e ch ~ h dhdh. (4.115) 



Equation for TCij is: 




+ C W- I = (4-116) 




Its implication for fjv is: 

« « \ 

= 0, (4.117) 

since the r.h.s. of ( |4.116|) does not contribute to determinant (the entries in 
the i-th row are proportional to those in the % — 1-th row). 

In particular case of 1 -matrix model c = 0, and we recognize the lowest 
Virasoro constraint L_ifjv = 0. Traditionally the name of string equation is 
given not to L_i-constraint itself, but to its ti-derivative: ^-(L_ifjv) = 0. 
For 2-matrix model (|4.117|) is the lowest (m = 1, n — 0) component of 



the Ward identities [W { ™^\t) - (-) m+n c n+1 W^\t)J f N = 0. Of course, 
there is also a similar equation with t «-> i. 

Multicomponent (conformal) models. The crucial feature of these mod- 
els is that intermatrix interaction, when rewritten in terms of eigenvalues, 

usually contains only differences h\ (X " > — . Thus there is usually covari- 

ance under simultaneous shift of all eigenvalues oh\ = const by a same 
constant. This gives rise to a string equation of the form 



J2L {a 1 ) )r^ = 0. (4.118) 
a 



See |39| for details. 



Generalized Kontsevich Model. 

In order to derive string equation one should act on partition function 
Z v {T k = itrA-*} = Cy X T v {L = V'(A)} with operator tr^g- = tr^^. 
We can rewrite the result of this action in terms of time-derivatives, 

tr^logZWT} = -£ («v^)£r^Z r{ T } . (4,19, 
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Alternatively we can use the fact that tr 



d 



E \r it 



1 



d 



, I = v'(\), 



j tr 1 V"(Xy) <9A 7 
and explicit expression for Z v in terms of eigenvalues (Miwa coordinates), 



Zy ~ e 



trV(A)-trAV'(A) 



/n^(A 7 ) 



det <p 7 {\s) 
A(A) 



det (p-y(\& 
A(A) 



(4.120) 



to get: 



( tr ^) l0gZ ^ T > = 

1 V"'{K) I V'iX^-V'iXs) 

2 (\/"(A))2 + 2 7 ^ A 7 -A 5 

d 



V"(\)V"{\ 5 ) 



(4.121) 



-trA + ^2 "H7~ log det 75 (^ 7 (/ 5 ). 
Comparison of these two expressions gives: 



v 



1 



sK"(A) 



r z. r - v v (A)A^JdT k 

1 y»(\J - V"(\ 5 ) l d_ 

2 ±1 A. - A* ' V"(\ 1 )V"(\ S ) m 



7 >5 ^7 ~~ ^ 

9 



(4.122) 



9 



. n log Z y + trA - lo § det 7 5<£ 7 (fc). 

Cii ^ 

One can show that the r.h.s. is equal to zero, and thus the string equation 
arises in the form 



0. 



(4.123) 



If potential is monomial, V p = ^ppj-, then r k = — ^-^5, 



and 



r v p i f = 



1 

P 



2 (A; + p)(T k+p + r k+P )^r + ^ £ Mp - k)T k T p _ k 



k>0 



P-1 

£• 

fe=i 



(4.124) 
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The technical idea behind the proof BO] is to represent 



9 ^^ z v{T k + ^}d\ 



log Z v = Res ^ JT l , (4-125) 

Oli Z v {lkf 

and make use of the second determinant representation in ( (4.120| ) both in 
the denominator and the numerator: 



<^5(A 7 ) (p n+ l(\. 



1) 



det 

d dX V V^( A ) Vn + i(X) 

logZv/ = Res — — (4.126) 



m to n"=i(A-A 7 ) det^(A 7 ) 

Now we recall that 

V5 7 (A) ~ A 7 - 1 (l + G(A- 1 )). (4.127) 
At some moment we shall need even more: in fact 
V5 7 (A) ~ A^ 1 (1 + 0(A~ 2 )), i.e. 

y 7 (A) = A 7-1 + c 7 A 7-2 + . . . , and c 7 = for any 7. (4.128) 

This is a rather delicate property of GKM, it follows from two facts: first, 
/ ynn (y'"Y\ 

if 1 = 1 + O \ — , - — , thus Ci = 0, and second, Kac-Scwarz operator 

\{v y ( v "y j 

A, defined in eq. ((4.106|) below, does not have contributions with zero-th 
power of A, thus c 7+ i = c 7 . (For example, if V(x) = \ + a>%, 



1 



j x i i g \{x A) 2 ^ x _ ^7 1 _|_ q . ^7 2 _|_ . dangerous terms with a 
simply do not show up in the expression for ip T ) 

After this comment we can come back to evaluation of ((4.126 ). The 



product in denominator, which arised from the Van-der-Monde determinant, 
is already proportional to A n : II 7 =i(A — A 7 ) = A n (1 + 0(A -1 )). Because of 
this and the asymptotic formulas ( (4.127| ), it is clear that if determinant in the 
numerator of ( (4.126Q is rewritten as linear combination of n x n determinants 



with the coefficients y 7 (A) from the last row, only items with 7 > n can 
contribute. There are two such items: 7 = n and 7 = n + 1. In the 
expansion of (n+1) x (n+1) determinant (p n+ i(X) is multiplied by det <p 7 (\s), 
which exactly cancels with determinant in denominator, and the relevant 
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contribution is 



Res rr» 7\ = Cn+1 + 2^ A - 

ll 7 =ll A - ^7 ) 7 



= c n+1 +trA. (4.129) 

The item with <p n (\) is 

det (yi(A 7 ) . . .y? w _i(A 7 ) ^"+i( A 7)) Rc „ ^n(A)ciA 
det (<^(A 7 ) . . . y>„_i(A 7 ) <^n(A 7 )) 68 n 7= i(A - A 7 ) • 1 ' j 

The remaining residue is just unity. Determinant in the numerator differs 
from the one in denominator by substitution of the colomn with entries 
<£> n (A 7 ) for that with <£> n+1 (A 7 ). 

At last we can return to eq. ( |4.122j ) and recall that J^ 7 (7) = <£> 7+1 (7), thus 

d 

logdet 7 5^ 5 (/ 7 ) = 
P (4.131) 

det (^l(/ 7 ) . . . 0n-l(ly) 0n+l{^)) 
det ((^l(/ 7 ) • • • <^ n _i(/ 7 ) 0n(lj)) 

what is just the same as ( |4.13(J[ ) , since (f>s differ from <pg by (^-independent 
factor of e v r ( x )~ xv 'W y V" '(A) . Thus we conclude that the r.h.s. of (refsteqend) 
is equal to — c n+ \, which actually vanishes, as was explained several lines 
above. 

Two things deserve paying attention in this derivation. First, it was ab- 
solutely crucial that we had g|r log Z v at the r.h.s. of ( |4.122| ) to make it 
vanishing, and therefore immedeately appears in the expression for the 

C_2 operator at the l.h.s. (this is the origin of ^-corrections in ( |4.124j ). Sec- 
ond, the result is both simple and natural, but the proof is full of technical 
details and looks somewhat artificial. It becomes even more involved, when 
the general formula ([4 . 1 3 6|) for T^-derivatives of Zy with 1 < k < p 



is 



derived, which plays an important role in the theory of GKM and its appli- 
cations to the theory of quantum gravity. The proof of the string equation is 
just a particular case of that formula, since using integral representation of 
ip(l) one can represent the r.h.s. of ( |4.131|) as -^(trX), where ( ) now stand 
for the average, defined by Kontsevich integral. Thus 

c on Zv (4-izo d Zv + (trA _ trf) (S Q (4 132) 

Oil 
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4.9 On the theory of GKM 

We remind that GKM is abbreviation for the Generalized Kontsevich model, 
This theory is the naturally broad collection of topics for a separate big 
section in these notes. However, we decided not to include such detailed 
presentation. This is because GKM theory seems too incomplete now. First, 
we believe that the natural invariant formulation - of which existing matrix 
integral is only a specific realization - is still lacking. Second, GKM is not 
yet generalized enough to fulfil its main purpose of incorporating infromation 
about all the models of 2d gravity (in fact it should include even more: the 
entire theory of integrable hierarchies and geometrical quantization). Third, 
though the whole approach is very conceptual and deep, many proofs, as avail- 
able nowadays, are still very technical and long. All this implies, that the 
proper view on the subject of GKM still needs to be found. At the moment 
we could describe two complementary approaches: one, starting from inte- 
gral representations, another - from the Duistermaat-Heckman (localization) 
theory and Fourier analysis on group manifolds. Though intimately related, 
these two approaches are still technically different in too many respects. The 
second one is more fundamental (since ordinary integrals arise from discrete 
sums either in special limits or in the cases of inifinite-dimensional algebras, 
and, more important, since integral representation is only one of many pos- 
sible ways to define the quantities of interest). However, many of the most 
important results, obtained in the first approach yet do not have their proper 
names and exact counterparts in the second one. We believe that this whole 
issue will be very much clarified in the near future and decided to postpone 
a detailed review till that time. What we can not avoid in these notes, is 
giving at least a list of topics, already included in the theory of GKM, and 
this is the purpose of the present subsection. 

Kontsevich model with V — was derived by Maxim Kontsevich \T2 
from the original definition of topological 2d gravity, given by EWitten [|9] 
in terms of generating functional for Chern classes of certain bundles over 
Riemann surfaces. Generalization of this reasoning (when more bundles are 
taken into consideration) leads to the theory of Landau- Ginzburg Gravity 
(LGG), which is believed to be the same as GKM, though not all the proofs 
are already available. 

38 Intermidiate results include the study of spherical approximation to LGG, which 
exhibits the structures, peculiar for " quasiclassical integrable hierarchies" (of which Bate- 
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The crucial feature of non-perturbative partition functions, as we dis- 
cussed at the neginning of section 2, is their intrinsic integrability. For 
2d gravity this general idea acquires a very concrete formulation: partition 
functions are usually just r-functions of conventional integrable hierarchies, 
moreover - for LGG, associated with minimal models, these are just ordinary 
multicomponent Toda hierarchies^ 3 ! 

M.Kontsevich found representation for generating functional in the form 
of matrix integral, i.e. formulated a matrix model, which later allowed to 
prove Witten's conjecture that the functional is in fact a r-function. The 
concept of GKM as a universal matrix model, including all the information 
about generic (eigenvalue?) matrix models and thus all the models of 2d(7) 
gravity was introduced in p0|| , and the analogue of Kontsevich model with 
arbitrary potential V(X), i.e. expression 

z v {T}\ Tk=brA . k = c v {k)- l F v (y\k)) ~ 

K k 

ietV 'W r dXe - t , vm+t ,v WX (4133) 



(2vr) n2 /2 e tr(AV"(A)-V(A)) J n 

was proposed as an intemediate step in this direction.^ This (still restricted) 
version of GKM is already enough to unify all the (p, l)-models of 2d grav- 
ity. In some sense, (p, g)-models with q ^ 1 are also included, but in a very 
non-transparent way (using analytical continuation), which does not even 
explicitly respect the p q symmetry. Partition function of such GKM, 
Z V {T} depends on two types of variables: time-variables T k and potential 
V. Formally these two types of variables are absolutely different, V being 
responsible for the choice of particular LGG model or, what is essentially the 
same, of particular reduction of Toda-lattice or KP hierarchy; while are 
parameters of the generating functional of all correlation functions in this 

man hierarchy to be briefly mentioned in section 5.2 below is an example), and which 
also arise in " quasiclassical approximation" to GKM. For some results in this direction see 
& 0,0,0, and references therein. 

39 One can say that this is natural: both such models and Toda hierarchies are associated 
with the level k = 1 Kac-Moody algebras and corresponding simplified versions of the 
WZNW model. However, too much still remains to be clarified about this "obvious" 
connection. 

40 We remind that detV"(A) is defined somewhat tricky, see s.2.5 above. The same 



matrix integral (4.133) was also considered in refs.M,|3ll,[p2| 
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particular model. But of course, since we deal with exact (non-perturbative) 
approach, there is almost no real difference between these types of depen- 
dencies - on the model (vacuum state) and on the T's: the model can be 
changed by non-infinitesimal shift of T-variables. Technically, in GKM this 



is reflected in the identity of the form |40 



Zv p { T } = f P ( r I T k + r k ) ■ r{f k + r k \ G p } (4.134) 

where r k = fc ^_ fc s Res (V (//)) p d\i provide a specific parametrization of 
potentials V, which is here assumed to be any polinomial of degree p, and f p 
is some simple function: 

f P (r | T fc + r k ) = exp -]- ^ A^ifi + n)(fj + rj); 

fs2 1 (p) 

Aj = Res (V'^)Y /p d(V'^)y /p = , (4.135) 

and Tq P ^ is a r-function of " quasiclassical hierarchy". What is important, 
G p (which define the shape of the r-function as function of T + r) and f p 
depend only on the degree p, but not on the other details of the shape of the 
potential. This is a deep formula. It accounts for two phenomena at once: 
First, it says that Z depends on the sum of T and r.Q Second, dependence 
on V is not quite smooth: when the degree of potential changes, the shapes of 
the functions / and r also changes abruptly. Another side of the same phe- 
nomenon is that partition function Zy{T}, which in principle is well defined 
as a matrix integral for all choices of V and L (thus T) at once, is in fact 
singular at some points: there are phase transitions, manifesting the switch 

41 In Miwa parametrization T k — jtr (V^'(A)) k ^ p . Throughout these notes we used 
different time- variables — ■^trA _fc , which are independent of the potential V, instead 
V-dependence of Z v - which we did not really study - was rather nontrivial. If expressed 
in terms of T, partition function Zy\T + r} = Zy{T} becomes almost independent of 
V: it changes - abruptly - only when the degree p of potential changes. This second type 
of description is of course in better accordance with the symmetries of particular model 
- which arc different in different "vacua" (for different p). Therefore these are variables 
T + r, rather than in T which arise naturally in Ward identites - as we saw in sections 
2.5 and 2.6. T's and T's are suited for different purposes: T's are nice when universality 
aspects of GKM are concerned, while T's arise when specific features of particular models 
(orbits, vacua) are considered. 
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from one LGG model to another. After a phase transition original integral 
expression becomes sonewhat symbolical: it defines partition function only 
in the sense of analytical continuation, and it is a separate problem to find 
an integral representation, adequate in the new phases. In practice, what is 
nicely decribed by the integral representation GKM in the form of eq. (|4.133|) , 
are (p, l)-models, with p + 1 being just the power of potential V(x). What 
is not yet found, is analogous representation for (p, g)-models with q ^ 1 (it 
can involve multiple matrix integrals, and universal model is supposed to be 
"matrix quantum mechanics in external fields"). 

Derivation of the crucial formula ( }4.134| )in any approach - starting from 



GKM in the form of either LGG or matrix integrals - is still very tedious. In 
matrix-model representation it relies upon identity |41J 



<)Z[ (trA fc - trX k ) = 

(4.136) 



dT k 



Cy 1 J (trA k - trX fc ) e ^v p (x)+trv;(A)x dx for 1 < k < 



P, 



which look trivial but are rather hard to derive. (A proof of the string equa- 
tion in GKM at the end of the previous subsection is the simplest example of 
this kind of exersices.) Certainly some simple derivation "in two lines" should 
exist, but it is not yet found. Formulas of this kind are very important for 
all aspects of GKM theory. Besides other things, they are just necessary to 
actually evaluate correlation functions in (p, l)-models of 2d gravity, of which 
Zy{T} is a generating functional. If instead of these "physical" questions, 
one asks about integrability theory, identites of this sort also play important 
role. For example, looking at ( |4.136| ) for a special k = p and special choice of 
potential - monomial V P {X) = XP +1 ,- ane can note that the r.h.s. vanishes: 
this is just a Ward identity, reflecting invariance under the shift of the inte- 
gration variable, 5X = const. This is a simplest version of a more general 
statement :0 

if VJX) = , then — - = for all n G Z + . (4.137) 

V + 1 oTpk 



42 This property was technically implicit in original Kontsevich's work p2[ for p+1 = 3, 
where it was related with certain combinatorial identities. A tricky proof, relying upon 
properties of r-functions, was given for any p in |3C[| . An example of straightforward proof, 
again for p+ 1 = 3, - just in terms of Kontsevich matrix integrals - can be found in ref. |99) . 
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Looking from the point of view of integrable hierarchies, one immediately 
recognizes ( 4.137 ) as an example of reduction condition ( [4.40 ). It corresponds 



to the so called p-reduction of KP-hierarhy, of which KdV (p = 2) and 
Boussinesq (p = 3) are the most celebrated examples. We refer to [^0| and 
(|(J for all details and references, the only thing to mention here is that the 



slightly weaker verison of the constraint ( |4.137| ) 

dZ v 



BT 



const, (4.138) 



where a n do not depend on any time variables, can be simply expressed in 
Miwa parametrization: it is just the statement that (^-functions in 

satisfy the p-reduction condition: 

aW) = EVj(a)- ( 4 - 139 ) 

5=1 

This a restrictive relation, because <£>'s are infinite seria in 1/A, while at th 
r.h.s. in ( |4.139| ) there is only a finite number of items. In GKM it is satisfied 
for monomial potential just as a corollary of Gross-Newman equation, or 
more exact, of the Ward identity for the integral 



Indeed, integral does not change under the shift Sx = const, and this implies: 
' x<- x (v\x) - V'(A) - e~ v ^ +v '^dx = 0, 



i.e. 



P+i 



£ kv k (y 7 +*-i(A) - A fc -V 7 (^)) - (7 - l)^ 7 -i = 0. (4.140) 



fe=i 



If only v p+ \ 7^ 0, this leads to an identity of the required form of ( |4.139| ). 



This description of reduction can be modified to allow for non-monomial po- 
tentials, making use of the concept of " equivalent hierarchies", see ||100|| , PC| : 



111 



in this framework the reduction condition is 



\/'(A)y? 7 (A)=£V 7 ^,(A), (4.141) 

8 

but classes of essentially different reductions are labeled by the degree of 
potential only. 

As we already discussed in the previous subsection, linear constraints like 
( |4.139| ) are not restrictive enough to fix the shape of the r-function (the point 
G in the universal module space) unambigously: string equation should be 
also imposed. If expressed in terms of (p's, string equation is just the property 

(p J+1 = Anpy, (4.142) 
where the Kac-Schwarz operator 

A- 1 9 1 V "'W I A f4143l 

It has obvious generalization of the form 

A = _» 1 | rym (4 144) 

where Q q (X) is a polinomial of degree q + 1 and ( f4.142|) is substituted be 

ip 1+q = Ap !q (fy. (4.145) 

This generalization is naturally related to the string equation in (p.q) models, 
see [^TJ and references therein. Generic (p, q) LGG model can be described 
by a system of constraints, 

(A p -V p )M = 0, 

A p , q { V } = 0, (4.146) 

where both operators V p and A Ptq are not uniquely fixed by choosing p and 
q, and there is also a freedom to change variables A — > /(A) and make 
a triangular transformation of basis tp 1 — > (p 1 + X^ <7 C^$(ps- Altogether 
the set of equations ( |4. 146| ) modulo these allowed transformations is finite - 
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(p — l)(q — 1) - dimensional, this is dimension of the module space of LGG 
models with given p and q. Kontsevich integral can be now used to establish 
duality transformation from (p, q) to (q, p) model ET[ : 



Here 



Z V>Q (A) = C^ Q (A) / dXe- tvSv ^ x ^ +tlV '^ Q '^Z Q y(X). (4.147) 

Jnxn 

S v , Q (x, A) = V{y)Q[\y)dy = f V\y)dQ'{y), (4.148) 



As usual, CV,q(A) is the quasiclassical approximation to the integral, and 



Zy ; g(A) = dct7 w A 7 | A ^ , where </? are solutions to ( [4.146|) with V p and ^4 P)(} 



A (A) 

defined by eqs.( |4.l4T| ) and ( |4. 144| ) respectively^ This relation does not 
provide any formula for Z Vpj Q q (A) unless q — 1. The case of q — 1 is 
distinguished because Z Ql y p is trivial. Indeed, the 1-reduction constraint 
A<y? 7 = y? 7+ i + E<5< 7 V y 6<Ps, implies that det^y^A^) = A(A) Yls^i(h), thus 
Zqi,v p — ex vJ2k a kTk, what is essentially the same as Zg 1 y p = 1,0 and 
( f4. 147| ) is just our old formula ( |4.133| ) for the (p, 1) version of GKM. (In 



fact Qi(X) ~ X 2 , and Zg 1 y p is nothing but Gaussian Kontsevich model. It 
is trivial when the "zero-time" N = 0, as we assume here.) Matrix model 
realization of Zv p ,Q q for q 7^ 1 is yet unknown. 

This is not the only important further generalization of GKM ( |4.133| ). 



Another one is implied by the formula for Ty in terms of eigenvalues from 
the section 3.3, 

n „ 

Fv ~ II / dx y e- v(x -' ) A 2 {x)I{x,l). (4.149) 

7=1 J 



43 Also expression for rfc-variables is now modified: 

1. TWT/'C/^-fc/JViT , 



For monomial V p and Q q r k = ~^S k . p+q . 

44 Since if! (A) = f + Ek>o b k^ k , logVi(A) = E k >of X ~ k > and the sum 
X^log^As) = J2k>o 0>2sK k ) = Y.k>o a k T k- Addition of any /mear combination 
of time-variables to log r does not essentially change r-function. For example, the ordi- 
nary integrable equations (like KdV or KP) are usually written in terms of variables like 



u 



-gjr logr, which are second derivatives of logT. 
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As it was already mentioned in section 3.3, the Itzykson-Zuber integral 



[ ^~J^ axu,l ~wS) (4150) 

is in fact a coadjoint orbit integral and has group theoretical interpretation: 
under certian conditions it turns into a character Xr(q) = Tr#g of the group 
GL(n). Here g = e L is considered as a group element, representation R is 
labeled by integer- valued parameters m 1; . . . ,m n - essentially the lengths of 
rows in the Young diagramm. Exact statement is: 

J(ro , .M = ^L = *M, (4.151) 
A(<j) A(m)A(j) <* R 

i.e. in order to get a character we should integrate over matrices X with 
integer-valued eigenvalues.^ Dimension <Ir of representation can be also 
expressed in terms of m- variables: da = A(m). As to the traces trX k = 
J2 7 x t — > J2^ m ^ which appear in the action of GKM, they are very similar 
to the k-th Casimir eigenvalue Ck(R) (though is not exactly the same). Thus 
we see that the integral in ( f4.149|) is in fact very similar to 

H U {9,g} = ExR(9)XR(9)e-^«=° VkCkW , (4.152) 

R 

evaluated at the point g — I. The only real difference is that instead of 
the integral we have a sum over discrete values of m (sum over all the rep- 
resentations, or a model of GL(n)). This " discretized" (quantum?) GKM 
is more general that the continuum one which can be obtained by various 
limiting procedures. It is now obvious that the theory of discretized GKM 
largely overlaps with that of 2d Yang-Mills theory. The simplest ingredient 

45 The ratio 

A(Q = yr 1^ -l s 

A(q) H e l ~i — e l * 

is the usual correction factor, which is the price for the possibility to reduce quantum- 
mechanical problem of motion on the orbit to a single matrix integral. The full problem 
of matrix quantum mechanics can and should be considered as a multi-matrix (in fact, 
infinite-matrix) generalization of GKM ( 4.133| ), which incorporates all the (p,q) LGG 
models. 
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of this theory is the classical result | 101|| , that GL(N) characters are in fact 



(singular) Toda-lattice and KP r-functions. Moreover, the entire sum at the 
r.h.s. of ( ^4.152| ) , if considered as a function of T k = ^trg k , T k = ^trg k is 
in fact a Toda-lattice r-function. There are also features parallel to (|4.134| ). 
We refer to 
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for a little more details about discretized GKM (see also a 
recent paper This is one more very important direction of the further 

investigation of GKM. 



4.10 1-Matrix model versus Toda-chain hierarchy 

At the end of this section we use an explicit example of dicrete 1-matrix model 
to illustrate, how a more familiar Lax description of integrable hierarchies 
arises from determinant formulas. This example will be also usefull in Section 
5.3 below, when one of the ways to take double-scaling continuum limit of the 
1-matrix model will be discussed. Lax representation appears usually after 
some coordinate system is chosen in the Grassmannian. In the example which 
we are now considering this system is introduced by the use of orthogonal 
polinomials. 

We already know from section 3.6, that partition function of 1-matrix 
model (which is a one-component model) is given by 

jV-1 N-l 

Z N = Deto^x^ | V) =n^ = ^ill Rf~\ (4.153) 

i=0 i=l 

where the last two representations are in terms of the norms of orthogonal 
polinomials 

(Qn I Qm) = e+»6 nm (4.154) 

and parameter of the 3-term relation 

hQ n (h) = Q n+1 (h) + c n Q n (h) + RnQn^h), 
Z x = = (1 | 1), R n = e K - K '\ 

Of course all the information is contained in the determinant formula together 
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with the rule, which defines time-dependence of TC{j = (h % \ h?) = TC{ + j: 

ml f f 

o, — rL i+k,j ~ rL i,j+ki UI 

k (4.155) 
dt k ~ n ^ 

(The possibility to express everything in terms of 7i{ with a single matrix 
index i is the feature of Toda-chain reduction of generic Toda-lattice hierar- 
chy.) 

However, in order to reveal the standard Lax representation we need 
to go into somewhat more involved considerations. Namely, we consider 
representation of two operators in the basis of orthogonal polinomials. First, 

n+k I | Lfc I \ n+k 
h k Qn(h) = £ V ' } Q m (h) = E l { nlQm{h) (4.156) 
m^O ( m I m ) rn^O 

(here the simplified notation is introduced for (n \ f(h) \ m) = (Q n \ f(h) \ 
Q m ) and 7^ = — — ; — — .) Second, 



° lk ™=0 \ m I m l m=0 (4.157) 

d(p n _ (n | h k I n) _ {k) 
dt k ~ (n\n) ~ lnn ' 

(These last relations arise from differentiation of orthogonality condition 



0{Q n I Qm) 



O , -Tim r\. 

ot k dt k 

= I Qm) + (Qn | ^) + (Qn I h" \ Q m ) 

otk otk 

by looking at the cases of m < n and m = n respectively.) 

From these relations one immediately derives the Lax-like formula: 

n— 1 n+k 

E 7#7ffi+ E (4-158) 



Z=m-fc Z=m+1 
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or, in a matrix form, 



()f [Rl (q) ,l { % (4.159) 



where 



(We remind that usually i?-matrix acts on a function f(h) = J2n=~oo fnh n 
according to the rule: Rf(h) = J2 n >i fnh n — J2 n <i fnh n with some "level" 
I.) These 7^ are not symmetric matrices, but one can also rewrite all the 
formulas above in terms of symmetric ones: 

£ W s eW^jM = j m j hk ] n) (4.161) 

J (m | m) {n \ n) 

From eqs. (|4.158|) one can easily deduce Toda-equations for <p n : 

d 2 6 n 



E-E 



d (n 


h l n) 


h l n) 


dt k (n 
(n \ h k \ 


n) 
m) (m 


(m 


m){n 


n) 



(4.162) 



\m>n m<nj 

where the r.h.s. can be expressed in terms of R m = e^" 1- *™- 1 . In particular, 
^-S- = Rn+i ~Rn = e^~^ - e^-^. (4.163) 

Let us also mention that in this formalism the Ward identities (Virasoro 
constraints) follow essentially from the relation 

where Hermitean conjugation is w.r.to the scalar product ( | ). For example, 
this relation implies, that 

(Qn | ^) = "(^ | Qn) - E MkiQn \ h^ 1 \ Q n ) . (4.165) 



fc>0 
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Now we note that -Jk is a polinomial of degree n — 1, thus (Qn | ) = 0. 
(In fact 

Also we recall that (Q n | | (5„) = (Qn I Qn) > an d obtain: 

£^-|^ = (4.166) 

fe>o ^fc-i 

for any n. This should be supplemented by relation |^ = 0„. In order to get 
the lowest Virasoro constraint (string equation), L_ 1 Z N = or L_i logZ/v = 
it is enough just to sum over n from to N — 1. 

For more details about 1-matrix model, Toda-chain hierarchy and appli- 
cation of the formalism of orthogonal polinomials in this context see [p6fl . 



5 Continuum limits of discrete matrix mod- 
els 

5.1 What is continuum limit 

Continuum limit of matrix models is, of course, the crucial issue for their 
physical applications whenever these models are interpreted as discrete (lat- 
tice) approximations to continuum theory. The very first thing to be kept 
in mind is that it is not the only possible view on matrix models. Another 
approach considers them as describing topological (and thus also in a certain 
sense "discrete") properties of the theory. Such models, when appearing in 
the field of, say, quantum gravity (which after all is a sort of a pure topological 
theory) do not require any continuum limit to be taken: their discrete nature 
(occurence of integer-valued matrix indices) reflects not the discrete approx- 
imation to the space-time (which does not really exist in quantum gravity), 
but rather the essential discreteness of the underlying structures: topology of 
the module spaces of geometries. Example of matrix models which allow for 
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this kind of interpretation - in terms of topology of module spaces of bundles 
over Riemann surfaces - is provided by Kontsevich models, and this is why 
they usually do not require any continuum limit and why we once called them 
"continuous matrix models" in the Introduction to these notes. The models 
which are usually interpreted in more traditional way - as lattice theories 
- are represented by our "discrete" models, the 1-matrix, conventional and 
"conformal" multimatrix models being included into this class. More sophis- 
ticated examples are provided by "c = l"-theories, Kazakov-Migdal model 
and, say, Wilson's QCD (and infmetely many other lattice theories). It is not 
a surprise that continuum limits of some discrete models provide the theories 
of Kontsevich type: this happens whenever continuum theory is supposed to 
have a kind of topological nature. This is usually the case for quantum grav- 
ity (which, as we said, is conceptually a topological theory in the "module 
space of geometries" - the notion which is already made more or less explicit 
in the 2d case) , but in principle this can be also true for many other theories, 
including exhaustive quantum theory of Yang- Mills fields (again there is al- 
ready considerable progress in this direction, as soon as 2d Yang- Mills model 
is concerned). There should not be confusion about the presence of gauge 
particles in dimensions greater than 2 (for Yang-Mills) and 3 (for gravity): 
there is no reason to prevent generic topological theory from possessing con- 
tinuum spectrum of excitations, though explicit analogue of Kontsevich-like 
description of such situations is not yet found (as we mentioned many times, 
it should probably rely upon non-eigenvalue models). 

We shall not discuss the non-trivial history of invention and understand- 
ing of all these notions (the crucial steps being discovery of the "multiscaling 
continuum limits" |19|],[|(J, which preserve integrable structure of discrete 
models in continuum case; hypothesis of equivalence of quantum and topo- 
logical 2d gravities || and its proof [^3| , |2|] , provided by discovery of Kont- 
sevich models as a peculiar and powerful tool for description of topology 
of the module spaces). Instead, following the mean line of these notes, we 
shall concentrate on intrinsic relation between (multiscaling) continuum lim- 
its and integrability: the notion of continuum limits is in fact built into the 
theory of integrable hierarchies and the underlying representation theory of 
Kac-Moody algebras. 

In the case of the eigenvalue models the central issue here is the interrela- 
tion between Toda-lattice and KP hierarchies, even its more narrow aspect: 
elimination of the zero-time N, present in the Toda-lattice case. From rep- 



119 



resentation theory (or conformal field theory, what is essentially the same) 
point of view the thing is that the zero-time (which labels the filling level 
of Dirac sea in the fermionic picture) is associated with the zero-modes of 
scalar field and its elimination is just the change of boundary conditions 
which eliminates zero-modes. The simplest example of this "twisting" pro- 
cedure is just transformation from periodic to antiperiodic scalars - it still 
preserves possibility to have fermionic description (where it looks like a switch 
from Ramond to Neveu-Schwarz sector), and thus does not take us out of 
the field of conventional integrable hierarchies. In representation theory one 
can interpret the same operation just as a switch from the homogeneous to 
principal representation, which are associated with the Toda-lattice and KP 
hierarchies respectively. 

This remarkably simple description is of course far from obvious, if one 
investigates continuum limit in naive way, without taking integrable structure 
into account explicitly, but just sending the number of degrees of freedom in 
discrete theory (i.e. the matrix size N) to infinity (together with the inverse 



lattice spacing, if any). We refer to the classical review [|T^] for discussion 
of what are the naive continuum limits in lattice gauge theories, i.e. what 
are the conditions for getting the second-order phase transitions, which allow 
for a continuum-like scaling behavious in the vicinity of the critical point, 
with critical exponents defining all the continuum physics, from quantum 
dimesnion of the space-time to spectrum of particles. The problem with 
naive continuum limits is that they can easily destroy integrable structure of 
the theory (the underlying hidden symmetries), unless special precaution is 
taken: the critical point (which is in fact a low-codimensional hypersurface 
in the infinite-dimensional space of parameters) should be approached from 
the certain directions, so that Ward identities are not explicitly broken. 

As soon as this word - Ward -dentities - is pronounced, we already get into 
the field of integrable systems and the issue can be discussed inside this field. 
The above-mentioned switch from periodic to antiperiodic fields is of course 
apparent if the discrete and continuos Virasoro constraints (represented by 
formulas (1.2) and (1.3) in the Introduction) are compared, but this is a 
posteriori information, because so far we interpreted "continuous Virasoro 
constraints" as the Ward-identities for the V = X 3 -Kontsevich model, and 
it still remains to be explained why Kontsevich model is indeed what arises 
after continuum limit is taken. The simplest approach to this problem is 
to make use of the identity between discrete 1-matrix model and Gaussian 
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Kontsevich model | 5"6fl , established in section 3.8 above. Then the X -model 
arises in the large- N limit just when the matrix integral is evaluated by the 
steepest descent method |3TSfl . We shall present this simple calculation in 
the last subsection below, but before we take a somewhat more direct (and 
complicated) approach in order to reveal at least some of ideas, underlying 
the entire theory of continuum limits. 

5.2 From Toda-chain to KdV 

We begin with the simplest existing example: continuum limit, in which the 
lowest equation of the " Volterra hierarchy" , 

dR n 



Rn{Rn+l ~ Rn-1 ), (5-1) 



dt 

turns into the lowest KdV equation: 

^- = -~r"' - 2rr'. (5.2) 

Volterra hierarchy is a reduction of Toda-chain hierarchy, with R n = e^™ - ^™ -1 , 
arising when all the odd-times t 2 k+i = and all 4> n are supposed to be inde- 
pendent of them. Therefore this hierarchy is clearly related to the discrete 
1-matrix model. We" 11 turn to the study of 1-matrix model in the next 
subsection, but here we just address the transformation from (|5.1|) to (|5.2j) 



103 , 26 



The basic idea of taking continuum limit is to change discrete "zero-time" 
n for continuum variable x (to be after all substituted by T\ of the continuous 
hierarchy). In other words, the idea is to consider a subset of functions R n , 
which satisfy Volterra equation and depend on n very smoothly, so that 
they can actually be substituted by a smooth function R(x). This is a very 
natural thing to do, of course, when one is interested in the large-n limit of 
the equation. Namely, one substitutes (|5.1j) by 

= -R(x)(R(x + e) - R(x - e)), (5.3) 

and take the limit e — > 0, which, after rescaling x — > ex, gives rise to "Bate- 
man equation" , 

= -R(x)R'(x). (5.4) 
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This is a very interesting equation (see ||104|| for description of the amusements 
of the related theory, which is in fact intimately related to the theory of 
jets). However, it is much simpler that KdV equation (for example, it is 
completely integrable in the most trivial sense of the word: entire set of 
solutions, satisfying any boundary conditions can be immediately written 
down, see ||104|| ). KdV equation can be considered as a sort of "quantization" 
of (|5.4|) (unfortunately this very interesting subject did not yet attract enough 
attention and is not studied well enough). 

Remarkably, Bateman equation is not the only possible limit of Volterra 
equation: a fine tuning procedure ("double-scaling limit") exist, which can 
provide less trivial - KdV - equation ||103 |. Indeed, imagine, that in contin- 
uum limit R n tends to a constant Rq, and the function r(x) arises only as 
scaling approximation to this constant: R(x) = -Ro(l + e s r(x)). Then the 
leading term at the r.h.s. of fl5"H| ) is eRR'(x) = —2e s r(x)(l + 0(e 2 1 e s )) ) and 
instead of ( [5 . 4| j we would get: 

- = -2eR r'(x)((l + O(e 2 ,e s )). 

This equation is even simpler that ( |5.4|) - it is just linear, but in fact it is too 
simple to preserve its form: by a simple change of variables^ 

x = x - 2eR t, (5.6) 
t = e 3 R t (5.7) 

it can be transformed into 

and terms at the r.h.s. also deserves beeing taken into account. Then we 
get: 

dr(x) 



dt 



-2eR ( 1 + e s r(x))(r'(x) + -e 2 r"'(x) + £>(e 4 ) 
V 6 



■2e.Ro r'(x) + -e 2 r"'(x) + e s rr'(x) + e 2 0(e 2 , e s 
6 



4(i 



This change of variables is implied by the relation: 

d , d [dt , di\ d fdx , d£\ d d , . 
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and, after the change of variables (|57 



^1 = -- r '"(x) - 2e s - 2 rr'(i) + 0{e\e s ). 
at 3 

It is now clear, that the choice s = 2 is distinguished (a critical point) and 
at this point we get: 

dr 1 d 3 r _ dr . o . 

2r— — , (5.8) 



dT 3 3dTf dTi 

where new notation 7\ and T 3 is introduced for x and t respectively. This is 
already the KdV equation (|5.2j), and our conclusion is: 

While the naive continuum limit of Volterra equation is just a simple 
Bateman equation, the scaling limit can be fine tuned so that KdV equation 
arises instead. The crucial ingredient of this adjustement is the change of 
time- variables {t} — > {T}, which involves singular parameter e. The proce- 
dure can be easily generalized to the entire Volterra hierarchy, and fine tuning 
allows to get the entire KdV hierarchy in the limit of e — > 0. Usually trans- 
formation to the "Kazakov variables" {T} (they are a little different from 
those originally introduced by V.Kazakov in [|19|]) from {t} is some linear 
triangular transformation. 

An important detail is that this procedure requires restriction to only 
even time- variables t2 m , m > 0. (If odd times are also involved, a pair 
of KdV hierarchies arises in the continuum limit - this is not a "minimal" 
case.) Thus "irreducible" realization of continuum limit requires reduction 
of original hierarchy. This can be also seen from the fact that the lowest 
KdV equation arises from the lowest Volterra equation, which is related to 
the second eqaution of Toda-chain hierarchy. 

Unfortunately this simple piece of theory (continuum limits in terms of 
hierarchies) has never been worked out in full details (for the entire Toda- 
lattice hierarchy, its multicomponent generalizations and their reductions). 
As we already mentioned, this theory will involve the general relation between 
homogeneous and principal representations of the (level k = 1) Kac-Moody 
algebras. 

5.3 Double-scaling limit of 1-matrix model 

Now we proceed to discussion of a slightly different approach to continuum 
limits, which is directly adjusted to the needs of matrix models. The naive 
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idea [pm,P9| is to forget about integrability and just look at the Ward identi- 
ties (Virasoro constraints in the 1-matrix case) and take a continuum limit of 
these identities. This approach makes close contact with the standard tech- 
nique of "loop equations" (Makeenko-Migdal equations ||105|| ) in the theory 



of matrix models, of which Virasoro and jy-constraints are just particular 
examples.^ 

However, carefull analysis of continuum limit of discrete Virasoro con- 
straints |28] makes it clear that the procedure is far less simple than one can 
think in advance (usually derivations are not very carefull and details are 
just "put under the carpet"). The crucial problem is that we want peculiar 
(double scaling) rather than naive limit, and, as we mentioned in the pre- 
vious subsection, this also requires a certain reduction (elimination of the 
odd-times t 2m +i)- If parity symmetry (w.r.to the change of H — > —H in the 
original matrix integral) is taken into account, one can easily throw away 

dZ 

first derivatives w.r.to the odd-times t 2m+ i, just because — 

Ob 

but this is no longer true as far as the second derivatives 



2m+l 

N 



d 2 Z, 



= 0, 

2fc+l=0 



dt 2m+l dt 2 i- L t2k+1=0 
are concerned, which appear in (the "quantum piece" of) the Virasoro con- 
straints (1.2). It is a highly non-trivial feature of loop equations (having its 
origin in their integrable structure!), that in continuum limit these terms can 
be in fact carefully eliminated. The thing is that the second derivatives of 
log Z N appear to be a local objects, in the sense that they depend only on 
Zfj with the difference | N — N \< m + I, which does not blow up as iV — > oo 

in continuum limit. Moreover, the differences — — — ^ — ^- — %r— al- 

dt 2m+ idt 2 i-i ot 2m dt 2 i 

most tend to zero, leaving some simple (though vitaly important) correction 
to arising continuous loop equations. This locality property allows one to 
get rid of these dangerous odd-time derivatives, substituting them just by 
second derivatives w.r.to the even-times. Since such substitution is possi- 



47 One of the puzzles in the theory of non-eigenvalue models is to identify group- 
theoretical meaning of generic loop equations: they are usually introduced as equations of 
motion rather than as Ward identities (see discussion at the beginning of Section 2 above) , 
and thus their implications are more obscure and technical means to deal with them are 
much more restricted. When group theory description will be found, it will very soon 
reveal the (generalized) integrable structure of non-eigenvalue models and it will be a big 
step forward in the whole theory. 
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ble only for logarithms of Zn, continuous constraints appear imposed on the 
square root of original partition function (or on the i-th power in the case of 
the p — 1-component conformal models). Another aspect of this trick to deal 
with the odd-time derivatives is that it makes the entire derivation depend- 
ing on the fact that the theory is integrable - this is what guarantees the 
above-mentioned locality. Since the way to reveal integrability, by looking 
at the loop equations themselves is yet not very well understood, the whole 
calculation becomes not quite self-contained (but of course, if we know ev- 
erything about integrable structure this is not a real drawback, this is just 
a limitation of particular approach, starting from the loop equations) . In 
particular, this is the only loophole, which is still not filled in the description 
of continuum limit of conformal (multi-component) matrix models, which in 
all other respects goes in exactly in parallel with the 1-component (1-matrix) 
case.0 

We shall now describe briefly the steps of this calculation for the 1-matrix 
model, refering for all the details to refs. [^] and . Our previous discussion 
already contain motivations for the main steps, so we do not need to go into 
detailed explanations. Manipulations below, involving Kazakov variables 
can look a little artificial, but we repeat that they can be interpreted as a 
switch from the Toda-type to KP-type hierarchies, which, as we already saw 
in the previous subsection, is naturally associated with the double-scaling 
continuum limit. 

We start from the discrete Virasoro constraints (1.2), rewritten in terms 
of a generating functional ("stress tensor" on the spectral plane): 

L-(z)Z N = 0, (5.9) 

where 

OO 1 

L_(*)= E L n z- n -* = -(j\z)) , (5.10) 

n>-l A 



48 It transforms discrete W^-constraints into continuum W^-constraints. which in their 



turn arise from the GKM with the appropriate potential |30[] , 55 . Unfortunately, since 
the GKM-inteprpretation of discrete multicomponent models (like the one existing in the 
1-matrix case, see s.3.8) is yet unknown, the direct way to take their continuum limit 
- like the one to be described in the next subsection for the 1-matrix case - is also yet 
unavailable. For more details about conformal matrix models, their integrable structure 
and continuum limits see refs.[09l. 



125 



and 



J(z) = d<j>(z) = £ Jnz" 71 - 1 ; 

n=— oo 



(5.11) 



^2 fc > fc>0 fc dt k 

d 1 

J_ fc = V2— ; J fc = — =Art fc , fe > 0; 
Ota: V2 

d 

— — Zjy — NZjy. 

Next, we need to reduce the original partition function: 

Zn{£} ► Z^ d {t even } = Zjv{todd = 0,t ev cn}- (5-12) 



dZ 

All odd virasoro generators L 2n+ i act trivially on Z™ d , since — 

C^2fc+1 

and we need to consider only L 2n . Introduce also^ 



= 0, 

t o dd=0 



ired / 



U cd (z) = \ (d^ ed (z 



Now we have two issues to be discussed separately. The first one is the 
change from t 2 k to Kazakov variables T2 m +i- The second is the difference 
between constraints imposed on Z red and Z. 



49 Note that </f> red (.z) ^ </>(-z)L _ n and similarly LJjrfC 2 ) 7^ ^2«L _ n : some factors of 

'''odd — " ' ''odd — " 

2 in ( |5.13| ) being responsible for this discreapancy. In fact L lcd are related to generators 
of the Virasoro constraints in the complex-maiiix model p8j , 



Z% = [ dMexp [ ^t2fcTr(MMt) 

•* \k>0 



and in continuum limit Z^ ~ \l %2n- 
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The simplest way to describe Kazakov variables is to introduce one more 
- antiperiodic - scalar field, 

V ^ fc> fc>0 K ~+~ 2 C J 2fe+l 

Here T and T are related by transfromation 

T 2fc+ i = f 2k+1 + e 2 — -f 2k _ x + 2eNS kfi . (5.15) 

fc + 2 

Impose now a relation: 

,90 rcd (2) = l l U- 1 &f>(u)U; 

z 2 = 1 + e 2 M, (5.16) 

and in continuum limit e is assumed to vanish. This is a relation which 
maps homogeneous representations into principal, but its invariant meaning 
(especially - from the point of view of conformal field theory) does not seem to 
be enough understood. Anyhow, this relations establishes a relation between 
t even and T. Namely, comparing the coefficients in front of the positive powers 
of u at both sides of this equation, we get: 

T — l c 2fc+l v^oo g m T(m+\) h ~~> C\- 

±2k+l - 2 t 2^ m >k ( m _ fc ) !r ( fc+ |)> ft ^ u ' 

g m = mt 2 m, m > 1; g = 2N. (5.17) 
Inverse transformation looks like 

a ~ 2 V ( ) k - m T2k + lT ( k + D (Kie) 



Now, 



' lg T(k + ^ 8 (51g) 



9*2* 2 ^ (k - m - l)!r(m + |) OT 2m+1 
and using this formula when comparing the negative powers of u we find: 

U = exp (j2 m ,n A mn T 2m+1 T 2n+ ij , 

_ ( )m +n r(m+§)r(n+|) 9f) ^ 

^ ~~ ^ e 2(m+n+l) m !n!(m+n+l)(m+n+2) • (U.^u; 



127 



The square of relation ( |5.16| ) is: 

.2 _ 1 
e 

or 



recA - -U- 1 (<9$) 2 (u)U, (5.21) 



E <^ 2P " 2 = V 1 ( E ] U. (5.22) 



p>0 e \">-l 



This equality implies that 

u n+l du 



U- 1 C 2n U = e A Y.L^i 



r 2p+2 
p>0 -°° Z 

n+1 

_ — 2n / \n+l— p/^yp r red 
p=0 

since 

M n+1 rfM / M n+1 rfM 1 F(-p) 



(5.23) 



e 4 



oo z 2 p+ 2 loo (1 + c 2 m)p +1 e 2n (n + 1 - p)W(-n - 1) 
H n+P+1 (n + 1)! 



e 2n p\(^ n _|_ _ 



c 2n ^n+l- 



Explicit expressions for the generators (which are harmonics of the stress 
tensor | (<9<l>) 2 (w) of antiperiodic field $(«)), are: 

1 <9 T 2 

£-2 = E(^ + ^+i-~ + f; 

fc>l 2 C-t2(fc-l)+l 4 



1 d 

^0 = E( fc + o) T 2fc+l 



fc>0 Z O r -'2fc+l 

1 d 



(5.24) 



fc>0 ^ ^?2(fe+n)4-l 

1 ^ <9 2 (-)" 

4 fc=0 ^2fc+1^2(n-fe-l)+l n 

So far everything what was done was just change of variables and all relations 
were exact for any e, no limits were taken. 
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Operators ( |5.24j ) are very similar to £,2n, arising in the "continuous Vi- 
rasoro constraints" (1.3), imposed on partition function of X 3 -Kontsevich 
model. There are, however, two discreapancies. 

First, appear in ( E>.24|) instead of ^ in generators in (1.3). One can 

argue that this difference is not really essential, since T 2k+ i and T 2k+ \ differ 
by terms, which are proportional to e 2 and thus vanish in the continuum 
limit e — > 0. (Note, however, that this reasoning can be applied only for 
every particular constraint £2 n Z = 0, n > —1, not to the entire generating 
functional, where different terms are summed, multiplied by different powers 
of e.) 

The second discreapancy is a litlle more serious: it is the occurence of 

(_)n+l 

an extra term — for all n > (this difference is present for n = 

16e 2n - 



as well, because Cq contains the item j^, which is lacking in ( 5.24| ).) This 



extra term can not be eliminated by just taking continuum limit: moreover, 
it blows up instead of vanishing when e —>■ 0. Remarkably enough, this term 
disappears when we turn to consideration of actual Virasoro constraints, not 
just a formal choice of time variables. It cancels completely with the other 
potential source of problem for the derivation of continuous Ward identites. 
We proceed now to this most sophisticated matter in this whole subsection. 

The thing is that, as we already mentioned before, the reduction of dis- 
crete Virasoro constraint L 2 „^tv — contains some non- vanishing terms with 
the odd-time derivatives: 



k>0 Ol2k+2n k =Q Ol2kOl2n-2k 



reel 



pp. n-1 o_ 

° j z ied 

xk=0 dt2kdt2n~2k k=0 dt2k+ldt2n-2k-l 



(5.25) 



We added an extra term with the second even-time derivatives to both sides 
of the identity, in order to get at the r.h.s. a combination, which has a chance 
to vanish in continuum limit. (This formula still needs to be corrected, see 
eq.fTgg) below.) 



In order to find rigorous reason for elimination of the terms at the r.h.s 
we need to address to explicit formulas from the last subsection of section 4 
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(no simpler way is known so far). The crucial formula which we need is: 



d 2 



dt k dti 



d (n 


h l 


n) 


dt k (n 


n) 



\in>n m<n/ 



(n 


h k 


m)(m 


h l 


n) 


(m 


m){n 


n) 



E-E 1 7\ 1 1\ > ( 5 - 26 ) 



and the most important feature of it is its .R-matrix structure (the fact that 
a difference occurs at the r.h.s.). This structure implies almost complete 
cancellation of terms, when we sum over n in order to get log Zn = Y^o" 1 <fin, 
leaving only a finite sum of the length independent of N: 



d 2 \ogZ N 



N-l 

E I E V 'T l'\J J ' v -7 j y\ i "' \- ( 5 - 27 ) 



(n 


h k 


n + j)(n + j 


u 


n) 


(n 


n) (n + j 


n + j) 



dtkdf o<j<min(fc,0 \n=N-j 

The finite sum at the r.h.s. can be expressed in terms of R n = e^" - ^"- 1 , 
which are exactly the quantites to satisfy equations of Volterra hierarchy 
and tending to constant (denoted by Rq in the previous section) in continuum 
limit. Locality property - the finiteness of the sum at the r.h.s in ( |5.26| ) - 
implies that this r.h.s. tend to a constant value as iV — >• oo. This constant 
does not completely cancels in the difference 

/ n 02 n-l a2 \ 

{^wit — ^dr^Ft — r° gz " (5 ' 28) 

and the remaining contributions appears to be exactly what necessary to 

r_\n+i 

cancel the dangerous term — which appeared in the difference between 

6 16e 2n w 



C n and C n . We refer to ref . ||28| for more details about these cancellations, 
and the only thing to discuss at the rest of this subsection is the difference 
between the r.h.s. of ( |5.25 ) and Q5.28Q . In the second expression the second 



derivatives are taken of log Z, while they are of z itself in the first one. Of 



d 2 log Z r £ d 1 d 2 Zf d 



J N 



dt2k+ldt2n-2k-l Z r § A dt2k+ldt2n-2k~l 

but this is not true for even derivatives. So, identity fl5.25|) yet needs to be 
transformed a little more in order to contain exactly (|5.26| ) at its r.h.s. If 
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this is achieved, the l.h.s. acquires additional contribution and turns into 

£ 2kt 2k dz " d + jr (2- d2z * d 1 dz » d dz » d ^ - 



k>0 9t2k+2n k=0 \ dt2kdt2n-2k Z™ d 8t2k C^2n-2fc / (5 29) 

As a result of all this reasoning we conclude that the double scaling con- 
tinuum limit of reduced 1-matrix model can be described by the following 
relation: 



lim ^Z]f{t even } = U~ X Z V _^{T}, (5.30) 

d.s. e— >0, N— >oo ' v ~ - 



3 



where factor U is defined in ( p. 20] ), relation between t and T- variables is 
given by (|5.17|) and Z X 3 {T} is X 3 -Kontsevich model. The motivation for 

V ~ 3 

this conclusion is that both sides of the equation satisfy the same continuous 
Virasoro constraints (1.3). 

This whole derivation can be straightforwardly generalized to the case of 
multiscaling limit in conformal matrix models and analogous relation con- 



tains roots of the p-th degree, see [39| for detailed discussion. 



5.4 From Gaussian to X s Kontsevich model 

We shall now abandon these complicated matters and give a simple illustra- 
tion of how the things can work, if expressed in the adequate terms. Namely, 
as alternative to the sophisticated procedure, involving explicit switch to 
Kazakov variables and the study of limits of Ward-identites (loop equations), 
we shall just use the equivalence of the discrete 1-matrix model and Gaus- 
sian Kontsevich model, proved in the section 3.8 above in order to take the 
continuum limit just of this simplest Kontsevich model. This procedure, sug- 



gested in ref.|36| appears to be just a kind of a standard eveluation of the 
integral in the large iV-limit by the steepest descent method. It is important 
here that GKM is not sensitive to the size of the matrix n in Kontsevich 
integral, therefore this limit, when expressed in terms of GKM, has nothing 
to do with the infinitely large matrices. 
Relation to be proved below is 

, to? ^{V} =?{V}> ( 5 - 31 ) 

d.s. N— >oo 1 ' L 1 
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where V(X) = \X 2 -N\ogX and V(X) = §X 3 . 

Very naively, what happens as N — > oo is that in the Kontsevich integral, 

JdXexpti ^X 2 + XlogX + AX^ (5.32) 

a stationary point arises at X = X , such that 

N 

X = — + A. (5.33) 

Expansion of this action in powers of X = 7 _1 (X — Xo) comes entirely from 
the logarithmic piece: 



s-s„ =£^-jv(io g (i + f)-£ 

4(l + 5)^ + J*»fh£)'- (5-34) 



In the continuum limit 7 should be adjusted in such a way, that quadratic 

/ N \ -1/2 

term is finite, i.e. 7 ~ (^1 + j^rj . Now, if A remains finite as N — > 00, 

X ~ x/iV, 7 ~ 1 and all the terms with k > 3 in the sum are damped as 
-f k NXo k ~ N l ~2. This is the naive continuum limit. However, it is clear, 
that one can usually ask A to behave more adequately - blow up together 
with growth of N - and fine tune the way in which it tends to infinity so that 
at last the first term with k = 3 also survives. For this purpose A and thus 

2 ( N \ N l 3 

Xq should scale in such a way, that both quantites 7 1 + — ^ and „ 

\ X Q J X 

remain finite. This requirement in the case of the latter expression means 

that 7 ~ X N-^ 3 and then 

N + X 2 

AT2/3 ' 

This is never finite, unless N + Xq — > as N — > 00. This in turn implies 
that Xo ~ i\/~N and A — > 2Xq ~ 2iy/N should be pure imaginary. One can 
also check that the terms with k > 3 in the sum ( |5.34| ) all tend to zero in 




this specific limit. Thus we are left with a model which has only cubic and 
quadratic terms in the action. By simple shift of variables quadratic term 
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can be changed for a linear one and we get a description of the theory in the 
vicinity of the stationary point in terms of an X 3 -Kontsevich model. 

In practice things are a little more complicated, because also reduction 
to even-times should be taken into account. However, this does not really 
add too many new problems. We need that only even times t 2 k = ^tr^r 
remain non-vanishing, while all the odd times t 2 k+i = ^FTT^ a 2 *^ 1 = ^' This 
obviously implies that the matrix A should be of block form: 



A = 7 (5.35) 




and, therefore, the matrix integration variable is also naturally decomposed 
into block form: 




X=[ v , ; ] . (5.36) 

Then 



\\ =x 2 /2-N\ og x} - j dXdyd 2 



(5.37) 

det(xy - z-zy) N e^ z \ 2+x2 ' 2+y2 ' 2 - MX+My \ 

To take the limit iV — > oo, one should assume certain scaling behaviour of X, 
y and Z. Moreover, our previous naive consideration gave us some feeling 
of the fine tuned scaling behaviour can look like. So we take 

X = 1 (ipi + x), 

y=^' + y)' (5 . 38) 

Z = lC, 
M = ^ 1 (ial + m) 

with some large real a, (3 and 7. If expressed through these variables, the 
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action becomes: 

tr (\z\ 2 + x 2 /2 + y 2 /2 - mx + My - N\og(xy - z^-zy) 

= 7 2 tr Q(i/3J + x) 2 + ±tr(i/3J - y) 2 + |^| 2 ) - 
— tr(ial + m)(2ij3I + x — y) — 

-Mr log /3 2 7 2 f 1 - ^ + ^ - ^(1 + 0(1/ (3))) = 



(5.39) 



= [2a(3 - (3 2 Y - 2N log /^tr I - 20 tr m+ (A) 
+* [Pi 2 -a + jj(tTX-tiy) + ^ W ~ T^J (tr x 2 + tr y 2 )+ (B) 



+ ( v 7 2 + ^Jtr|C| 2 - (C) 
— tr mx + tr my + — — tr(a; 3 — y 3 ) + (D) 

+0{N/^) + o(\C\ 2 j)j. (E) 

We want to adjust the scaling behaviour of a, (5 and 7 in such a way that 
only the terms in the line (D) survive. This goal is achieved in several steps. 

The line (A) describes normalization of functional integral, it does not 
contain x and y. Thus, it is not of interest for us at the moment. 

Two terms in the line (B) are eliminated by adjustment of a and 7: 

9 N 2N , . 

7 2 = ^,« = T . (5.40) 

As we shall see soon, 7 2 = N/ f3 2 is large in the limit of N — > 00 . Thus, the 
term (C) implies that the fluctuations of £-field are severely suppressed, and 
this is what makes the terms of the second type in the line (E) negligible. 
More general, this is the reason for the integral Z^ v y to split into a product 
of two independent integrals leading to the square of partition function in 
the limit iV — > 00 (this splitting is evident as, if Z can be neglected, the only 

fx Z \ 

mixing term logdet ( ^ j turns into \ogXy = logX + logy). 
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Thus, we remain with a single free parameter (3 which can be adjusted so 
that 

(3 3 

— — > const as N oo 
N 

{i.e. (3 ~ iV 1/3 , 7 2 ~ N 1/3 , a ~ iV 2/3 ), (5.41) 

making the terms in the last line (E) vanishing and the third term in the 
line (D) finite. 

This proves the statement (|5.31|) in a rather straightforward way. Unfor- 
tunately no generalization of this procedure for other discrete models is found 
so far, the main problem beeing identification of GKM-type realizations of 
other (for example, conformal) discrete matrix models. 



6 Conclusion 

We came to the end of our brief review of the facts, that are already known 
about the relation between matrix models and integrable hierarchies. There 
are still several topics, which are already discussed in the literature, but not 
presented in these notes. 

First of all, we did not discuss the relation between matrix models and 
theories of topological (Landau-Ginzburg) gravity (LGG). This field is fastly 
developing during the last months and will be soon ready to inclusion in this 
kind of reviews. The list of things which are already clarified enough, in- 
cludes realization of the Ward identities in the form of "recursion relations" 
for topological gravity |J. Also the relation between quasiclassical hierar- 
chies, arising in the spherical approximation to topological theories [PB| , to 
integrable structure of Generalized Kontsevich model is more or less under- 
stood [ Hfl . Of special importance is the chapter of this theory, which pro- 



vides matrix-model description of module spaces, associated with Riemann 
surfaces [ |106| |,p2fl. What still deserves better understanding is axionatic con- 



struction of topological gravity, similar to remarkably simple construction of 
topological LG models (before they are coupled to 2d gravity) in terms of 
the Grothendieck residues and chiral rings ||107|| : see ||108|| for a very nice 



presentation of the latter case and |l7j for the first big steps towards similar 



construction in the former case. Also relation to the theory of non-conformal 
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LG models [|109|| deserves being clarified. A piece which is essentially lacking 
so far is the clear description of minimal (p, g)-models, coupled to 2d grav- 
ity in the case of p ^ 1. In this situation Generalized Kontsevich model is 
known to describe nothing more but duality transformation between (p, q) 



and (q,p) models j|T|, rather than the models themselves. This subject is 



also connected with the theory of Kac-Schwarz operator ||110|| . The work in 
this direction is extremely important for the understanding of unification of 
various string models and of essential symmetries of the future string field 
theory (in particular generic BRST and Batalin-Vilkovisky symmetries are 
very close analogues of the complete sets of the Ward identites, as described 
in the general framework in the beginning of section 2). All these things 
would constitute a natural next section of these notes, but we choosed to 
wait a little longer untill further clarification is achieved in this fragment of 
the theory. 

Second, we did not touch at all physical interpretations of matrix models, 
which include quantum gravity, Yang-Mills theory and many other possible 
applications. This should be a subject of very different reviews, for which 
the whole content of these notes is just a piece of techniques involved in the 
study of physical phenomena. 

Third, the biggest terra incognita in this branch of science, which re- 
mained beyond the scope of these notes, is the theory of non-eigenvalue ma- 
trix models, which are related to physical theories in space-time dimensions 
d > 2. It is indeed a terra incognita, at least from the point of view of the 
semi-rigorous analysis, which we are reviewing. The recent breakthrough in 
this field is due to appearence of the Kazakov-Migdal model p5| (see also the 
lates review [ |1 1 1|| and refernces therein), which for the first time opened the 
possibility to treat a wide class of non-eigenvalue models by exact methods 
of localization theory (other names for this field, which in fact is growing 
up into generic theory of integability are Duistermaat-Heckman theorem or 
Fourier analysis on group manifolds). The work in this direction is, how- 
ever, only at the early stages and this is why we decided not to present the 
first non-systematized results in these notes. A part of it which is very close 
to beeing satisfactrily understood is the "boudary model" of c = 1 string 
("d = 2 dilaton gravity") - a very important one from the point of view of 
general string theory. For the present state of knowledge about this model 
see [[11211 , and its relation to integrability theory is partly revealed in refs. 
PI and gg. 
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In the domain, which was actually reviewed the weakest points are the 
theory of continuum limits and that of the multicomponent hierarchies. 
These theories, when developed, can also help to move in the most impor- 
tant direction, which was many time mentioned in the text above: towards 
creation of more general theory of integrability. The next natural step, when 
approached from this side should be generalization of conventional integrable 
hierarchies, which would lift the restriction to level k = 1 simply-laced Kac- 
Moody algebras and unitary representations. The emerging theory will of 
course have much to do with both localization theory and non-eigenvalue 
matrix models, and when it is created we shall find ourselves at a new level 
of understanding, which will be one step closer to the goal of construction of 
the entire building of string theory (mathematical physics) and will probably 
provide us with unexpected new means for investigation of the features of 
the real physical world around us. 
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